Changeset 514 in Sophya for trunk/SophyaLib/NTools/difeq.h

Timestamp:

Oct 25, 1999, 6:43:04 PM (26 years ago)

Author:

ansari

Message:

elimination des OVector/OMatrix au profit des TVector/TMatrix cmv 25/10/99

File:

• trunk/SophyaLib/NTools/difeq.h (modified) (12 diffs)

trunk/SophyaLib/NTools/difeq.h

@@ -5 +5 @@
 #include "machdefs.h"
 #include "pexceptions.h"
-#include "cvector.h"
-
-namespace PlanckDPC {class GeneralFunction;}
+#include "tvector.h"
+
+namespace PlanckDPC {
+
+class GeneralFunction;
 
 // <summary> fonction pour equadifs </summary>
@@ -28 +30 @@
   
   // Calcule les valeurs des derivees fpi a partir des valeurs des fonctions fi
-  virtual void ComputeV(OVector& fpi, OVector const& fi)
+  virtual void ComputeV(Vector& fpi, Vector const& fi)
        { Compute(fpi(0), fi(0)); }
 
@@ -44 +46 @@
   // Pour ajuster vecteur de depart quand il y a des fonctions a usage
   // interne...
-  virtual void AdjustStart(OVector& /*start*/, double /*tstart*/) 
+  virtual void AdjustStart(Vector& /*start*/, double /*tstart*/) 
     {}   
 protected:
@@ -87 +89 @@
   // Implementation de Compute qui va utiliser la fonction fournie
   // au constructeur.
-  virtual void ComputeV(OVector& fpi, OVector const& fi);
+  virtual void ComputeV(Vector& fpi, Vector const& fi);
 
   // Implementation de AdjustStart qui gere la fonction a usage interne.
-  virtual void AdjustStart(OVector& start, double tstart);
+  virtual void AdjustStart(Vector& start, double tstart);
 protected:
   DIFEQFCNT1 mFcn;
@@ -102 +104 @@
   DiffEqFcn2(DIFEQFCN2);
 
-  virtual void ComputeV(OVector& fpi, OVector const& fi);
+  virtual void ComputeV(Vector& fpi, Vector const& fi);
 protected:
   DIFEQFCN2 mFcn;
@@ -121 +123 @@
   // Implementation de Compute qui va utiliser la fonction fournie
   // au constructeur.
-  virtual void ComputeV(OVector& fpi, OVector const& fi);
+  virtual void ComputeV(Vector& fpi, Vector const& fi);
 
   // Implementation de AdjustStart qui gere la fonction a usage interne.
-  virtual void AdjustStart(OVector& start, double tstart);
+  virtual void AdjustStart(Vector& start, double tstart);
 protected:
   DIFEQFCNT2 mFcn;
@@ -130 +132 @@
 
 // Cas y'' = f(y',y) avec des 3-vecteurs
-typedef void(*DIFEQFCNV)(OVector&, OVector const&, OVector const&);
+typedef void(*DIFEQFCNV)(Vector&, Vector const&, Vector const&);
 
 // <summary> y'' = f(y',y,t) </summary>
 // Cas y'' = f(y',y,t), on fournit la fonction f, sous la forme
-// double f(OVector), et ca construit la bonne DiffEqFunction
+// double f(Vector), et ca construit la bonne DiffEqFunction
 class DiffEqFcnV : public DiffEqFunction {
 public:
-  // Constructeur, on fournit une fonction (OVector)->double
+  // Constructeur, on fournit une fonction (Vector)->double
   // qui donne y'' en fonction du vecteur (t, y, y')
   DiffEqFcnV(DIFEQFCNV);
@@ -143 +145 @@
   // Implementation de Compute qui va utiliser la fonction fournie
   // au constructeur.
-  virtual void ComputeV(OVector& fpi, OVector const& fi);
+  virtual void ComputeV(Vector& fpi, Vector const& fi);
 protected:
   DIFEQFCNV mFcn;
-  OVector tmp1, tmp2, tmp3;
+  Vector tmp1, tmp2, tmp3;
 };
 
@@ -176 +178 @@
   // Change les conditions initiales. Notation chainee possible.
   // <group>
-  DiffEqSolver& StartV(OVector const& yi, double t);
+  DiffEqSolver& StartV(Vector const& yi, double t);
   // si NFunc == 1
   DiffEqSolver& Start1(double        yi, double t);
@@ -184 +186 @@
   // Lance la resolution, avec ou sans conservation de n valeurs intermediaires
   // <group>
-  virtual void SolveV(OVector& yf, double tf);
+  virtual void SolveV(Vector& yf, double tf);
   // si NFunc == 1
   virtual void Solve1(double& yf, double tf);  
   virtual void Solve(double* yf, double tf);
-  virtual void SolveArr(OMatrix&  y, double* t, double tf, int n)=0;
+  virtual void SolveArr(Matrix&  y, double* t, double tf, int n)=0;
   // si NFunc == 1
   virtual void SolveArr1(double*  y, double* t, double tf, int n);   
@@ -198 +200 @@
   bool mOwnFunc;
 
-  OVector  mYStart;
+  Vector  mYStart;
   double  mXStart;
   double  mStep;
@@ -215 +217 @@
 
   // Implementation de RK4
-  virtual void SolveArr(OMatrix& y, double* t, double tf, int n);
+  virtual void SolveArr(Matrix& y, double* t, double tf, int n);
 
 protected:
   // Un pas RK4
-  void RKStep(OVector& newY, OVector const& y0, double dt);
+  void RKStep(Vector& newY, Vector const& y0, double dt);
   // Vecteurs utilises en interne, pour ne pas les reallouer.
-  OVector k1, k2, k3, k4;
-};
-
-
+  Vector k1, k2, k3, k4;
+};
+
+
+} // Fin du namespace
 
 #endif