Changeset 307 in Sophya for trunk/SophyaLib/NTools/matrix.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/matrix.cc (modified) (3 diffs)

trunk/SophyaLib/NTools/matrix.cc

@@ -1 @@
-// $Id: matrix.cc,v 1.3 1999-05-03 16:55:21 ansari Exp $
+// $Id: matrix.cc,v 1.4 1999-05-19 15:58:02 ansari Exp $
 
 #include "machdefs.h"
@@ -983 +983 @@
 //////////////////////////////////////////////////////////
 //++
-Matrix* Matrix::FitResidus(GeneralFit& gfit)
+Matrix Matrix::FitResidus(GeneralFit& gfit
+               ,double xorg,double yorg,double dx,double dy)
 //
 //      Retourne une classe contenant les residus du fit ``gfit''.
-//      On suppose que x=j (colonnes) et y=i (lignes) pour m(i,j).
-//--
-{
-if(NCol()<=0||NRows()<=0) return NULL;
+//      Les coordonnees de l'element (i,j) sont :
+//      (i,j) -> x = xorg + j*dx , y = yorg + i*dy
+//--
+{
+if(NCol()<=0||NRows()<=0)
+  throw(SzMismatchError("Matrix::FitResidus: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
-if(f==NULL) return NULL;
+if(f==NULL)
+  throw(NullPtrError("Matrix::FitResidus: NULL pointer\n"));
 Vector par = gfit.GetParm();
-Matrix* m = new Matrix(*this);
+Matrix m(*this);
 for(int i=0;i<NRows();i++) for(int j=0;j<NCol();j++) {
-  double x[2] = {(double)j,(double)i};
-  (*m)(i,j) -= f->Value(x,par.Data());
+  double x[2] = {xorg+j*dx,yorg+i*dy};
+  m(i,j) -= f->Value(x,par.Data());
 }
 return m;
@@ -1002 +1006 @@
 
 //++
-Matrix* Matrix::FitFunction(GeneralFit& gfit)
+Matrix Matrix::FitFunction(GeneralFit& gfit
+               ,double xorg,double yorg,double dx,double dy)
 //
 //      Retourne une classe contenant la fonction du fit ``gfit''.
-//      On suppose que x=j (colonnes) et y=i (lignes) pour m(i,j).
-//--
-{
-if(NCol()<=0||NRows()<=0) return NULL;
+//      Les coordonnees de l'element (i,j) sont :
+//      (i,j) -> x = xorg + j*dx , y = yorg + i*dy
+//--
+{
+if(NCol()<=0||NRows()<=0)
+  throw(SzMismatchError("Matrix::FitFunction: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
-if(f==NULL) return NULL;
+if(f==NULL)
+  throw(NullPtrError("Matrix::FitFunction: NULL pointer\n"));
 Vector par = gfit.GetParm();
-Matrix* m = new Matrix(*this);
+Matrix m(*this);
 for(int i=0;i<NRows();i++) for(int j=0;j<NCol();j++) {
-  double x[2] = {(double)j,(double)i};
-  (*m)(i,j) = f->Value(x,par.Data());
+  double x[2] = {xorg+j*dx,yorg+i*dy};
+  m(i,j) = f->Value(x,par.Data());
 }
 return m;