fixed files form Math #71

Spirals-Team · Mar 7, 2017 · 8939597 · 8939597
1 parent 975d277
commit 8939597
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 0 deletions.
diff --git a/projects/Math/71/org/apache/commons/math/ode/nonstiff/EmbeddedRungeKuttaIntegrator.java b/projects/Math/71/org/apache/commons/math/ode/nonstiff/EmbeddedRungeKuttaIntegrator.java
@@ -297,6 +297,10 @@ public double integrate(final FirstOrderDifferentialEquations equations,
                   // it is so small (much probably exactly 0 due to limited accuracy)
                   // that the code above would fail handling it.
                   // So we set up an artificial 0 size step by copying states
+                  interpolator.storeTime(stepStart);
+                  System.arraycopy(y, 0, yTmp, 0, y0.length);
+                  hNew     = 0;
+                  stepSize = 0;
                   loop     = false;
               } else {
                   // reject the step to match exactly the next switch time

diff --git a/projects/Math/71/org/apache/commons/math/ode/nonstiff/RungeKuttaIntegrator.java b/projects/Math/71/org/apache/commons/math/ode/nonstiff/RungeKuttaIntegrator.java
@@ -177,6 +177,9 @@ public double integrate(final FirstOrderDifferentialEquations equations,
                 // it is so small (much probably exactly 0 due to limited accuracy)
                 // that the code above would fail handling it.
                 // So we set up an artificial 0 size step by copying states
+                interpolator.storeTime(stepStart);
+                System.arraycopy(y, 0, yTmp, 0, y0.length);
+                stepSize = 0;
                 loop     = false;
             } else {
                 // reject the step to match exactly the next switch time