vtkAmoebaMinimizer.h

Go to the documentation of this file.

/*=========================================================================
 
  Program:   Visualization Toolkit
  Module:    vtkAmoebaMinimizer.h
 
  Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
  All rights reserved.
  See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
 
     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notice for more information.
 
=========================================================================*/
#ifndef vtkAmoebaMinimizer_h
#define vtkAmoebaMinimizer_h
 
#include "vtkCommonMathModule.h" // For export macro
#include "vtkObject.h"
 
class VTKCOMMONMATH_EXPORT vtkAmoebaMinimizer : public vtkObject
{
public:
  static vtkAmoebaMinimizer *New();
  vtkTypeMacro(vtkAmoebaMinimizer,vtkObject);
  void PrintSelf(ostream& os, vtkIndent indent);
 
  void SetFunction(void (*f)(void *), void *arg);
 
  void SetFunctionArgDelete(void (*f)(void *));
 
 
  void SetParameterValue(const char *name, double value);
  void SetParameterValue(int i, double value);
 
 
  void SetParameterScale(const char *name, double scale);
  double GetParameterScale(const char *name);
  void SetParameterScale(int i, double scale);
  double GetParameterScale(int i) { return this->ParameterScales[i]; };
 
 
  double GetParameterValue(const char *name);
  double GetParameterValue(int i) { return this->ParameterValues[i]; };
 
  const char *GetParameterName(int i) { return this->ParameterNames[i]; };
 
  int GetNumberOfParameters() { return this->NumberOfParameters; };
 
  void Initialize();
 
  virtual void Minimize();
 
  virtual int Iterate();
 
 
  vtkSetMacro(FunctionValue,double);
  double GetFunctionValue() { return this->FunctionValue; };
 
 
  vtkSetClampMacro(ContractionRatio,double,0.5,1.0);
  vtkGetMacro(ContractionRatio,double);
 
 
  vtkSetClampMacro(ExpansionRatio,double,1.0,2.0);
  vtkGetMacro(ExpansionRatio,double);
 
 
  vtkSetMacro(Tolerance,double);
  vtkGetMacro(Tolerance,double);
 
 
  vtkSetMacro(ParameterTolerance,double);
  vtkGetMacro(ParameterTolerance,double);
 
 
  vtkSetMacro(MaxIterations,int);
  vtkGetMacro(MaxIterations,int);
 
 
  vtkGetMacro(Iterations,int);
 
 
  vtkGetMacro(FunctionEvaluations,int);
 
  void EvaluateFunction();
 
protected:
  vtkAmoebaMinimizer();
  ~vtkAmoebaMinimizer();
 
//BTX
  void (*Function)(void *);
  void (*FunctionArgDelete)(void *);
  void *FunctionArg;
//ETX
 
  int NumberOfParameters;
  char **ParameterNames;
  double *ParameterValues;
  double *ParameterScales;
  double FunctionValue;
 
  double ContractionRatio;
  double ExpansionRatio;
 
  double Tolerance;
  double ParameterTolerance;
  int MaxIterations;
  int Iterations;
  int FunctionEvaluations;
 
private:
// specific to amoeba simplex minimization
//BTX
  double **AmoebaVertices;
  double *AmoebaValues;
  double *AmoebaSum;
  double AmoebaSize;
  double AmoebaHighValue;
  int AmoebaNStepsNoImprovement;
 
  void InitializeAmoeba();
  void GetAmoebaParameterValues();
  void TerminateAmoeba();
  double TryAmoeba(double sum[], int high, double fac);
  int PerformAmoeba();
  int CheckParameterTolerance();
//ETX
 
  vtkAmoebaMinimizer(const vtkAmoebaMinimizer&);  // Not implemented.
  void operator=(const vtkAmoebaMinimizer&);  // Not implemented.
};
 
#endif