AstroToolbox/ZeroFinder.h

   1 // ------------------------------------------------------------------------ //
   2 // This source file is part of the 'ESA Advanced Concepts Team's                        //
   3 // Space Mechanics Toolbox' software.                                       //
   4 //                                                                          //
   5 // The source files are for research use only,                              //
   6 // and are distributed WITHOUT ANY WARRANTY. Use them on your own risk.     //
   7 //                                                                          //
   8 // Copyright (c) 2004-2007 European Space Agency                            //
   9 // ------------------------------------------------------------------------ //
  10
  11 #ifndef ZEROFINDER_H
  12 #define ZEROFINDER_H
  13
  14 /// Namespace
  15 namespace ZeroFinder
  16 {
  17
  18         /// Class for one dimensional functions
  19         //  with some parameters
  20         /**  The ()-operator with one double argument
  21          *  has to be overloaded for a derived class
  22          *  The return value is the ordinate computed for
  23          *  the abscissa-argument.
  24          */
  25         class Function1D
  26         {
  27         public:
  28                 //virtual double Compute(double x)=0;
  29                 virtual double operator()(double x)=0;
  30                 // parameters
  31                 double p1,p2;
  32                 void SetParameters(double a, double b);
  33                 virtual ~Function1D() {}
  34         };
  35
  36         class Function1D_7param
  37         {
  38         public:
  39                 //virtual double Compute(double x)=0;
  40                 virtual double operator()(double x)=0;
  41                 // parameters
  42                 double p1,p2,p3,p4,p5,p6,p7;
  43                 void SetParameters(double a, double b, double c, double d, double e, double f, double g);
  44                 virtual ~Function1D_7param() {}
  45         };
  46
  47
  48         class FZero
  49         {
  50                 private:
  51                         double a, c; // lower and upper bound
  52
  53                 public:
  54                         FZero(double a, double b); // constructor
  55                         // fzero procedure
  56                         double FindZero(Function1D& f);
  57                         double FindZero7(Function1D_7param& f);
  58                         void SetBounds(double a, double b);
  59         };
  60
  61
  62 }
  63 #endif