T
- the type of the field elementspublic abstract class RungeKuttaFieldIntegrator<T extends RealFieldElement<T>> extends AbstractFieldIntegrator<T> implements FieldButcherArrayProvider<T>
These methods are explicit Runge-Kutta methods, their Butcher arrays are as follows :
0 | c2 | a21 c3 | a31 a32 ... | ... cs | as1 as2 ... ass-1 |-------------------------- | b1 b2 ... bs-1 bs
EulerFieldIntegrator
,
ClassicalRungeKuttaFieldIntegrator
,
GillFieldIntegrator
,
MidpointFieldIntegrator
Modifier | Constructor and Description |
---|---|
protected |
RungeKuttaFieldIntegrator(Field<T> field,
String name,
T step)
Simple constructor.
|
Modifier and Type | Method and Description |
---|---|
protected abstract org.apache.commons.math3.ode.nonstiff.RungeKuttaFieldStepInterpolator<T> |
createInterpolator(boolean forward,
T[][] yDotK,
FieldODEStateAndDerivative<T> globalPreviousState,
FieldODEStateAndDerivative<T> globalCurrentState,
FieldEquationsMapper<T> mapper)
Create an interpolator.
|
protected T |
fraction(int p,
int q)
Create a fraction.
|
FieldODEStateAndDerivative<T> |
integrate(FieldExpandableODE<T> equations,
FieldODEState<T> initialState,
T finalTime)
Integrate the differential equations up to the given time.
|
T[] |
singleStep(FirstOrderFieldDifferentialEquations<T> equations,
T t0,
T[] y0,
T t)
Fast computation of a single step of ODE integration.
|
acceptStep, addEventHandler, addEventHandler, addStepHandler, clearEventHandlers, clearStepHandlers, computeDerivatives, getCurrentSignedStepsize, getCurrentStepStart, getEquations, getEvaluations, getEvaluationsCounter, getEventHandlers, getField, getMaxEvaluations, getName, getStepHandlers, getStepSize, getStepStart, initIntegration, isLastStep, resetOccurred, sanityChecks, setIsLastStep, setMaxEvaluations, setStateInitialized, setStepSize, setStepStart
clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait
getA, getB, getC
protected RungeKuttaFieldIntegrator(Field<T> field, String name, T step)
field
- field to which the time and state vector elements belongname
- name of the methodstep
- integration stepprotected T fraction(int p, int q)
p
- numeratorq
- denominatorprotected abstract org.apache.commons.math3.ode.nonstiff.RungeKuttaFieldStepInterpolator<T> createInterpolator(boolean forward, T[][] yDotK, FieldODEStateAndDerivative<T> globalPreviousState, FieldODEStateAndDerivative<T> globalCurrentState, FieldEquationsMapper<T> mapper)
forward
- integration direction indicatoryDotK
- slopes at the intermediate pointsglobalPreviousState
- start of the global stepglobalCurrentState
- end of the global stepmapper
- equations mapper for the all equationspublic FieldODEStateAndDerivative<T> integrate(FieldExpandableODE<T> equations, FieldODEState<T> initialState, T finalTime) throws NumberIsTooSmallException, DimensionMismatchException, MaxCountExceededException, NoBracketingException
This method solves an Initial Value Problem (IVP).
Since this method stores some internal state variables made
available in its public interface during integration (FirstOrderFieldIntegrator.getCurrentSignedStepsize()
), it is not thread-safe.
integrate
in interface FirstOrderFieldIntegrator<T extends RealFieldElement<T>>
equations
- differential equations to integrateinitialState
- initial state (time, primary and secondary state vectors)finalTime
- target time for the integration
(can be set to a value smaller than t0
for backward integration)finalTime
if
integration reached its target, but may be different if some FieldEventHandler
stops it at some point.NumberIsTooSmallException
- if integration step is too smallMaxCountExceededException
- if the number of functions evaluations is exceededNoBracketingException
- if the location of an event cannot be bracketedDimensionMismatchException
public T[] singleStep(FirstOrderFieldDifferentialEquations<T> equations, T t0, T[] y0, T t)
This method is intended for the limited use case of very fast computation of only one step without using any of the rich features of general integrators that may take some time to set up (i.e. no step handlers, no events handlers, no additional states, no interpolators, no error control, no evaluations count, no sanity checks ...). It handles the strict minimum of computation, so it can be embedded in outer loops.
This method is not used at all by the integrate(FieldExpandableODE,
FieldODEState, RealFieldElement)
method. It also completely ignores the step set at
construction time, and uses only a single step to go from t0
to t
.
As this method does not use any of the state-dependent features of the integrator, it should be reasonably thread-safe if and only if the provided differential equations are themselves thread-safe.
equations
- differential equations to integratet0
- initial timey0
- initial value of the state vector at t0t
- target time for the integration
(can be set to a value smaller than t0
for backward integration)t
Copyright © 2003–2016 The Apache Software Foundation. All rights reserved.