patternsearch.h

Go to the documentation of this file.

#ifndef PATTERNSEARCH_H
#define PATTERNSEARCH_H

#include <vec.h>
#include <cirbuffarr.h>
#include <typeop.h>
#include <print.h>

template< typename EXP >
class patternsearch
{
public:

  typedef typename typeop<EXP>::Tbare::Ttype T;
  typedef typename typeop<EXP>::Tbare::XItype XI;

  EXP exp;

  T hstep;

  patternsearch( EXP exp_, T hstep_=1);

  cirbuffarr< T > xi0;

  vec< XI, T > xivec;

  void reset();

  void reset( T const * x0 )
    { exp.reset(x0); reset(); }


  void operator ++ ();

  bool const operator ! () const
    { return exp.valid; }
  
};


/*
 *---------------------------------------------------------------
 * Implementation
 */


template< typename EXP >
void patternsearch<EXP>::reset()
{ 
  exp.reset(); 
  xi0.push(exp.xi); 
  exp.move(); 
  xi0.push(exp.xi); 
}

template< typename EXP >
void patternsearch<EXP>::operator ++ ()
{
  for (bool loop=true; loop;)
  {
    xi0.push(exp.xi);
    xivec -= xi0[1];
    xivec *= hstep;
    xivec += xi0[0];
  
    exp.movelocal();

    loop = exp.evaluate();
  }
  xi0.copyto(exp.xi);
}


template< typename EXP >
patternsearch<EXP>::patternsearch( EXP exp_, T hstep_ )
  : exp(exp_), hstep(hstep_), xi0(3,exp.N), xivec(exp.xi,exp.N)
{
}


#endif