Changeset 301 in Sophya for trunk/SophyaLib/NTools/tmatrix.cc

Timestamp:

May 18, 1999, 2:23:15 PM (26 years ago)

Author:

ansari

Message:

cmv 18/5/99

File:

• trunk/SophyaLib/NTools/tmatrix.cc (modified) (6 diffs)

trunk/SophyaLib/NTools/tmatrix.cc

@@ -1 @@
-// $Id: tmatrix.cc,v 1.5 1999-05-17 17:16:43 ansari Exp $
+// $Id: tmatrix.cc,v 1.6 1999-05-18 12:23:13 ansari Exp $
 //                         C.Magneville          04/99
 #include "machdefs.h"
@@ -66 +66 @@
 TMatrix<T> a; a.Clone(*this); a.SetTemp(true);
 a.mNr = mNc; a.mNc = mNr;
-{for(int i=0; i<mNr; i++)
-  for(int j=0; j<mNc; j++) {
+{for(uint_4 i=0; i<mNr; i++)
+  for(uint_4 j=0; j<mNc; j++) {
     a(j,i) = (*this)(i,j);
 }}
@@ -187 +187 @@
 //////////////////////////////////////////////////////////
 //**** Residus des fits
-TMatrix<r_8> TMatrix<r_8>::FitResidus(GeneralFit& gfit)
-//      Retourne une classe contenant les residus du fit ``gfit''.
-//      On suppose que x=j (colonnes) et y=i (lignes) pour m(i,j).
+TMatrix<r_8> TMatrix<r_8>::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).
+// Les coordonnees de l'element (i,j) sont :
+//    (i,j) -> x = xorg+j*dx , y = yorg+i*dy
 {
 if(NCols()<=0||NRows()<=0)
@@ -202 +205 @@
 {for(int i=0;i<npar;i++) par[i] = gfit.GetParm(i);}
 TMatrix<r_8> m(*this);
-for(int i=0;i<NRows();i++) for(int j=0;j<NCols();j++) {
-  double x[2] = {(double)j,(double)i}; // attention x=col y=row!!
+for(uint_4 i=0;i<NRows();i++) for(uint_4 j=0;j<NCols();j++) {
+  double x[2] = {xorg+j*dx,yorg+i*dy};
   m(i,j) -= f->Value(x,par);
 }
@@ -210 +213 @@
 }
 
-TMatrix<r_8> TMatrix<r_8>::FitFunction(GeneralFit& gfit)
-//      Retourne une classe contenant la fonction du fit ``gfit''.
-//      On suppose que x=j (colonnes) et y=i (lignes) pour m(i,j).
+TMatrix<r_8> TMatrix<r_8>::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).
+// Les coordonnees de l'element (i,j) sont :
+//    (i,j) -> x = xorg + j*dx , y = yorg + i*dy
 {
 if(NCols()<=0||NRows()<=0)
@@ -225 +231 @@
 {for(int i=0;i<npar;i++) par[i] = gfit.GetParm(i);}
 TMatrix<r_8> m(*this);
-for(int i=0;i<NRows();i++) for(int j=0;j<NCols();j++) {
-  double x[2] = {(double)j,(double)i};
+for(uint_4 i=0;i<NRows();i++) for(uint_4 j=0;j<NCols();j++) {
+  double x[2] = {xorg+j*dx,yorg+i*dy};
   m(i,j) = f->Value(x,par);
 }