iterative_solver_der.h

Go to the documentation of this file.

#ifndef ITERATIVE_SOLVER_DER_H
#define ITERATIVE_SOLVER_DER_H
#include <iterative_solver.h>
#include <cost_func_diff.h>

namespace FullPhysics {

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

class IterativeSolverDer : 
    public IterativeSolver {

public:


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  IterativeSolverDer(int max_cost_function_calls, bool vrbs)
    : IterativeSolver(max_cost_function_calls, vrbs)
  {}


  virtual ~IterativeSolverDer() {}


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual std::vector< blitz::Array<double, 1> > gradient_at_accepted_points() const
  { return Gradient_at_accepted_points; }


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void print(std::ostream& Os) const 
  { Os << "IterativeSolverDer"; }


protected:

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  void record_gradient_at_accepted_point(const blitz::Array<double, 1>& gradient)
  { Gradient_at_accepted_points.push_back(gradient); }


private:

  std::vector< blitz::Array<double, 1> > Gradient_at_accepted_points;

};
}
#endif