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 2007/11/07 21:52:54 UTC
svn commit: r592894 - in /commons/proper/math/trunk/src:
java/org/apache/commons/math/ode/ test/org/apache/commons/math/ode/
Author: luc
Date: Wed Nov 7 12:52:53 2007
New Revision: 592894
URL: http://svn.apache.org/viewvc?rev=592894&view=rev
Log:
renamed the RungeKuttaFehlbergIntegrator base class to EmbeddedRungKuttaIntegrator
since it is more general than the single method designed by Fehlberg
renamed some misleading private fields in the corresponding derived classes
Added:
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/EmbeddedRungeKuttaIntegrator.java
- copied, changed from r592878, commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaFehlbergIntegrator.java
Removed:
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaFehlbergIntegrator.java
Modified:
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/AbstractStepInterpolator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54Integrator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54StepInterpolator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853Integrator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853StepInterpolator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DummyStepInterpolator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/GraggBulirschStoerIntegrator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54Integrator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54StepInterpolator.java
commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaStepInterpolator.java
commons/proper/math/trunk/src/test/org/apache/commons/math/ode/DormandPrince54IntegratorTest.java
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/AbstractStepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/AbstractStepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/AbstractStepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/AbstractStepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -70,7 +70,7 @@
* instance in order to initialize the internal arrays. This
* constructor is used only in order to delay the initialization in
* some cases. As an example, the {@link
- * RungeKuttaFehlbergIntegrator} uses the prototyping design pattern
+ * EmbeddedRungeKuttaIntegrator} uses the prototyping design pattern
* to create the step interpolators by cloning an uninitialized
* model and latter initializing the copy.
*/
@@ -296,7 +296,7 @@
/**
* Finalize the step.
- * <p>Some Runge-Kutta-Fehlberg integrators need fewer functions
+ * <p>Some embedded Runge-Kutta integrators need fewer functions
* evaluations than their counterpart step interpolators. These
* interpolators should perform the last evaluations they need by
* themselves only if they need them. This method triggers these
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54Integrator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54Integrator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54Integrator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54Integrator.java Wed Nov 7 12:52:53 2007
@@ -21,7 +21,7 @@
* This class implements the 5(4) Dormand-Prince integrator for Ordinary
* Differential Equations.
- * <p>This integrator is an embedded Runge-Kutta-Fehlberg integrator
+ * <p>This integrator is an embedded Runge-Kutta integrator
* of order 5(4) used in local extrapolation mode (i.e. the solution
* is computed using the high order formula) with stepsize control
* (and automatic step initialization) and continuous output. This
@@ -44,15 +44,15 @@
*/
public class DormandPrince54Integrator
- extends RungeKuttaFehlbergIntegrator {
+ extends EmbeddedRungeKuttaIntegrator {
private static final String methodName = "Dormand-Prince 5(4)";
- private static final double[] c = {
+ private static final double[] staticC = {
1.0/5.0, 3.0/10.0, 4.0/5.0, 8.0/9.0, 1.0, 1.0
};
- private static final double[][] a = {
+ private static final double[][] staticA = {
{1.0/5.0},
{3.0/40.0, 9.0/40.0},
{44.0/45.0, -56.0/15.0, 32.0/9.0},
@@ -61,7 +61,7 @@
{35.0/384.0, 0.0, 500.0/1113.0, 125.0/192.0, -2187.0/6784.0, 11.0/84.0}
};
- private static final double[] b = {
+ private static final double[] staticB = {
35.0/384.0, 0.0, 500.0/1113.0, 125.0/192.0, -2187.0/6784.0, 11.0/84.0, 0.0
};
@@ -84,7 +84,7 @@
public DormandPrince54Integrator(double minStep, double maxStep,
double scalAbsoluteTolerance,
double scalRelativeTolerance) {
- super(true, c, a, b, new DormandPrince54StepInterpolator(),
+ super(true, staticC, staticA, staticB, new DormandPrince54StepInterpolator(),
minStep, maxStep, scalAbsoluteTolerance, scalRelativeTolerance);
}
@@ -100,7 +100,7 @@
public DormandPrince54Integrator(double minStep, double maxStep,
double[] vecAbsoluteTolerance,
double[] vecRelativeTolerance) {
- super(true, c, a, b, new DormandPrince54StepInterpolator(),
+ super(true, staticC, staticA, staticB, new DormandPrince54StepInterpolator(),
minStep, maxStep, vecAbsoluteTolerance, vecRelativeTolerance);
}
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54StepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54StepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54StepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince54StepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -35,7 +35,7 @@
* {@link #reinitialize} method should be called before using the
* instance in order to initialize the internal arrays. This
* constructor is used only in order to delay the initialization in
- * some cases. The {@link RungeKuttaFehlbergIntegrator} uses the
+ * some cases. The {@link EmbeddedRungeKuttaIntegrator} uses the
* prototyping design pattern to create the step interpolators by
* cloning an uninitialized model and latter initializing the copy.
*/
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853Integrator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853Integrator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853Integrator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853Integrator.java Wed Nov 7 12:52:53 2007
@@ -21,7 +21,7 @@
* This class implements the 8(5,3) Dormand-Prince integrator for Ordinary
* Differential Equations.
- * <p>This integrator is an embedded Runge-Kutta-Fehlberg integrator
+ * <p>This integrator is an embedded Runge-Kutta integrator
* of order 8(5,3) used in local extrapolation mode (i.e. the solution
* is computed using the high order formula) with stepsize control
* (and automatic step initialization) and continuous output. This
@@ -52,19 +52,19 @@
*/
public class DormandPrince853Integrator
- extends RungeKuttaFehlbergIntegrator {
+ extends EmbeddedRungeKuttaIntegrator {
private static final String methodName = "Dormand-Prince 8 (5, 3)";
private static final double sqrt6 = Math.sqrt(6.0);
- private static final double[] c = {
+ private static final double[] staticC = {
(12.0 - 2.0 * sqrt6) / 135.0, (6.0 - sqrt6) / 45.0, (6.0 - sqrt6) / 30.0,
(6.0 + sqrt6) / 30.0, 1.0/3.0, 1.0/4.0, 4.0/13.0, 127.0/195.0, 3.0/5.0,
6.0/7.0, 1.0, 1.0
};
- private static final double[][] a = {
+ private static final double[][] staticA = {
// k2
{(12.0 - 2.0 * sqrt6) / 135.0},
@@ -130,7 +130,7 @@
};
- private static final double[] b = {
+ private static final double[] staticB = {
104257.0/1920240.0,
0.0,
0.0,
@@ -176,7 +176,7 @@
public DormandPrince853Integrator(double minStep, double maxStep,
double scalAbsoluteTolerance,
double scalRelativeTolerance) {
- super(true, c, a, b,
+ super(true, staticC, staticA, staticB,
new DormandPrince853StepInterpolator(),
minStep, maxStep, scalAbsoluteTolerance, scalRelativeTolerance);
}
@@ -193,7 +193,7 @@
public DormandPrince853Integrator(double minStep, double maxStep,
double[] vecAbsoluteTolerance,
double[] vecRelativeTolerance) {
- super(true, c, a, b,
+ super(true, staticC, staticA, staticB,
new DormandPrince853StepInterpolator(),
minStep, maxStep, vecAbsoluteTolerance, vecRelativeTolerance);
}
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853StepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853StepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853StepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DormandPrince853StepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -39,7 +39,7 @@
* {@link #reinitialize} method should be called before using the
* instance in order to initialize the internal arrays. This
* constructor is used only in order to delay the initialization in
- * some cases. The {@link RungeKuttaFehlbergIntegrator} uses the
+ * some cases. The {@link EmbeddedRungeKuttaIntegrator} uses the
* prototyping design pattern to create the step interpolators by
* cloning an uninitialized model and latter initializing the copy.
*/
@@ -95,11 +95,11 @@
}
/** Reinitialize the instance
- * Some Runge-Kutta-Fehlberg integrators need fewer functions
+ * Some embedded Runge-Kutta integrators need fewer functions
* evaluations than their counterpart step interpolators. So the
* interpolator should perform the last evaluations they need by
- * themselves. The {@link RungeKuttaFehlbergIntegrator
- * RungeKuttaFehlbergIntegrator} abstract class calls this method in
+ * themselves. The {@link EmbeddedRungeKuttaIntegrator
+ * EmbeddedRungeKuttaIntegrator} abstract class calls this method in
* order to let the step interpolator perform the evaluations it
* needs. These evaluations will be performed during the call to
* <code>doFinalize</code> if any, i.e. only if the step handler
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DummyStepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DummyStepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DummyStepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/DummyStepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -44,7 +44,7 @@
* should be called before using the instance in order to initialize
* the internal arrays. This constructor is used only in order to delay
* the initialization in some cases. As an example, the {@link
- * RungeKuttaFehlbergIntegrator} uses the prototyping design pattern
+ * EmbeddedRungeKuttaIntegrator} uses the prototyping design pattern
* to create the step interpolators by cloning an uninitialized
* model and latter initializing the copy.
*/
Copied: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/EmbeddedRungeKuttaIntegrator.java (from r592878, commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaFehlbergIntegrator.java)
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/EmbeddedRungeKuttaIntegrator.java?p2=commons/proper/math/trunk/src/java/org/apache/commons/math/ode/EmbeddedRungeKuttaIntegrator.java&p1=commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaFehlbergIntegrator.java&r1=592878&r2=592894&rev=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaFehlbergIntegrator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/EmbeddedRungeKuttaIntegrator.java Wed Nov 7 12:52:53 2007
@@ -18,7 +18,7 @@
package org.apache.commons.math.ode;
/**
- * This class implements the common part of all Runge-Kutta-Fehlberg
+ * This class implements the common part of all embedde Runge-Kutta
* integrators for Ordinary Differential Equations.
* <p>These methods are embedded explicit Runge-Kutta methods with two
@@ -50,11 +50,11 @@
* evaluation is saved. For an <i>fsal</i> method, we have cs = 1 and
* asi = bi for all i.</p>
- * @version $Id: RungeKuttaFehlbergIntegrator.java 1705 2006-09-17 19:57:39Z luc $
+ * @version $Id: EmbeddedRungeKuttaIntegrator.java 1705 2006-09-17 19:57:39Z luc $
*/
-public abstract class RungeKuttaFehlbergIntegrator
+public abstract class EmbeddedRungeKuttaIntegrator
extends AdaptiveStepsizeIntegrator {
/** Build a Runge-Kutta integrator with the given Butcher array.
@@ -70,7 +70,7 @@
* @param scalAbsoluteTolerance allowed absolute error
* @param scalRelativeTolerance allowed relative error
*/
- protected RungeKuttaFehlbergIntegrator(boolean fsal,
+ protected EmbeddedRungeKuttaIntegrator(boolean fsal,
double[] c, double[][] a, double[] b,
RungeKuttaStepInterpolator prototype,
double minStep, double maxStep,
@@ -107,7 +107,7 @@
* @param vecAbsoluteTolerance allowed absolute error
* @param vecRelativeTolerance allowed relative error
*/
- protected RungeKuttaFehlbergIntegrator(boolean fsal,
+ protected EmbeddedRungeKuttaIntegrator(boolean fsal,
double[] c, double[][] a, double[] b,
RungeKuttaStepInterpolator prototype,
double minStep, double maxStep,
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/GraggBulirschStoerIntegrator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/GraggBulirschStoerIntegrator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/GraggBulirschStoerIntegrator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/GraggBulirschStoerIntegrator.java Wed Nov 7 12:52:53 2007
@@ -31,7 +31,7 @@
* integration, in order to minimize computation cost. It is
* particularly well suited when a very high precision is needed. The
* limit where this method becomes more efficient than high-order
- * Runge-Kutta-Fehlberg methods like {@link DormandPrince853Integrator
+ * embedded Runge-Kutta methods like {@link DormandPrince853Integrator
* Dormand-Prince 8(5,3)} depends on the problem. Results given in the
* Hairer, Norsett and Wanner book show for example that this limit
* occurs for accuracy around 1e-6 when integrating Saltzam-Lorenz
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54Integrator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54Integrator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54Integrator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54Integrator.java Wed Nov 7 12:52:53 2007
@@ -21,7 +21,7 @@
* This class implements the 5(4) Higham and Hall integrator for
* Ordinary Differential Equations.
- * <p>This integrator is an embedded Runge-Kutta-Fehlberg integrator
+ * <p>This integrator is an embedded Runge-Kutta integrator
* of order 5(4) used in local extrapolation mode (i.e. the solution
* is computed using the high order formula) with stepsize control
* (and automatic step initialization) and continuous output. This
@@ -32,15 +32,15 @@
*/
public class HighamHall54Integrator
- extends RungeKuttaFehlbergIntegrator {
+ extends EmbeddedRungeKuttaIntegrator {
private static final String methodName = "Higham-Hall 5(4)";
- private static final double[] c = {
+ private static final double[] staticC = {
2.0/9.0, 1.0/3.0, 1.0/2.0, 3.0/5.0, 1.0, 1.0
};
- private static final double[][] a = {
+ private static final double[][] staticA = {
{2.0/9.0},
{1.0/12.0, 1.0/4.0},
{1.0/8.0, 0.0, 3.0/8.0},
@@ -49,11 +49,11 @@
{1.0/12.0, 0.0, 27.0/32.0, -4.0/3.0, 125.0/96.0, 5.0/48.0}
};
- private static final double[] b = {
+ private static final double[] staticB = {
1.0/12.0, 0.0, 27.0/32.0, -4.0/3.0, 125.0/96.0, 5.0/48.0, 0.0
};
- private static final double[] e = {
+ private static final double[] staticE = {
-1.0/20.0, 0.0, 81.0/160.0, -6.0/5.0, 25.0/32.0, 1.0/16.0, -1.0/10.0
};
@@ -69,7 +69,7 @@
public HighamHall54Integrator(double minStep, double maxStep,
double scalAbsoluteTolerance,
double scalRelativeTolerance) {
- super(false, c, a, b, new HighamHall54StepInterpolator(),
+ super(false, staticC, staticA, staticB, new HighamHall54StepInterpolator(),
minStep, maxStep, scalAbsoluteTolerance, scalRelativeTolerance);
}
@@ -85,7 +85,7 @@
public HighamHall54Integrator(double minStep, double maxStep,
double[] vecAbsoluteTolerance,
double[] vecRelativeTolerance) {
- super(false, c, a, b, new HighamHall54StepInterpolator(),
+ super(false, staticC, staticA, staticB, new HighamHall54StepInterpolator(),
minStep, maxStep, vecAbsoluteTolerance, vecRelativeTolerance);
}
@@ -117,9 +117,9 @@
double error = 0;
for (int j = 0; j < y0.length; ++j) {
- double errSum = e[0] * yDotK[0][j];
- for (int l = 1; l < e.length; ++l) {
- errSum += e[l] * yDotK[l][j];
+ double errSum = staticE[0] * yDotK[0][j];
+ for (int l = 1; l < staticE.length; ++l) {
+ errSum += staticE[l] * yDotK[l][j];
}
double yScale = Math.max(Math.abs(y0[j]), Math.abs(y1[j]));
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54StepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54StepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54StepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/HighamHall54StepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -35,7 +35,7 @@
* {@link AbstractStepInterpolator#reinitialize} method should be called
* before using the instance in order to initialize the internal arrays. This
* constructor is used only in order to delay the initialization in
- * some cases. The {@link RungeKuttaFehlbergIntegrator} uses the
+ * some cases. The {@link EmbeddedRungeKuttaIntegrator} uses the
* prototyping design pattern to create the step interpolators by
* cloning an uninitialized model and latter initializing the copy.
*/
Modified: commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaStepInterpolator.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaStepInterpolator.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaStepInterpolator.java (original)
+++ commons/proper/math/trunk/src/java/org/apache/commons/math/ode/RungeKuttaStepInterpolator.java Wed Nov 7 12:52:53 2007
@@ -22,11 +22,10 @@
import java.io.IOException;
/** This class represents an interpolator over the last step during an
- * ODE integration for Runge-Kutta and Runge-Kutta-Fehlberg
- * integrators.
+ * ODE integration for Runge-Kutta and embedded Runge-Kutta integrators.
*
* @see RungeKuttaIntegrator
- * @see RungeKuttaFehlbergIntegrator
+ * @see EmbeddedRungeKuttaIntegrator
*
* @version $Id: RungeKuttaStepInterpolator.java 1705 2006-09-17 19:57:39Z luc $
*
@@ -41,7 +40,7 @@
* instance in order to initialize the internal arrays. This
* constructor is used only in order to delay the initialization in
* some cases. The {@link RungeKuttaIntegrator} and {@link
- * RungeKuttaFehlbergIntegrator} classes uses the prototyping design
+ * EmbeddedRungeKuttaIntegrator} classes uses the prototyping design
* pattern to create the step interpolators by cloning an
* uninitialized model and latter initializing the copy.
*/
@@ -97,7 +96,7 @@
* than their counterpart step interpolators. So the interpolator
* should perform the last evaluations they need by themselves. The
* {@link RungeKuttaIntegrator RungeKuttaIntegrator} and {@link
- * RungeKuttaFehlbergIntegrator RungeKuttaFehlbergIntegrator}
+ * EmbeddedRungeKuttaIntegrator EmbeddedRungeKuttaIntegrator}
* abstract classes call this method in order to let the step
* interpolator perform the evaluations it needs. These evaluations
* will be performed during the call to <code>doFinalize</code> if
Modified: commons/proper/math/trunk/src/test/org/apache/commons/math/ode/DormandPrince54IntegratorTest.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/test/org/apache/commons/math/ode/DormandPrince54IntegratorTest.java?rev=592894&r1=592893&r2=592894&view=diff
==============================================================================
--- commons/proper/math/trunk/src/test/org/apache/commons/math/ode/DormandPrince54IntegratorTest.java (original)
+++ commons/proper/math/trunk/src/test/org/apache/commons/math/ode/DormandPrince54IntegratorTest.java Wed Nov 7 12:52:53 2007
@@ -143,7 +143,7 @@
double scalAbsoluteTolerance = Math.pow(10.0, i);
double scalRelativeTolerance = 0.01 * scalAbsoluteTolerance;
- RungeKuttaFehlbergIntegrator integ =
+ EmbeddedRungeKuttaIntegrator integ =
new DormandPrince54Integrator(minStep, maxStep,
scalAbsoluteTolerance, scalRelativeTolerance);
TestProblemHandler handler = new TestProblemHandler(pb, integ);