commons-commits mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From l..@apache.org
Subject svn commit: r781158 - /commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java
Date Tue, 02 Jun 2009 20:43:59 GMT
Author: luc
Date: Tue Jun  2 20:43:59 2009
New Revision: 781158

URL: http://svn.apache.org/viewvc?rev=781158&view=rev
Log:
fixed event handling during multistep restart phase,
the interpolator used for regular integration and the interpolator used during restart went
out of sync

Modified:
    commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java

Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java?rev=781158&r1=781157&r2=781158&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java
(original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java
Tue Jun  2 20:43:59 2009
@@ -357,6 +357,9 @@
                     // convert to Nordsieck representation
                     scaled    = convertToNordsieckLow();
                     nordsieck = convertToNordsieckHigh(scaled);
+                    interpolator.reinitialize(stepSize, scaled, nordsieck);
+                    interpolator.storeTime(stepStart);
+
 
                 } else {
 
@@ -373,11 +376,11 @@
                     }
                     nordsieck.walkInOptimizedOrder(new NordsieckUpdater(scaled, end, coefficients.c1));
                     scaled = end;
+                    interpolator.reinitialize(stepSize, scaled, nordsieck);
 
                 }
 
                 System.arraycopy(y, 0, yTmp, 0, n);
-                interpolator.reinitialize(stepSize, scaled, nordsieck);
 
             }
 



Mime
View raw message