A research-oriented implementation and theoretical study of the L-BFGS and L-BFGS-B optimization algorithms.
$$ B_{k+1}s_k = y_k $$
Second-order methods provide fast convergence but require Hessian computation. L-BFGS approximates curvature efficiently using limited memory.