Changeset 508 in Sophya for trunk/SophyaLib/NTools/linfit.cc

Timestamp:

Oct 25, 1999, 12:36:22 PM (26 years ago)

Author:

ansari

Message:

Vector/Matrix OVector/OMatrix cmv 25/10/99

File:

• trunk/SophyaLib/NTools/linfit.cc (modified) (9 diffs)

trunk/SophyaLib/NTools/linfit.cc

@@ -5 +5 @@
 #include "cvector.h"
 
-double LinFit(const Vector& x, const Vector& y, int nf, double (*f)(int, double), 
-           Vector& c)
+double LinFit(const OVector& x, const OVector& y, int nf, double (*f)(int, double), 
+           OVector& c)
 {
   int n = x.NElts();
   if (n != y.NElts()) THROW(sizeMismatchErr);
   
-  Matrix fx(nf, n);
+  OMatrix fx(nf, n);
   for (int i=0; i<nf; i++)
     for (int j=0; j<n; j++)
@@ -21 +21 @@
 
 
-double LinFit(const Matrix& fx, const Vector& y, Vector& c)
+double LinFit(const OMatrix& fx, const OVector& y, OVector& c)
 {
   int n = y.NElts();
@@ -28 +28 @@
   int nf = fx.NRows();
 
-  Matrix a(nf,nf);
+  OMatrix a(nf,nf);
 
   for (int j=0; j<nf; j++)
@@ -36 +36 @@
   c = fx * y;
 
-  if (Matrix::GausPiv(a,c) == 0) THROW(singMatxErr);
+  if (OMatrix::GausPiv(a,c) == 0) THROW(singMatxErr);
 
   double xi2 = 0;
@@ -52 +52 @@
 
 
-double LinFit(const Vector& x, const Vector& y, const Vector& errY2, int nf,
-           double (*f)(int, double), Vector& c, Vector& errC)
+double LinFit(const OVector& x, const OVector& y, const OVector& errY2, int nf,
+           double (*f)(int, double), OVector& c, OVector& errC)
 {
   int n = x.NElts();
   if (n != y.NElts()) THROW(sizeMismatchErr);
   
-  Matrix fx(nf, n);
+  OMatrix fx(nf, n);
   for (int i=0; i<nf; i++)
     for (int j=0; j<n; j++)
@@ -67 +67 @@
 
 
-double LinFit(const Matrix& fx, const Vector& y, const Vector& errY2,
-           Vector& c, Vector& errC)
+double LinFit(const OMatrix& fx, const OVector& y, const OVector& errY2,
+           OVector& c, OVector& errC)
 {
   int n = y.NElts();
@@ -76 +76 @@
   int nf = fx.NRows();
 
-  Matrix a(nf,nf);
+  OMatrix a(nf,nf);
 
   c.Realloc(nf);
@@ -89 +89 @@
     }
  
-  Matrix d(nf,nf+1);
+  OMatrix d(nf,nf+1);
   for (int k=0; k<nf; k++) {
     double x=0;
@@ -99 +99 @@
   }
 
-  if (Matrix::GausPiv(a,d) == 0) THROW(singMatxErr);
+  if (OMatrix::GausPiv(a,d) == 0) THROW(singMatxErr);
 
   for (int l=0; l<nf; l++) {