Changeset 307 in Sophya for trunk/SophyaLib/NTools/cvector.cc

Timestamp:

May 19, 1999, 5:58:06 PM (26 years ago)

Author:

ansari

Message:

FIO_... + grosses modifs cmv 19/5/99

File:

• trunk/SophyaLib/NTools/cvector.cc (modified) (2 diffs)

trunk/SophyaLib/NTools/cvector.cc

@@ -171 +171 @@
 //////////////////////////////////////////////////////////
 //++
-Vector* Vector::FitResidus(GeneralFit& gfit)
+Vector Vector::FitResidus(GeneralFit& gfit,double xorg,double dx)
 //
 //      Retourne une classe contenant les residus du fit ``gfit''.
-//      On suppose que x=i.
-//--
-{
-if(NElts()<=0) return NULL;
+//      La coordonnee de l'element (i) est -> x = xorg + i*dx
+//--
+{
+if(NElts()<=0)
+  throw(SzMismatchError("Vector::FitResidus: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
-if(f==NULL) return NULL;
+if(f==NULL)
+  throw(NullPtrError("Vector::FitResidus: NULL pointer\n"));
 Vector par = gfit.GetParm();
-Vector* v = new Vector(*this);
+Vector v(*this);
 for(int i=0;i<NElts();i++) {
-  double x = (double)i;
-  (*v)(i) -= f->Value(&x,par.Data());
+  double x = xorg+i*dx;
+  v(i) -= f->Value(&x,par.Data());
 }
 return v;
@@ -190 +192 @@
 
 //++
-Vector* Vector::FitFunction(GeneralFit& gfit)
+Vector Vector::FitFunction(GeneralFit& gfit,double xorg,double dx)
 //
 //      Retourne une classe contenant la fonction du fit ``gfit''.
-//      On suppose que x=i.
-//--
-{
-if(NElts()<=0) return NULL;
+//      La coordonnee de l'element (i) est -> x = xorg + i*dx
+//--
+{
+if(NElts()<=0)
+  throw(SzMismatchError("Vector::FitResidus: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
-if(f==NULL) return NULL;
+if(f==NULL)
+  throw(NullPtrError("Vector::FitResidus: NULL pointer\n"));
 Vector par = gfit.GetParm();
-Vector* v = new Vector(*this);
+Vector v(*this);
 for(int i=0;i<NElts();i++) {
-  double x = (double)i;
-  (*v)(i) = f->Value(&x,par.Data());
+  double x = xorg+i*dx;
+  v(i) = f->Value(&x,par.Data());
 }
 return v;