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

Timestamp:

Aug 28, 2003, 3:29:03 PM (22 years ago)

Author:

cmv

Message:

Quelques ameliorations mineures cmv 28/8/2003

File:

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

trunk/SophyaLib/NTools/integ.cc

@@ -117 +117 @@
 
 Integrator::~Integrator()
-{delete mFunc;}
+{if(mFunc) delete mFunc;}
 
 //++
@@ -222 +222 @@
   Limits(xmin,xmax);
   return Value();
+}
+
+void
+Integrator::Print(int lp)
+{
+ cout<<"Integrator between "<<mXMin<<" and "<<mXMax
+     <<" in "<<mNStep<<" steps"<<endl;
+ if(lp>1) cout<<"mFunc="<<mFunc<<"mGFF="<<mGFF
+              <<"mReqPrec="<<mReqPrec<<"mDX="<<mDX<<endl;
 }
 
@@ -332 +341 @@
 //      fonction aux zéros des polynomes de Legendre. Avec les poids
 //      qui vont bien, GL d'ordre n est exacte pour des polynomes de 
-//      degré <= 2n+1.
+//      degré <= 2n+1 (monome le + haut x^(2*n-1).
 //      Impossible de demander une précision donnée.
 //
@@ -408 +417 @@
 GLInteg::~GLInteg()
 {
-  delete[] mXPos;
-  delete[] mWeights;
-  mXPos = mWeights = NULL;
+ if(mXPos)    delete[] mXPos;     mXPos    = NULL;
+ if(mWeights) delete[] mWeights;  mWeights = NULL;
 }
 
@@ -416 +424 @@
 GLInteg::InvalWeights()
 {
-  delete[] mXPos;
-  delete[] mWeights;
-  mXPos = mWeights = NULL;
+ if(mXPos)    delete[] mXPos;     mXPos    = NULL;
+ if(mWeights) delete[] mWeights;  mWeights = NULL;
 }
 
@@ -425 +432 @@
 GLInteg::LimitsChanged()
 {
-  ComputeBounds();
+ ComputeBounds();
 }
 
@@ -431 +438 @@
 GLInteg::StepsChanged()
 {
-  ComputeBounds();
+ ComputeBounds();
 }
 
@@ -437 +444 @@
 GLInteg::AddBound(double x)
 {
-  if (x<=mXMin || x>=mXMax) THROW(rangeCheckErr); 
-  // On introduira les classes d'exections apres reflexion et de maniere systematique (Rz+cmv)
-  // if (x<=mXMin || x>=mXMax) throw range_error("GLInteg::AddBound  bound outside interval");
-  mBounds.insert(x);
-  return *this;
+ if (x<=mXMin || x>=mXMax) THROW(rangeCheckErr); 
+ // On introduira les classes d'exections apres reflexion et de maniere systematique (Rz+cmv)
+ // if (x<=mXMin || x>=mXMax) throw range_error("GLInteg::AddBound  bound outside interval");
+ mBounds.insert(x);
+ return *this;
 }
 
@@ -448 +455 @@
 GLInteg::SetOrder(int order)
 {
-  mOrder = order;
-  InvalWeights();
-  return *this;
+ mOrder = order;
+ InvalWeights();
+ return *this;
 }
 
@@ -458 +465 @@
 GLInteg::Value()
 {
-  if (!mXPos) ComputeWeights();
-  double s = 0;
-  set<double>::iterator i=mBounds.begin();
-  set<double>::iterator j=i; j++;
-  while (j != mBounds.end()) {
-    double s1 = 0;
-    double x1 = *i;
-    double x2 = *j;
-    for (int k=0; k<mOrder; k++)
-       s1 += mWeights[k] * FVal(x1 + (x2-x1)*mXPos[k]);
-    s += s1*(x2-x1);
-    i++; j++;
-  }
-  return s;
+ if (!mXPos) ComputeWeights();
+ double s = 0;
+ set<double>::iterator i=mBounds.begin();
+ set<double>::iterator j=i; j++;
+ while (j != mBounds.end()) {
+   double s1 = 0;
+   double x1 = *i;
+   double x2 = *j;
+   for(int k=0; k<mOrder; k++)
+      s1 += mWeights[k] * FVal(x1 + (x2-x1)*mXPos[k]);
+   s += s1*(x2-x1);
+   i++; j++;
+ }
+ return s;
 }
 
@@ -485 +492 @@
 GLInteg::ComputeWeights()
 {
-   const double EPS_gauleg = 3.0e-11;
+   const double EPS_gauleg = 5.0e-12;
 
    int m=(mOrder+1)/2;
@@ -519 +526 @@
 }
 
+void
+GLInteg::Print(int lp)
+{
+ Integrator::Print(lp);
+ cout<<"GLInteg order="<<mOrder<<endl;
+ if(lp>0 && mOrder>0) {
+   for(int i=0;i<mOrder;i++)
+     cout<<" ("<<mXPos[i]<<","<<mWeights[i]<<")";
+   cout<<endl;
+ }
+}