Member LSLOpt::Implementation::BFGSStorage< Scalar, OutputFunction >::update (const Vector< Scalar > &s, const Vector< Scalar > &y, const Vector< Scalar > &g): Here we could also introduce a check if the product B * H really is the identity.
Member LSLOpt::Implementation::BoxConstrainedAlgorithm< Problem, Scalar, Storage >::search_direction (Problem &problem, Storage &storage, const Vector< Scalar > &x, const Vector< Scalar > &g, const Scalar &check_epsilon, const Scalar &machine_epsilon, OutputFunction &output_function): It must be checked if this is something that just happens or if this points to some subtle bug.
Member LSLOpt::Implementation::LBFGSStorage< Scalar, OutputFunction >::calculate_vBv (const Vector< Scalar > &v): There is a formulation that calculates this in $\Theta(m^2)$ .
page Main Page: Citation