Changeset 302 in Sophya for trunk/SophyaLib/NTools

Timestamp:

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

Author:

ansari

Message:

cmv 18/5/99

Location:

trunk/SophyaLib/NTools

Files:

• tmatrix.cc (modified) (4 diffs)
• tmatrix.h (modified) (1 diff)
• tvector.cc (modified) (4 diffs)

trunk/SophyaLib/NTools/tmatrix.cc

@@ -1 @@
-// $Id: tmatrix.cc,v 1.6 1999-05-18 12:23:13 ansari Exp $
+// $Id: tmatrix.cc,v 1.7 1999-05-18 12:58:24 ansari Exp $
 //                         C.Magneville          04/99
 #include "machdefs.h"
@@ -9 +9 @@
 #include "tmatrix.h"
 #include "objfio.h"
-#include "generalfit.h"
 
 using namespace PlanckDPC;
@@ -166 +165 @@
 ////////////////////////////////////////////////////////////////
 //**** Pour inversion
+r_8 TMatrix<r_8>::GausPiv(TMatrix<r_8>& a, TMatrix<r_8>& b)
+// Pivot de Gauss
+// * Attention: egcs impose que cette fonction soit mise dans le .cc
+//              avant ::Inverse() (car Inverse() l'utilise)
+{
+Matrix A(a);
+Matrix B(b);
+return (r_8) Matrix::GausPiv(A,B);
+}
+
 TMatrix<r_8> TMatrix<r_8>::Inverse() const
 // Inversion
@@ -177 +186 @@
 }
 
-double TMatrix<r_8>::GausPiv(TMatrix<r_8>& a, TMatrix<r_8>& b)
-// Pivot de Gauss
-{
-Matrix A(a);
-Matrix B(b);
-return Matrix::GausPiv(A,B);
-}
-
+#include "generalfit.h"
 //////////////////////////////////////////////////////////
 //**** Residus des fits

trunk/SophyaLib/NTools/tmatrix.h

@@ -39 +39 @@
 
   // Informations pointeur/data
-  inline int NRows() const {return mNr;}
-  inline int NCols() const {return mNc;}
+  inline uint_4 NRows() const {return mNr;}
+  inline uint_4 NCols() const {return mNc;}
   inline T const& operator()(uint_4 r,uint_4 c) const
                             {return *(mNDBlock.Begin()+r*mNc+c);}

trunk/SophyaLib/NTools/tvector.cc

@@ -1 @@
-// $Id: tvector.cc,v 1.2 1999-05-18 12:23:15 ansari Exp $
+// $Id: tvector.cc,v 1.3 1999-05-18 12:58:26 ansari Exp $
 //                         C.Magneville          04/99
 #include "machdefs.h"
@@ -9 +9 @@
 #include "tvector.h"
 #include "objfio.h"
-#include "generalfit.h"
 
 using namespace PlanckDPC;
@@ -51 +50 @@
 }
 
+#include "generalfit.h"
 //////////////////////////////////////////////////////////
 //**** Residus des fits
@@ -68 +68 @@
 {for(int i=0;i<npar;i++) par[i] = gfit.GetParm(i);}
 TVector<r_8> v(*this);
-for(int i=0;i<NElts();i++) {
+for(uint_4 i=0;i<NElts();i++) {
   double x = xorg+i*dx;
   v(i) -= f->Value(&x,par);