Class DormandPrince853StepInterpolator
- java.lang.Object
-
- org.apache.commons.math3.ode.sampling.AbstractStepInterpolator
-
- org.apache.commons.math3.ode.nonstiff.RungeKuttaStepInterpolator
-
- org.apache.commons.math3.ode.nonstiff.DormandPrince853StepInterpolator
-
- All Implemented Interfaces:
java.io.Externalizable,java.io.Serializable,StepInterpolator
class DormandPrince853StepInterpolator extends RungeKuttaStepInterpolator
This class represents an interpolator over the last step during an ODE integration for the 8(5,3) Dormand-Prince integrator.- Since:
- 1.2
- See Also:
DormandPrince853Integrator
-
-
Field Summary
Fields Modifier and Type Field Description private static doubleB_01Propagation weights, element 1.private static doubleB_06Propagation weights, element 6.private static doubleB_07Propagation weights, element 7.private static doubleB_08Propagation weights, element 8.private static doubleB_09Propagation weights, element 9.private static doubleB_10Propagation weights, element 10.private static doubleB_11Propagation weights, element 11.private static doubleB_12Propagation weights, element 12.private static doubleC14Time step for stage 14 (interpolation only).private static doubleC15Time step for stage 15 (interpolation only).private static doubleC16Time step for stage 16 (interpolation only).private static double[][]DInterpolation weights.private static doubleK14_01Internal weights for stage 14, element 1.private static doubleK14_06Internal weights for stage 14, element 6.private static doubleK14_07Internal weights for stage 14, element 7.private static doubleK14_08Internal weights for stage 14, element 8.private static doubleK14_09Internal weights for stage 14, element 9.private static doubleK14_10Internal weights for stage 14, element 10.private static doubleK14_11Internal weights for stage 14, element 11.private static doubleK14_12Internal weights for stage 14, element 12.private static doubleK14_13Internal weights for stage 14, element 13.private static doubleK15_01Internal weights for stage 15, element 1.private static doubleK15_06Internal weights for stage 15, element 6.private static doubleK15_07Internal weights for stage 15, element 7.private static doubleK15_08Internal weights for stage 15, element 8.private static doubleK15_09Internal weights for stage 15, element 9.private static doubleK15_10Internal weights for stage 15, element 10.private static doubleK15_11Internal weights for stage 15, element 11.private static doubleK15_12Internal weights for stage 15, element 12.private static doubleK15_13Internal weights for stage 15, element 13.private static doubleK15_14Internal weights for stage 15, element 14.private static doubleK16_01Internal weights for stage 16, element 1.private static doubleK16_06Internal weights for stage 16, element 6.private static doubleK16_07Internal weights for stage 16, element 7.private static doubleK16_08Internal weights for stage 16, element 8.private static doubleK16_09Internal weights for stage 16, element 9.private static doubleK16_10Internal weights for stage 16, element 10.private static doubleK16_11Internal weights for stage 16, element 11.private static doubleK16_12Internal weights for stage 16, element 12.private static doubleK16_13Internal weights for stage 16, element 13.private static doubleK16_14Internal weights for stage 16, element 14.private static doubleK16_15Internal weights for stage 16, element 15.private static longserialVersionUIDSerializable version identifier.private double[][]vVectors for interpolation.private booleanvectorsInitializedInitialization indicator for the interpolation vectors.private double[][]yDotKLastLast evaluations.-
Fields inherited from class org.apache.commons.math3.ode.nonstiff.RungeKuttaStepInterpolator
integrator, previousState, yDotK
-
Fields inherited from class org.apache.commons.math3.ode.sampling.AbstractStepInterpolator
currentState, h, interpolatedDerivatives, interpolatedPrimaryDerivatives, interpolatedPrimaryState, interpolatedSecondaryDerivatives, interpolatedSecondaryState, interpolatedState, interpolatedTime
-
-
Constructor Summary
Constructors Constructor Description DormandPrince853StepInterpolator()Simple constructor.DormandPrince853StepInterpolator(DormandPrince853StepInterpolator interpolator)Copy constructor.
-
Method Summary
All Methods Instance Methods Concrete Methods Modifier and Type Method Description protected voidcomputeInterpolatedStateAndDerivatives(double theta, double oneMinusThetaH)Compute the state and derivatives at the interpolated time.protected StepInterpolatordoCopy()Really copy the finalized instance.protected voiddoFinalize()Really finalize the step.voidreadExternal(java.io.ObjectInput in)voidreinitialize(AbstractIntegrator integrator, double[] y, double[][] yDotK, boolean forward, EquationsMapper primaryMapper, EquationsMapper[] secondaryMappers)Reinitialize the instancevoidstoreTime(double t)Store the current step time.voidwriteExternal(java.io.ObjectOutput out)-
Methods inherited from class org.apache.commons.math3.ode.nonstiff.RungeKuttaStepInterpolator
shift
-
Methods inherited from class org.apache.commons.math3.ode.sampling.AbstractStepInterpolator
copy, finalizeStep, getCurrentTime, getGlobalCurrentTime, getGlobalPreviousTime, getInterpolatedDerivatives, getInterpolatedSecondaryDerivatives, getInterpolatedSecondaryState, getInterpolatedState, getInterpolatedTime, getPreviousTime, isForward, readBaseExternal, reinitialize, setInterpolatedTime, setSoftCurrentTime, setSoftPreviousTime, writeBaseExternal
-
-
-
-
Field Detail
-
serialVersionUID
private static final long serialVersionUID
Serializable version identifier.- See Also:
- Constant Field Values
-
B_01
private static final double B_01
Propagation weights, element 1.- See Also:
- Constant Field Values
-
B_06
private static final double B_06
Propagation weights, element 6.- See Also:
- Constant Field Values
-
B_07
private static final double B_07
Propagation weights, element 7.- See Also:
- Constant Field Values
-
B_08
private static final double B_08
Propagation weights, element 8.- See Also:
- Constant Field Values
-
B_09
private static final double B_09
Propagation weights, element 9.- See Also:
- Constant Field Values
-
B_10
private static final double B_10
Propagation weights, element 10.- See Also:
- Constant Field Values
-
B_11
private static final double B_11
Propagation weights, element 11.- See Also:
- Constant Field Values
-
B_12
private static final double B_12
Propagation weights, element 12.- See Also:
- Constant Field Values
-
C14
private static final double C14
Time step for stage 14 (interpolation only).- See Also:
- Constant Field Values
-
K14_01
private static final double K14_01
Internal weights for stage 14, element 1.- See Also:
- Constant Field Values
-
K14_06
private static final double K14_06
Internal weights for stage 14, element 6.- See Also:
- Constant Field Values
-
K14_07
private static final double K14_07
Internal weights for stage 14, element 7.- See Also:
- Constant Field Values
-
K14_08
private static final double K14_08
Internal weights for stage 14, element 8.- See Also:
- Constant Field Values
-
K14_09
private static final double K14_09
Internal weights for stage 14, element 9.- See Also:
- Constant Field Values
-
K14_10
private static final double K14_10
Internal weights for stage 14, element 10.- See Also:
- Constant Field Values
-
K14_11
private static final double K14_11
Internal weights for stage 14, element 11.- See Also:
- Constant Field Values
-
K14_12
private static final double K14_12
Internal weights for stage 14, element 12.- See Also:
- Constant Field Values
-
K14_13
private static final double K14_13
Internal weights for stage 14, element 13.- See Also:
- Constant Field Values
-
C15
private static final double C15
Time step for stage 15 (interpolation only).- See Also:
- Constant Field Values
-
K15_01
private static final double K15_01
Internal weights for stage 15, element 1.- See Also:
- Constant Field Values
-
K15_06
private static final double K15_06
Internal weights for stage 15, element 6.- See Also:
- Constant Field Values
-
K15_07
private static final double K15_07
Internal weights for stage 15, element 7.- See Also:
- Constant Field Values
-
K15_08
private static final double K15_08
Internal weights for stage 15, element 8.- See Also:
- Constant Field Values
-
K15_09
private static final double K15_09
Internal weights for stage 15, element 9.- See Also:
- Constant Field Values
-
K15_10
private static final double K15_10
Internal weights for stage 15, element 10.- See Also:
- Constant Field Values
-
K15_11
private static final double K15_11
Internal weights for stage 15, element 11.- See Also:
- Constant Field Values
-
K15_12
private static final double K15_12
Internal weights for stage 15, element 12.- See Also:
- Constant Field Values
-
K15_13
private static final double K15_13
Internal weights for stage 15, element 13.- See Also:
- Constant Field Values
-
K15_14
private static final double K15_14
Internal weights for stage 15, element 14.- See Also:
- Constant Field Values
-
C16
private static final double C16
Time step for stage 16 (interpolation only).- See Also:
- Constant Field Values
-
K16_01
private static final double K16_01
Internal weights for stage 16, element 1.- See Also:
- Constant Field Values
-
K16_06
private static final double K16_06
Internal weights for stage 16, element 6.- See Also:
- Constant Field Values
-
K16_07
private static final double K16_07
Internal weights for stage 16, element 7.- See Also:
- Constant Field Values
-
K16_08
private static final double K16_08
Internal weights for stage 16, element 8.- See Also:
- Constant Field Values
-
K16_09
private static final double K16_09
Internal weights for stage 16, element 9.- See Also:
- Constant Field Values
-
K16_10
private static final double K16_10
Internal weights for stage 16, element 10.- See Also:
- Constant Field Values
-
K16_11
private static final double K16_11
Internal weights for stage 16, element 11.- See Also:
- Constant Field Values
-
K16_12
private static final double K16_12
Internal weights for stage 16, element 12.- See Also:
- Constant Field Values
-
K16_13
private static final double K16_13
Internal weights for stage 16, element 13.- See Also:
- Constant Field Values
-
K16_14
private static final double K16_14
Internal weights for stage 16, element 14.- See Also:
- Constant Field Values
-
K16_15
private static final double K16_15
Internal weights for stage 16, element 15.- See Also:
- Constant Field Values
-
D
private static final double[][] D
Interpolation weights. (beware that only the non-null values are in the table)
-
yDotKLast
private double[][] yDotKLast
Last evaluations.
-
v
private double[][] v
Vectors for interpolation.
-
vectorsInitialized
private boolean vectorsInitialized
Initialization indicator for the interpolation vectors.
-
-
Constructor Detail
-
DormandPrince853StepInterpolator
public DormandPrince853StepInterpolator()
Simple constructor. This constructor builds an instance that is not usable yet, thereinitialize(org.apache.commons.math3.ode.AbstractIntegrator, double[], double[][], boolean, org.apache.commons.math3.ode.EquationsMapper, org.apache.commons.math3.ode.EquationsMapper[])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. TheEmbeddedRungeKuttaIntegratoruses the prototyping design pattern to create the step interpolators by cloning an uninitialized model and latter initializing the copy.
-
DormandPrince853StepInterpolator
DormandPrince853StepInterpolator(DormandPrince853StepInterpolator interpolator)
Copy constructor.- Parameters:
interpolator- interpolator to copy from. The copy is a deep copy: its arrays are separated from the original arrays of the instance
-
-
Method Detail
-
doCopy
protected StepInterpolator doCopy()
Really copy the finalized instance.This method is called by
AbstractStepInterpolator.copy()after the step has been finalized. It must perform a deep copy to have an new instance completely independent for the original instance.- Specified by:
doCopyin classAbstractStepInterpolator- Returns:
- a copy of the finalized instance
-
reinitialize
public void reinitialize(AbstractIntegrator integrator, double[] y, double[][] yDotK, boolean forward, EquationsMapper primaryMapper, EquationsMapper[] secondaryMappers)
Reinitialize the instanceSome 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
RungeKuttaIntegratorandEmbeddedRungeKuttaIntegratorabstract classes call this method in order to let the step interpolator perform the evaluations it needs. These evaluations will be performed during the call todoFinalizeif any, i.e. only if the step handler either calls thefinalizeStepmethod or thegetInterpolatedStatemethod (for an interpolator which needs a finalization) or if it clones the step interpolator.- Overrides:
reinitializein classRungeKuttaStepInterpolator- Parameters:
integrator- integrator being usedy- reference to the integrator array holding the state at the end of the stepyDotK- reference to the integrator array holding all the intermediate slopesforward- integration direction indicatorprimaryMapper- equations mapper for the primary equations setsecondaryMappers- equations mappers for the secondary equations sets
-
storeTime
public void storeTime(double t)
Store the current step time.- Overrides:
storeTimein classAbstractStepInterpolator- Parameters:
t- current time
-
computeInterpolatedStateAndDerivatives
protected void computeInterpolatedStateAndDerivatives(double theta, double oneMinusThetaH) throws MaxCountExceededExceptionCompute the state and derivatives at the interpolated time. This is the main processing method that should be implemented by the derived classes to perform the interpolation.- Specified by:
computeInterpolatedStateAndDerivativesin classAbstractStepInterpolator- Parameters:
theta- normalized interpolation abscissa within the step (theta is zero at the previous time step and one at the current time step)oneMinusThetaH- time gap between the interpolated time and the current time- Throws:
MaxCountExceededException- if the number of functions evaluations is exceeded
-
doFinalize
protected void doFinalize() throws MaxCountExceededExceptionReally finalize the step. The default implementation of this method does nothing.- Overrides:
doFinalizein classAbstractStepInterpolator- Throws:
MaxCountExceededException- if the number of functions evaluations is exceeded
-
writeExternal
public void writeExternal(java.io.ObjectOutput out) throws java.io.IOException- Specified by:
writeExternalin interfacejava.io.Externalizable- Overrides:
writeExternalin classRungeKuttaStepInterpolator- Throws:
java.io.IOException
-
readExternal
public void readExternal(java.io.ObjectInput in) throws java.io.IOException, java.lang.ClassNotFoundException- Specified by:
readExternalin interfacejava.io.Externalizable- Overrides:
readExternalin classRungeKuttaStepInterpolator- Throws:
java.io.IOExceptionjava.lang.ClassNotFoundException
-
-