Optim.jl implements an only_fg! function that allows calculations made in the objective function, f, to be passed to the gradient, g. By constructing an fg! function, you can avoid redundant calculations that are required by both f and g.
Does Optimization.jl have something like this?