simplex.H

/* simplex.H */
/* $Id: simplex.H,v 1.3 2000/01/26 13:32:39 hoh Exp $ */
#ifndef __simplex_H
#define __simplex_H

#include <math.h>
#include <mix/wenv.h>
#include <mix/geometry.H>
#include <mix/DoubleAr.H>
#include <mix/IndexAr.H>

#define DEFAULT_LESSER_VALUE 0.37
#define DEFAULT_ENLARGE_VALUE 1.5

#define DEFAULT_A   0.3
#define DEFAULT_B   1.03
#define DEFAULT_G   1.6

class simplex2D{
protected:
  DoubleAr daValues;
  IndexAr iaSorted;
  aPoints2D apCurrent;
  double A,B,G;
  double LESSER_VALUE;
  double ENLARGE_VALUE;
  int iLastReflected;
  double dStep;
  double dSMin;
  double dBValue;

  virtual double value(const point2D&);
  Logic new_top();
  void change_rib(double dChangeCoef);
  void prepare(const point2D &p1,const double dE1,const point2D &p2,
           const double dE2);
public:
  simplex2D(double le=DEFAULT_LESSER_VALUE,double ee=DEFAULT_ENLARGE_VALUE,
        double _a=DEFAULT_A,double _b=DEFAULT_B,double _g=DEFAULT_G);
  virtual void minimize (double &dEps,const double &dMStep,point2D &pRes);

};

#endif /* simplex.H */
/*  Melnikov George Yu. 24-Dec-1999  */ 
/*  Melnikov George Yu. 27-Dec-1999  */ 
/*  Melnikov George Yu. 10-Jan-2000  */