#ifndef _LBFGS_H_
#define _LBFGS_H_

// template<class FuncGrad>
// std::vector<double>
// perform_LBFGS(FuncGrad func_grad, const std::vector<double> & x0);

std::vector<double> perform_LBFGS(
    double (*func_grad)(const std::vector<double> &, std::vector<double> &),
    const std::vector<double> &x0);

std::vector<double> perform_OWLQN(
    double (*func_grad)(const std::vector<double> &, std::vector<double> &),
    const std::vector<double> &x0, const double C);

// const int    LBFGS_M = 7;
const int LBFGS_M = 10;

#endif