Open roflmaostc opened 9 months ago
This error is from Optim's alloc_H
method which relies on being able to transpose the optimization variable. Replacing that with Adam actually works so it's not really an issue but more of a limitation of the methods
res = solve(problem, OptimizationOptimisers.Adam(), maxiters=20)
I'm not sure if I understand but Optim.jl can optimize on 3D (ND even) arrays.
I raised this issue since I always use Optim.jl and it works but it failed with Optimization.jl
Interesting, do you use out-of-place derivative evaluations there?
Usually Optim.optimize(Optim.only_fg!(fg!), x, opt=LBFGS())
but all of it works.
3D optimization with
LBFGS
fails. 2D arrays work because adjoint is probably defined.results in: