Hi! I've got following output: ``` 21.085590 seconds (3.80 M allocations: 105.695 MiB, 0.04% gc time)
Status: success
Candidate solution
Minimizer: [4.09e-01, 5.23e+00, 6.28e+00, ...]
Minimum: 8.416426e-01
Found with
Algorithm: Fminbox with L-BFGS
Initial Point: [4.54e-01, 6.00e+00, 5.54e+00, ...]
Convergence measures
|x - x'| = 0.00e+00 ≰ -1.0e+00
|x - x'|/|x'| = 0.00e+00 ≰ -1.0e+00
|f(x) - f(x')| = 0.00e+00 ≰ -1.0e+00
|f(x) - f(x')|/|f(x')| = 0.00e+00 ≰ -1.0e+00
|g(x)| = 1.53e-03 ≰ 1.0e-06
Work counters
Seconds run: 21 (vs limit Inf)
Iterations: 6
f(x) calls: 11450
∇f(x) calls: 11450```
I recently made a post to the Discourse forum about trying to find a way to access the approximated Hessian from a converged solution to use as the preconditioner on a future minimization of a slightly altered initial state; I'm using BFGS to bifurcate from previously minimizers. After looking at Optim.jl, libLBFGS, and LBFGS++ I couldn't find a interface, so I imagine that it is most likely not available in Optim.jl either. I wanted to double check here if that's something viable with the current interface, if there would be any interest in adding it to the package, or if I should try to extend this on my own in Julia somehow.
So far, the only package outside of Julia that implements this is SciPy:
https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.fmin_bfgs.html
state
object
state
should just be added to MultivariateOptimizationResults
?
optimize(d::D, initial_x::Tx, method::M,
options::Options = Options(;default_options(method)...),
state = initial_state(method, options, d, initial_x))
state
it'll get mutated and you can get what you want from it after
Hi All, new here ... We have been using Optim.jl and LBFGS for a while, and for our particular problems the combination of alphaguess = LineSearches.InitialQuadratic() and linesearch = LineSearches.MoreThuente() works best. I thought I would try that for Rosenbrock on the front page, and it beats BFGS pretty handily, at least in terms of operation count, almost 2x. Thought that was worth mentioning.
Where is a pointer to the slack channel please? Cheers