linesearch.hpp

来自「FreeFem++可以生成高质量的有限元网格。可以用于流体力学」· HPP 代码 · 共 49 行

#ifndef LINE_SEARCH_HH
#define LINE_SEARCH_HH

#include <list>
#include "defs.hpp"
#include "NRJ.hpp"
#include <limits.h>
#include <float.h>

#define MAX_IT_LS	1000

/* 
  This is the base class for one dimensional
  line search algorithms. Some optimization
  algorithms need to find the optimimal point
  on one searching directions.
  Classes derived from this class inherit features of 
  LineSearch. This class cannot be used directly.
*/

// Attention :
// pas de copie
// pas d'op閞ateur =

template <class P, class V,class M,class VM,class R>
class LineSearch { 
public:
	typedef R Real;
	typedef VM VMat;
	typedef M Mat;
	typedef P Param;
	typedef V Vect;
	typedef  tNRJ<Param,Vect,VMat,Real> NRJ;

private:
  //ne sert que si on calcule le gradient de mani鑢e num閞ique
  //dans ce cas, on lui alloue de la m閙oire
  //sinon, ce n'est qu'un pointeur vers le gradient de nrj
  Vect* grad;
  int ref; //ref vaut 1 quand on a cr殚 le gradient de mani鑢e num閞ique

  void numericalGradient(const Param &);


protected:
  //maxinum number of iterations allowed;
  int iterMax;
  //the number iterations so far
  //