Changeset 2283 in Sophya for trunk

Timestamp:

Nov 29, 2002, 3:03:45 PM (23 years ago)

Author:

plaszczy

Message:

ajout de classfunc et changements dans integ pour que ca colle

Location:

trunk/SophyaLib/NTools

Files:

• Makefile (modified) (1 diff)
• classfunc.h (added)
• integ.cc (modified) (6 diffs)
• integ.h (modified) (4 diffs)

trunk/SophyaLib/NTools/Makefile

@@ -211 +211 @@
  $(INC)perrors.h fmath.h \
  $(INC)sopemtx.h
-$(OBJ)integ.o: integ.cc integ.h $(INC)machdefs.h \
+$(OBJ)integ.o: integ.cc integ.h classfunc.h $(INC)machdefs.h \
  $(INC)pexceptions.h generalfit.h \
  $(INC)tvector.h \

trunk/SophyaLib/NTools/integ.cc

@@ -54 +54 @@
 {}
 
+Integrator::Integrator(fun f, double xmin, double xmax)
+//
+//      Constructeur à partir de la fonction double->double, et des
+//      bornes d'intégration.
+//--
+: mFunc(new Function(f)), mGFF(NULL), mGFFParm(NULL),
+  mNStep(50), mDX(-1), mReqPrec(-1),
+  mXMin(xmin), mXMax(xmax)
+{}
+
+
 //++
 Integrator::Integrator(FUNC f)
@@ -66 +77 @@
 
 //++
+Integrator::Integrator(fun f)
+// 
+//      Constructeur sans spécifier les bornes. Elles sont positionnées
+//      à [0,1], et on pourra les modifier plus tard.
+//--
+: mFunc(new Function(f)), mGFF(NULL), mGFFParm(NULL),
+  mNStep(50), mDX(-1), mReqPrec(-1),
+  mXMin(0), mXMax(1)
+{}
+
+//++
 Integrator::Integrator(GeneralFunction* gff, double* par, double xmin, double xmax)
 //
@@ -95 +117 @@
 
 Integrator::~Integrator()
-{}
+{delete mFunc;}
 
 //++
@@ -193 +215 @@
   DBASSERT( mFunc || (mGFF && mGFFParm) );
   return mFunc ? (*mFunc)(x) : mGFF->Value(&x, mGFFParm); 
-}
+  }
 
 double
@@ -236 +258 @@
 : Integrator(f, xmin, xmax)
 {}
+//++
+TrpzInteg::TrpzInteg(fun f, double xmin, double xmax)
+//
+//--
+: Integrator(f, xmin, xmax)
+{}
 
 //++
 TrpzInteg::TrpzInteg(FUNC f)
+//
+//--
+: Integrator(f)
+{}
+
+TrpzInteg::TrpzInteg(fun f)
 //
 //--
@@ -331 +365 @@
   NStep(1);
 }
-
+GLInteg::GLInteg(fun f, double xmin, double xmax)
+//
+//--
+: Integrator(f, xmin, xmax), mXPos(NULL), mWeights(NULL), mOrder(8)
+{
+  NStep(1);
+}
 //++
 GLInteg::GLInteg(FUNC f)
+//
+//--
+: Integrator(f), mXPos(NULL), mWeights(NULL), mOrder(8)
+{
+  NStep(1);
+}
+GLInteg::GLInteg(fun f)
 //
 //--

trunk/SophyaLib/NTools/integ.h

@@ -10 +10 @@
 #include "pexceptions.h"
 #include <set>
+#include "classfunc.h"
 
 namespace SOPHYA {
@@ -17 +18 @@
 class Integrator {
 public:
-  typedef double(*FUNC)(double);
+  //typedef double(*FUNC)(double);
+  typedef ClassFunc *FUNC;
+  typedef double (*fun)(double);
 
   Integrator();
   Integrator(FUNC, double xmin, double xmax);
   Integrator(FUNC);
+
+  Integrator(fun, double xmin, double xmax);
+  Integrator(fun);
+
   Integrator(GeneralFunction*, double* par, double xmin, double xmax);
   Integrator(GeneralFunction*, double* par);
@@ -60 +67 @@
   TrpzInteg(FUNC, double xmin, double xmax);
   TrpzInteg(FUNC);
+  TrpzInteg(fun, double xmin, double xmax);
+  TrpzInteg(fun);
   TrpzInteg(GeneralFunction*, double* par, double xmin, double xmax);
   TrpzInteg(GeneralFunction*, double* par);
@@ -73 +82 @@
   GLInteg(FUNC, double xmin, double xmax);
   GLInteg(FUNC);
+  GLInteg(fun, double xmin, double xmax);
+  GLInteg(fun);
   GLInteg(GeneralFunction*, double* par, double xmin, double xmax);
   GLInteg(GeneralFunction*, double* par);