//////////////////////////////////////////////////////////////////////
// InnerOptimizationWrapperLBFGS.hpp
//
// Inner optimization algorithm.
//////////////////////////////////////////////////////////////////////

#ifndef INNEROPTIMIZATIONWRAPPERLBFGS_HPP
#define INNEROPTIMIZATIONWRAPPERLBFGS_HPP

#include "OptimizationWrapper.hpp"
#include "LBFGS.hpp"

template<class RealT>
class LBFGS;

template<class RealT>
class OptimizationWrapper;

template<class RealT>
class InnerOptimizationWrapper;

//////////////////////////////////////////////////////////////////////
// class InnerOptimizationWrapperLBFGS
//////////////////////////////////////////////////////////////////////

template<class RealT>
class InnerOptimizationWrapperLBFGS : public LBFGS<RealT>, public InnerOptimizationWrapper<RealT>
{
    RealT log_base;
    
public:
    InnerOptimizationWrapperLBFGS(OptimizationWrapper<RealT> *optimization_wrapper,
                                  const std::vector<int> &units,
                                  const std::vector<RealT> &C);
    
    RealT ComputeFunction(const std::vector<RealT> &x);
    void ComputeGradient(std::vector<RealT> &g, const std::vector<RealT> &x);
    void Report(int iteration, const std::vector<RealT> &x, RealT f, RealT step_size);
    void Report(const std::string &s);
    RealT Minimize(std::vector<RealT> &x0);
};

#include "InnerOptimizationWrapperLBFGS.ipp"

#endif