public interface StepInterpolator extends Externalizable
The various ODE integrators provide objects implementing this interface to the step handlers. These objects are often custom objects tightly bound to the integrator internal algorithms. The handlers can use these objects to retrieve the state vector at intermediate times between the previous and the current grid points (this feature is often called dense output).
One important thing to note is that the step handlers may be so tightly bound to the integrators that they often
share some internal state arrays. This imply that one should never use a direct reference to a step
interpolator outside of the step handler, either for future use or for use in another thread. If such a need arise,
the step interpolator must be copied using the dedicated copy()
method.
FirstOrderIntegrator
,
SecondOrderIntegrator
,
StepHandler
Modifier and Type | Method and Description |
---|---|
StepInterpolator |
copy()
Copy the instance.
|
double |
getCurrentTime()
Get the current grid point time.
|
double[] |
getInterpolatedDerivatives()
Get the derivatives of the state vector of the interpolated point.
|
double[] |
getInterpolatedSecondaryDerivatives(int index)
Get the interpolated secondary derivatives corresponding to the secondary equations.
|
double[] |
getInterpolatedSecondaryState(int index)
Get the interpolated secondary state corresponding to the secondary equations.
|
double[] |
getInterpolatedState()
Get the state vector of the interpolated point.
|
double |
getInterpolatedTime()
Get the time of the interpolated point.
|
double |
getPreviousTime()
Get the previous grid point time.
|
boolean |
isForward()
Check if the natural integration direction is forward.
|
void |
setInterpolatedTime(double time)
Set the time of the interpolated point.
|
readExternal, writeExternal
double getPreviousTime()
double getCurrentTime()
double getInterpolatedTime()
setInterpolatedTime(double)
has not been called, it returns
the current grid point time.void setInterpolatedTime(double time)
Setting the time outside of the current step is now allowed, but should be used with care since the accuracy of the interpolator will probably be very poor far from this step. This allowance has been added to simplify implementation of search algorithms near the step endpoints.
Setting the time changes the instance internal state. If a specific state must be preserved, a copy of the
instance must be created using copy()
.
time
- time of the interpolated pointdouble[] getInterpolatedState()
The returned vector is a reference to a reused array, so it should not be modified and it should be copied if it needs to be preserved across several calls.
getInterpolatedTime()
MaxCountExceededException
- if the number of functions evaluations is exceededgetInterpolatedDerivatives()
double[] getInterpolatedDerivatives()
The returned vector is a reference to a reused array, so it should not be modified and it should be copied if it needs to be preserved across several calls.
getInterpolatedTime()
MaxCountExceededException
- if the number of functions evaluations is exceededgetInterpolatedState()
double[] getInterpolatedSecondaryState(int index)
The returned vector is a reference to a reused array, so it should not be modified and it should be copied if it needs to be preserved across several calls.
index
- index of the secondary set, as returned by ExpandableStatefulODE.addSecondaryEquations()MaxCountExceededException
- if the number of functions evaluations is exceededgetInterpolatedState()
,
getInterpolatedDerivatives()
,
getInterpolatedSecondaryDerivatives(int)
,
setInterpolatedTime(double)
double[] getInterpolatedSecondaryDerivatives(int index)
The returned vector is a reference to a reused array, so it should not be modified and it should be copied if it needs to be preserved across several calls.
index
- index of the secondary set, as returned by ExpandableStatefulODE.addSecondaryEquations()MaxCountExceededException
- if the number of functions evaluations is exceededgetInterpolatedState()
,
getInterpolatedDerivatives()
,
getInterpolatedSecondaryState(int)
,
setInterpolatedTime(double)
boolean isForward()
This method provides the integration direction as specified by the integrator itself, it avoid some nasty problems in degenerated cases like null steps due to cancellation at step initialization, step control or discrete events triggering.
StepInterpolator copy()
The copied instance is guaranteed to be independent from the original one. Both can be used with different settings for interpolated time without any side effect.
MaxCountExceededException
- if the number of functions evaluations is exceeded
during step finalizationsetInterpolatedTime(double)
Copyright © 2023 CNES. All rights reserved.