Changeset 3098 in Sophya for trunk/SophyaLib/NTools/integ.cc

Timestamp:

Nov 2, 2006, 3:32:37 PM (19 years ago)

Author:

ansari

Message:

correction bug lie au delete de ClassFunc ds Integrator / petite amelioration interface Integrator - Reza 02/11/2006

File:

• trunk/SophyaLib/NTools/integ.cc (modified) (16 diffs)

trunk/SophyaLib/NTools/integ.cc

@@ -17 +17 @@
 
   On fournit une fonction double f(double) au constructeur, ou
+  une classe-fonction (ClassFuc) ou 
   une GeneralFunction avec des paramètres définis.
   L'objet Integrator est convertible en valeur double qui est la valeur
@@ -28 +29 @@
 //!  Constructeur par défaut. L'objet n'est pas utilisable en l'état.
 Integrator::Integrator()
-: mFunc(NULL), mGFF(NULL), mGFFParm(NULL),
+: mFunc(NULL), mClFun(NULL), mGFF(NULL), mGFFParm(NULL),
   mNStep(50), mDX(-1), mReqPrec(-1),
   mXMin(0), mXMax(1)
+{}
+
+//! Constructeur à partir d'une classe-fonction, et des bornes d'intégration.
+Integrator::Integrator(ClassFunc const & f, double xmin, double xmax)
+: mFunc(NULL), mClFun(&f), mGFF(NULL), mGFFParm(NULL),
+  mNStep(50), mDX(-1), mReqPrec(-1),
+  mXMin(xmin), mXMax(xmax)
 {}
 
 //! Constructeur à partir de la fonction double->double, et des bornes d'intégration.
 Integrator::Integrator(FUNC f, double xmin, double xmax)
-: mFunc(f), mGFF(NULL), mGFFParm(NULL),
+: mFunc(f), mClFun(NULL), mGFF(NULL), mGFFParm(NULL),
   mNStep(50), mDX(-1), mReqPrec(-1),
   mXMin(xmin), mXMax(xmax)
 {}
 
-//! Constructeur à partir de la fonction double->double, et des bornes d'intégration.
-Integrator::Integrator(fun f, double xmin, double xmax)
-: mFunc(new Function(f)), mGFF(NULL), mGFFParm(NULL),
-  mNStep(50), mDX(-1), mReqPrec(-1),
-  mXMin(xmin), mXMax(xmax)
-{}
-
+
+/*!
+  Constructeur a partir d'une classe-fonction 
+  sans spécifier les bornes. Elles sont positionnées
+  à [0,1], et on pourra les modifier plus tard.
+*/
+Integrator::Integrator(ClassFunc const & f)
+: mFunc(NULL), mClFun(&f), mGFF(NULL), mGFFParm(NULL),
+  mNStep(50), mDX(-1), mReqPrec(-1),
+  mXMin(0), mXMax(1)
+{}
 
 /*!
@@ -53 +65 @@
 */
 Integrator::Integrator(FUNC f)
-: mFunc(f), mGFF(NULL), mGFFParm(NULL),
-  mNStep(50), mDX(-1), mReqPrec(-1),
-  mXMin(0), mXMax(1)
-{}
-
-/*!
-  Constructeur sans spécifier les bornes. Elles sont positionnées
-  à [0,1], et on pourra les modifier plus tard.
-*/
-Integrator::Integrator(fun f)
-: mFunc(new Function(f)), mGFF(NULL), mGFFParm(NULL),
+: mFunc(f), mClFun(NULL), mGFF(NULL), mGFFParm(NULL),
   mNStep(50), mDX(-1), mReqPrec(-1),
   mXMin(0), mXMax(1)
@@ -74 +76 @@
 */
 Integrator::Integrator(GeneralFunction* gff, double* par, double xmin, double xmax)
-: mFunc(NULL), mGFF(gff), mGFFParm(par),
+: mFunc(NULL), mClFun(NULL), mGFF(gff), mGFFParm(par),
   mNStep(50), mDX(-1), mReqPrec(-1),
   mXMin(xmin), mXMax(xmax)
@@ -88 +90 @@
 */
 Integrator::Integrator(GeneralFunction* gff, double* par)
-: mFunc(NULL), mGFF(gff), mGFFParm(par),
+: mFunc(NULL), mClFun(NULL), mGFF(gff), mGFFParm(par),
   mNStep(50), mDX(-1), mReqPrec(-1),
   mXMin(0), mXMax(1)
@@ -96 +98 @@
 
 Integrator::~Integrator()
-{if(mFunc) delete mFunc;}
+{
+}
+
 
 /*! 
@@ -121 +125 @@
 {
   mDX = d;
-  mNStep = mReqPrec = -1;
+  mNStep = -1;
+  mReqPrec = -1.;
   StepsChanged();
   return *this;
@@ -154 +159 @@
 //! Spécifie la fonction à intégrer, sous forme double f(double).
 Integrator& 
-Integrator::Func(FUNC f)
-{
-  mFunc = f;
+Integrator::SetFunc(ClassFunc const & f)
+{
+  mFunc = NULL;
+  mClFun = &f;
   mGFFParm = NULL;
   mFunc = NULL;
@@ -163 +169 @@
 }
 
+//! Spécifie la fonction à intégrer, sous forme double f(double).
+Integrator& 
+Integrator::SetFunc(FUNC f)
+{
+  mFunc = f;
+  mClFun = NULL;
+  mGFFParm = NULL;
+  mFunc = NULL;
+  FuncChanged();
+  return *this;
+}
+
 /*!
   \brief Spécifie la fonction à intégrer, sous forme de GeneralFunction 
@@ -168 +186 @@
 */
 Integrator& 
-Integrator::Func(GeneralFunction* gff, double* par)
+Integrator::SetFunc(GeneralFunction* gff, double* par)
 {
   mGFF = gff;
@@ -181 +199 @@
 Integrator::FVal(double x) const
 {
-  DBASSERT( mFunc || (mGFF && mGFFParm) );
-  return mFunc ? (*mFunc)(x) : mGFF->Value(&x, mGFFParm); 
-  }
+  DBASSERT( mFunc || mClFun || (mGFF && mGFFParm) );
+  if (mFunc) return mFunc(x);
+  else return mClFun ? (*mClFun)(x) : mGFF->Value(&x, mGFFParm); 
+}
 
 double
@@ -197 +216 @@
  cout<<"Integrator between "<<mXMin<<" and "<<mXMax
      <<" in "<<mNStep<<" steps"<<endl;
- if(lp>1) cout<<"mFunc="<<mFunc<<"mGFF="<<mGFF
+ if(lp>1) cout<<"mFunc="<<mFunc<<"mClFun="<<mClFun<<"mGFF="<<mGFF
               <<"mReqPrec="<<mReqPrec<<"mDX="<<mDX<<endl;
 }
@@ -220 +239 @@
 {}
 
+TrpzInteg::TrpzInteg(ClassFunc const & f, double xmin, double xmax)
+: Integrator(f, xmin, xmax)
+{}
+
 TrpzInteg::TrpzInteg(FUNC f, double xmin, double xmax)
 : Integrator(f, xmin, xmax)
 {}
 
-TrpzInteg::TrpzInteg(fun f, double xmin, double xmax)
-: Integrator(f, xmin, xmax)
+TrpzInteg::TrpzInteg(ClassFunc const & f)
+: Integrator(f)
 {}
 
 TrpzInteg::TrpzInteg(FUNC f)
-: Integrator(f)
-{}
-
-TrpzInteg::TrpzInteg(fun f)
 : Integrator(f)
 {}
@@ -301 +320 @@
 
 
+GLInteg::GLInteg(ClassFunc const & f, double xmin, double xmax)
+: Integrator(f, xmin, xmax), mXPos(NULL), mWeights(NULL), mOrder(8)
+{
+  NStep(1);
+}
+
 GLInteg::GLInteg(FUNC f, double xmin, double xmax)
 : Integrator(f, xmin, xmax), mXPos(NULL), mWeights(NULL), mOrder(8)
@@ -307 +332 @@
 }
 
-GLInteg::GLInteg(fun f, double xmin, double xmax)
-: Integrator(f, xmin, xmax), mXPos(NULL), mWeights(NULL), mOrder(8)
+GLInteg::GLInteg(ClassFunc const & f)
+: Integrator(f), mXPos(NULL), mWeights(NULL), mOrder(8)
 {
   NStep(1);
@@ -314 +339 @@
 
 GLInteg::GLInteg(FUNC f)
-: Integrator(f), mXPos(NULL), mWeights(NULL), mOrder(8)
-{
-  NStep(1);
-}
-
-GLInteg::GLInteg(fun f)
 : Integrator(f), mXPos(NULL), mWeights(NULL), mOrder(8)
 {