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

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.cc (modified) (20 diffs)

trunk/SophyaLib/NTools/difeq.cc

@@ -106 +106 @@
 //++
 DiffEqSolver&
-DiffEqSolver::StartV(OVector const& yi, double t)
+DiffEqSolver::StartV(Vector const& yi, double t)
 //
 //      Spécifie le point de départ de l'intégration de l'équadif.
@@ -169 +169 @@
 
 //++
-// virtual void DiffEqSolver::SolveArr(OMatrix&  y, double* t, double tf, int n)=0
+// virtual void DiffEqSolver::SolveArr(Matrix&  y, double* t, double tf, int n)=0
 //      Lance la résolution de l'équadif, jusqu'à la valeur
 //      tf du temps. N valeurs intermédiaires sont retournées dans 
@@ -189 +189 @@
 //++
 void
-DiffEqSolver::SolveV(OVector& yf, double tf)
+DiffEqSolver::SolveV(Vector& yf, double tf)
 //
 //      Lance la résolution de l'équadif, jusqu'à la valeur
@@ -201 +201 @@
 {
   double t;
-  OMatrix m(1, mFunc->NFuncReal());
+  Matrix m(1, mFunc->NFuncReal());
   SolveArr(m, &t, tf, 1);
   yf.Realloc(mFunc->NFunc());
@@ -219 +219 @@
   ASSERT(mFunc->NFunc() == 1);
   double t;
-  OMatrix m(1,mFunc->NFuncReal());
+  Matrix m(1,mFunc->NFuncReal());
   SolveArr(m, &t, tf, 1);
   yf = m(0,0);
@@ -238 +238 @@
 {
   double t;
-  OMatrix m(1, mFunc->NFuncReal());
+  Matrix m(1, mFunc->NFuncReal());
   SolveArr(m, &t, tf, 1);
   for (int i=0; i<mFunc->NFunc(); i++)
@@ -258 +258 @@
 {
   ASSERT(mFunc->NFunc() == 1);
-  OMatrix m(n, mFunc->NFuncReal());
+  Matrix m(n, mFunc->NFuncReal());
   SolveArr(m, t, tf, n);
   for (int i=0; i<n; i++)
@@ -281 +281 @@
 //--
 {
-   OMatrix m(n, mFunc->NFuncReal());
+   Matrix m(n, mFunc->NFuncReal());
    SolveArr(m, t, tf, n);
    for (int i=0; i<n; i++)
@@ -318 +318 @@
 
 //++
-//  virtual void ComputeV(OVector& fpi, OVector const& fi)
+//  virtual void ComputeV(Vector& fpi, Vector const& fi)
 //      Calcule les valeurs des dérivées fpi à partir des valeurs 
 //      des fonctions fi. A redéfinir.
@@ -340 +340 @@
 
 //++
-// virtual void AdjustStart(OVector& start, double tstart)
+// virtual void AdjustStart(Vector& start, double tstart)
 //      Pour ajuster le vecteur de départ quand il y a des 
 //      fonctions à usage interne...
@@ -405 +405 @@
 
 void
-DiffEqFcnT1::ComputeV(OVector& fpi, OVector const& fi)
+DiffEqFcnT1::ComputeV(Vector& fpi, Vector const& fi)
 {
   fpi(0) = (*mFcn)(fi(0), fi(1));
@@ -412 +412 @@
 
 void
-DiffEqFcnT1::AdjustStart(OVector& start, double tstart)
+DiffEqFcnT1::AdjustStart(Vector& start, double tstart)
 {
   start.Realloc(2);
@@ -446 +446 @@
 
 void
-DiffEqFcn2::ComputeV(OVector& fpi, OVector const& fi)
+DiffEqFcn2::ComputeV(Vector& fpi, Vector const& fi)
 {
   fpi(0) = fi(1);
@@ -481 +481 @@
 
 void
-DiffEqFcnT2::ComputeV(OVector& fpi, OVector const& fi)
+DiffEqFcnT2::ComputeV(Vector& fpi, Vector const& fi)
 {
   fpi(0) = fi(1);
@@ -489 +489 @@
 
 void
-DiffEqFcnT2::AdjustStart(OVector& start, double tstart)
+DiffEqFcnT2::AdjustStart(Vector& start, double tstart)
 {
   start.Realloc(3);
@@ -505 +505 @@
 //      un vecteur de dimension 3.
 //      On fournit une fonction de type
-//|     typedef void(*DIFEQFCNV)(OVector& y2, OVector const& y1, 
-//|                              OVector const& y);
+//|     typedef void(*DIFEQFCNV)(Vector& y2, Vector const& y1, 
+//|                              Vector const& y);
 //      qui retourne y'' en fonction de y' et y.
 //      Note : le système résolu est alors en fait
@@ -527 +527 @@
 
 void
-DiffEqFcnV::ComputeV(OVector& fpi, OVector const& fi)
+DiffEqFcnV::ComputeV(Vector& fpi, Vector const& fi)
 {
   fpi(0) = fi(3);
@@ -582 +582 @@
 
 void
-RK4DiffEq::SolveArr(OMatrix& y, double* t, double tf, int n)
+RK4DiffEq::SolveArr(Matrix& y, double* t, double tf, int n)
 {
   //TIMEF;
@@ -592 +592 @@
    double dx = (tf - mXStart)/(n*nStep);
    
-   OVector yt = mYStart;
+   Vector yt = mYStart;
 
    k1.Realloc(mFunc->NFuncReal());
@@ -610 +610 @@
 
 void
-RK4DiffEq::RKStep(OVector& newY, OVector const& y0, double dt)
+RK4DiffEq::RKStep(Vector& newY, Vector const& y0, double dt)
 {
     mFunc->ComputeV(k1, y0);
     k1 *= dt;
-    mFunc->ComputeV(k2, y0 + k1/2);
+    mFunc->ComputeV(k2, y0 + k1/2.);
     k2 *= dt;
-    mFunc->ComputeV(k3, y0 + k2/2);
+    mFunc->ComputeV(k3, y0 + k2/2.);
     k3 *= dt;
     mFunc->ComputeV(k4, y0 + k3);
     k4 *= dt;
     
-    newY = y0 + (k1 + k2*2. + k3*2. + k4)/6;
-}
-
+    newY = y0 + (k1 + k2*2. + k3*2. + k4)/6.;
+}
+