Oliver J. Woodford
Oliver J. Woodford
I could not find any mention of doing within vector reductions, or equally transposing blocks of vectors to make such reductions vectorizable. If I am summing a long list of...
Basic arithmetic involving arrays and scalars is broken (see #691). The goal of this PR is to fix it. I'm new to Julia, and unfortunately understanding the internals of Symbolics.jl...
The following does not work: ```julia using Symbolics function proj(x) return x[1:end-1] ./ x[end] end @variables x[1:3] s println(proj(x)) # (broadcast(/, x[1:2], Ref(x[3])))[1:2] println(proj(x * s)) # ERROR: MethodError: no...
I'm getting the following test set [error](https://github.com/JuliaSymbolics/Symbolics.jl/runs/7813755491?check_suite_focus=true#step:7:427) in #692: `isequal(s + -s - x[2] + x[1], x[1] - x[2])` returns `false` Clearly these two expressions are equal. Indeed, if I...
If I have a function `f(x::Real)::Real` I have two options for computing `f(x), f'(x), f''(x)`: ```julia using ForwardDiff, DiffResults, StaticArrays, BenchmarkTools function computehessian(f, x::AbstractArray) result = DiffResults.HessianResult(x) result = ForwardDiff.hessian!(result,...
I am seeing some memory corruption in a function that I believe should be working. The code is: ``` using LoopVectorization function At_mul_A!(C, A) @turbo for n in indices((C, A),...
Below I create identical non-Hermitian matrices of types `Matrix` and `SMatrix`, then call `cholesky()` with `check=false`, i.e. don't throw errors. ``` using LinearAlgebra, StaticArrays mat = randn(10, 10) v =...
Sometimes the sparse structure of a problem can change during non-linear optimization. It would be useful to be able to update the symbolic analysis of the matrix in-place, to avoid...
It's great that you've decided to maintain this tool :heart: I am interested in doing all the memory layout randomization from within an LLVM compiled application of my own. I...