Changeset 508 in Sophya for trunk/SophyaLib/NTools/cvector.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/cvector.cc (modified) (14 diffs)

trunk/SophyaLib/NTools/cvector.cc

@@ -6 +6 @@
 
 //++
-// Class        Vector
+// Class        OVector
 // Lib  Outils++
 // include      cvector.h
@@ -16 +16 @@
 //++
 // Links        Parents
-// Matrix
+// OMatrix
 //--
 
@@ -24 +24 @@
 
 //++
-Vector::Vector(int n)
+OVector::OVector(int n)
 //
 //      Constructeur, n = nombre d'éléments.
 //--
-: Matrix(n, 1)
+: OMatrix(n, 1)
 {END_CONSTRUCTOR}
 
 
 //++
-Vector::Vector(int n, double* values)
+OVector::OVector(int n, double* values)
 //
 //      Constructeur, à partir des valeurs. Pas d'allocation.
 //--
-: Matrix(n, 1, values)
+: OMatrix(n, 1, values)
 {END_CONSTRUCTOR}
 
 
 //++
-Vector::Vector(const Vector& v)
+OVector::OVector(const OVector& v)
 //
 //      Constructeur par copie.
 //--
-: Matrix(v)
+: OMatrix(v)
 {END_CONSTRUCTOR}
 
 
 //++
-Vector::Vector(const Matrix& a)
+OVector::OVector(const OMatrix& a)
 //
 //      Constructeur par "copie" à partir d'une matrice, qui doit avoir une seule colonne.
@@ -58 +58 @@
 //      sizeMismatchErr    si la matrice a plus d'une colonne
 //--
-: Matrix(a.nr,1)      // Matrix(a) ne marche pas, pourquoi ???? bug gcc ???
+: OMatrix(a.nr,1)      // OMatrix(a) ne marche pas, pourquoi ???? bug gcc ???
 {
   if (a.nc != 1) THROW(sizeMismatchErr);
@@ -66 +66 @@
 
 //++
-Vector::Vector(const TVector<r_8>& v)
+OVector::OVector(const TVector<r_8>& v)
 //
 //      Constructeur par "copie" a partir d'un TVector<r_8>.
 //      Attention, les donnees sont partagees.
 //--
-: Matrix(v)
+: OMatrix(v)
 {
 }
@@ -87 +87 @@
 
 //++
-Vector& Vector::operator = (double x)
+OVector& OVector::operator = (double x)
 //
 //      Affectation à partir d'un scalaire. Tous les éléments prennent cette valeur.
@@ -98 +98 @@
 
 //++
-double operator* (const Vector& v1, const Vector& v2)
+double operator* (const OVector& v1, const OVector& v2)
 //
 //      Produit scalaire entre deux vecteurs.
@@ -117 +117 @@
 
 //++
-// Vector operator* (const Vector& v, double b)
+// OVector operator* (const OVector& v, double b)
 //      multiplication par un scalaire.
 //--
 
 //++
-// Vector operator* (const Matrix& a, const Vector& b)
+// OVector operator* (const OMatrix& a, const OVector& b)
 //      produit matrice*vecteur.
 //--
@@ -144 +144 @@
 
 //++
-double LinSolveInPlace(Matrix& a, Vector& b)
+double LinSolveInPlace(OMatrix& a, OVector& b)
 //
 //      Résolution du système A*C = B
@@ -152 +152 @@
   if (a.NCol() != a.NRows()) THROW(sizeMismatchErr);
 
-  return Matrix::GausPiv(a,b);
-}
-
-
-//++
-double LinSolve(const Matrix& a, const Vector& b, Vector& c)
+  return OMatrix::GausPiv(a,b);
+}
+
+
+//++
+double LinSolve(const OMatrix& a, const OVector& b, OVector& c)
 //
 //      Résolution du système A*C = B, avec C retourné dans B
@@ -165 +165 @@
   if (a.NCol() != a.NRows()) THROW(sizeMismatchErr);
   c = b;
-  Matrix a1(a);
-  return Matrix::GausPiv(a1,c);
+  OMatrix a1(a);
+  return OMatrix::GausPiv(a1,c);
 }
 
 //////////////////////////////////////////////////////////
 //++
-Vector Vector::FitResidus(GeneralFit& gfit,double xorg,double dx)
+OVector OVector::FitResidus(GeneralFit& gfit,double xorg,double dx)
 //
 //      Retourne une classe contenant les residus du fit ``gfit''.
@@ -178 +178 @@
 {
 if(NElts()<=0)
-  throw(SzMismatchError("Vector::FitResidus: size mismatch\n"));
+  throw(SzMismatchError("OVector::FitResidus: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
 if(f==NULL)
-  throw(NullPtrError("Vector::FitResidus: NULL pointer\n"));
-Vector par = gfit.GetParm();
-Vector v(*this);
+  throw(NullPtrError("OVector::FitResidus: NULL pointer\n"));
+OVector par = gfit.GetParm();
+OVector v(*this);
 for(int i=0;i<NElts();i++) {
   double x = xorg+i*dx;
@@ -192 +192 @@
 
 //++
-Vector Vector::FitFunction(GeneralFit& gfit,double xorg,double dx)
+OVector OVector::FitFunction(GeneralFit& gfit,double xorg,double dx)
 //
 //      Retourne une classe contenant la fonction du fit ``gfit''.
@@ -199 +199 @@
 {
 if(NElts()<=0)
-  throw(SzMismatchError("Vector::FitResidus: size mismatch\n"));
+  throw(SzMismatchError("OVector::FitResidus: size mismatch\n"));
 GeneralFunction* f = gfit.GetFunction();
 if(f==NULL)
-  throw(NullPtrError("Vector::FitResidus: NULL pointer\n"));
-Vector par = gfit.GetParm();
-Vector v(*this);
+  throw(NullPtrError("OVector::FitResidus: NULL pointer\n"));
+OVector par = gfit.GetParm();
+OVector v(*this);
 for(int i=0;i<NElts();i++) {
   double x = xorg+i*dx;