You are viewing a plain text version of this content. The canonical link for it is here.
Posted to commits@commons.apache.org by lu...@apache.org on 2009/06/02 22:43:59 UTC
svn commit: r781158 -
/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/nonstiff/AdamsIntegrator.java
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);
}