Package org.hipparchus.ode

Class AbstractFieldIntegrator<T extends CalculusFieldElement<T>>

java.lang.Object

org.hipparchus.ode.AbstractFieldIntegrator<T>

Type Parameters:

T - the type of the field elements

All Implemented Interfaces:

FieldODEIntegrator<T>

Direct Known Subclasses:

AdaptiveStepsizeFieldIntegrator, FixedStepRungeKuttaFieldIntegrator

public abstract class AbstractFieldIntegrator<T extends CalculusFieldElement<T>>
extends Object
implements FieldODEIntegrator<T>

Base class managing common boilerplate for all integrators.

Constructor Summary

Constructors

Modifier	Constructor	Description
protected	AbstractFieldIntegrator(Field<T> field, String name)	Build an instance.

Method Summary

Modifier and Type	Method	Description
protected FieldODEStateAndDerivative<T>	acceptStep(AbstractFieldODEStateInterpolator<T> interpolator, T tEnd)	Accept a step, triggering events and step handlers.
void	addEventDetector(FieldODEEventDetector<T> detector)	Add an event detector to the integrator.
void	addStepEndHandler(FieldODEStepEndHandler<T> handler)	Add a handler for step ends to the integrator.
void	addStepHandler(FieldODEStepHandler<T> handler)	Add a step handler to this integrator.
void	clearEventDetectors()	Remove all the event handlers that have been added to the integrator.
void	clearStepEndHandlers()	Remove all the handlers for step ends that have been added to the integrator.
void	clearStepHandlers()	Remove all the step handlers that have been added to the integrator.
T[]	computeDerivatives(T t, T[] y)	Compute the derivatives and check the number of evaluations.
T	getCurrentSignedStepsize()	Get the current signed value of the integration stepsize.
protected FieldExpandableODE<T>	getEquations()	Get the differential equations to integrate.
int	getEvaluations()	Get the number of evaluations of the differential equations function.
protected Incrementor	getEvaluationsCounter()	Get the evaluations counter.
List<FieldODEEventDetector<T>>	getEventDetectors()	Get all the event detectors that have been added to the integrator.
Field<T>	getField()	Get the field to which state vector elements belong.
int	getMaxEvaluations()	Get the maximal number of functions evaluations.
String	getName()	Get the name of the method.
List<FieldODEStepEndHandler<T>>	getStepEndHandlers()	Get all the handlers for step ends that have been added to the integrator.
List<FieldODEStepHandler<T>>	getStepHandlers()	Get all the step handlers that have been added to the integrator.
protected T	getStepSize()	Get the current step size.
FieldODEStateAndDerivative<T>	getStepStart()	Get the state at step start time ti.
protected void	incrementEvaluations(int nTimes)	Increment evaluations of derivatives.
protected FieldODEStateAndDerivative<T>	initIntegration(FieldExpandableODE<T> eqn, FieldODEState<T> s0, T t)	Prepare the start of an integration.
protected boolean	isLastStep()	Check if this step is the last one.
protected boolean	resetOccurred()	Check if a reset occurred while last step was accepted.
protected void	sanityChecks(FieldODEState<T> initialState, T t)	Check the integration span.
protected void	setIsLastStep(boolean isLastStep)	Set the last state flag.
void	setMaxEvaluations(int maxEvaluations)	Set the maximal number of differential equations function evaluations.
protected void	setStateInitialized(boolean stateInitialized)	Set the stateInitialized flag.
protected void	setStepSize(T stepSize)	Set the current step size.
protected void	setStepStart(FieldODEStateAndDerivative<T> stepStart)	Set current step start.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Methods inherited from interface org.hipparchus.ode.FieldODEIntegrator

integrate

Constructor Detail

AbstractFieldIntegrator

protected AbstractFieldIntegrator(Field<T> field,
                                  String name)

Build an instance.

Parameters:

field - field to which the time and state vector elements belong

name - name of the method

Method Detail

getField

public Field<T> getField()

Get the field to which state vector elements belong.

Returns:

field to which state vector elements belong

getName

public String getName()

Get the name of the method.

Specified by:

getName in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

name of the method

addStepHandler

public void addStepHandler(FieldODEStepHandler<T> handler)

Add a step handler to this integrator.

The handler will be called by the integrator for each accepted step.

Specified by:

addStepHandler in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Parameters:

handler - handler for the accepted steps

See Also:

FieldODEIntegrator.getStepHandlers(), FieldODEIntegrator.clearStepHandlers()

getStepHandlers

public List<FieldODEStepHandler<T>> getStepHandlers()

Get all the step handlers that have been added to the integrator.

Specified by:

getStepHandlers in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

an unmodifiable collection of the added events handlers

See Also:

FieldODEIntegrator.addStepHandler(FieldODEStepHandler), FieldODEIntegrator.clearStepHandlers()

clearStepHandlers

public void clearStepHandlers()

Remove all the step handlers that have been added to the integrator.

Specified by:

clearStepHandlers in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

See Also:

FieldODEIntegrator.addStepHandler(FieldODEStepHandler), FieldODEIntegrator.getStepHandlers()

addEventDetector

public void addEventDetector(FieldODEEventDetector<T> detector)

Add an event detector to the integrator.

Specified by:

addEventDetector in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Parameters:

detector - event detector

See Also:

FieldODEIntegrator.getEventDetectors(), FieldODEIntegrator.clearEventDetectors()

getEventDetectors

public List<FieldODEEventDetector<T>> getEventDetectors()

Get all the event detectors that have been added to the integrator.

Specified by:

getEventDetectors in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

an unmodifiable collection of the added events detectors

See Also:

FieldODEIntegrator.addEventDetector(FieldODEEventDetector), FieldODEIntegrator.clearEventDetectors()

clearEventDetectors

public void clearEventDetectors()

Remove all the event handlers that have been added to the integrator.

Specified by:

clearEventDetectors in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

See Also:

FieldODEIntegrator.addEventDetector(FieldODEEventDetector), FieldODEIntegrator.getEventDetectors()

addStepEndHandler

public void addStepEndHandler(FieldODEStepEndHandler<T> handler)

Add a handler for step ends to the integrator.

The stepEndOccurred(state, forward) method of the handler will be called at each step end.

Specified by:

addStepEndHandler in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Parameters:

handler - handler for step ends

See Also:

FieldODEIntegrator.getStepEndHandlers(), FieldODEIntegrator.clearStepEndHandlers()

getStepEndHandlers

public List<FieldODEStepEndHandler<T>> getStepEndHandlers()

Get all the handlers for step ends that have been added to the integrator.

Specified by:

getStepEndHandlers in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

an unmodifiable list of the added step end handlers

See Also:

FieldODEIntegrator.addStepEndHandler(FieldODEStepEndHandler), FieldODEIntegrator.clearStepEndHandlers()

clearStepEndHandlers

public void clearStepEndHandlers()

Remove all the handlers for step ends that have been added to the integrator.

Specified by:

clearStepEndHandlers in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

See Also:

FieldODEIntegrator.addStepEndHandler(FieldODEStepEndHandler), FieldODEIntegrator.getStepEndHandlers()

getCurrentSignedStepsize

public T getCurrentSignedStepsize()

Get the current signed value of the integration stepsize.

This method can be called during integration (typically by the object implementing the differential equations problem) if the signed value of the current stepsize that is tried is needed.

The result is undefined if the method is called outside of calls to integrate.

Specified by:

getCurrentSignedStepsize in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

current signed value of the stepsize

setMaxEvaluations

public void setMaxEvaluations(int maxEvaluations)

Set the maximal number of differential equations function evaluations.

The purpose of this method is to avoid infinite loops which can occur for example when stringent error constraints are set or when lots of discrete events are triggered, thus leading to many rejected steps.

Specified by:

setMaxEvaluations in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Parameters:

maxEvaluations - maximal number of function evaluations (negative values are silently converted to maximal integer value, thus representing almost unlimited evaluations)

getMaxEvaluations

public int getMaxEvaluations()

Get the maximal number of functions evaluations.

Specified by:

getMaxEvaluations in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

maximal number of functions evaluations

getEvaluations

public int getEvaluations()

Get the number of evaluations of the differential equations function.

The number of evaluations corresponds to the last call to the integrate method. It is 0 if the method has not been called yet.

Specified by:

getEvaluations in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

number of evaluations of the differential equations function

initIntegration

protected FieldODEStateAndDerivative<T> initIntegration(FieldExpandableODE<T> eqn,
                                                        FieldODEState<T> s0,
                                                        T t)

Prepare the start of an integration.

Parameters:

eqn - equations to integrate

s0 - initial state vector

t - target time for the integration

Returns:

initial state with derivatives added

getEquations

protected FieldExpandableODE<T> getEquations()

Get the differential equations to integrate.

Returns:

differential equations to integrate

getEvaluationsCounter

protected Incrementor getEvaluationsCounter()

Get the evaluations counter.

Returns:

evaluations counter

computeDerivatives

public T[] computeDerivatives(T t,
                              T[] y)
                       throws MathIllegalArgumentException,
                              MathIllegalStateException,
                              NullPointerException

Compute the derivatives and check the number of evaluations.

Parameters:

t - current value of the independent time variable

y - array containing the current value of the state vector

Returns:

state completed with derivatives

Throws:

MathIllegalArgumentException - if arrays dimensions do not match equations settings

MathIllegalStateException - if the number of functions evaluations is exceeded

NullPointerException - if the ODE equations have not been set (i.e. if this method is called outside of a call to integrate

incrementEvaluations

protected void incrementEvaluations(int nTimes)

Increment evaluations of derivatives.

Parameters:

nTimes - number of evaluations to increment

setStateInitialized

protected void setStateInitialized(boolean stateInitialized)

Set the stateInitialized flag.

This method must be called by integrators with the value false before they start integration, so a proper lazy initialization is done automatically on the first step.

Parameters:

stateInitialized - new value for the flag

acceptStep

protected FieldODEStateAndDerivative<T> acceptStep(AbstractFieldODEStateInterpolator<T> interpolator,
                                                   T tEnd)
                                            throws MathIllegalArgumentException,
                                                   MathIllegalStateException

Accept a step, triggering events and step handlers.

Parameters:

interpolator - step interpolator

tEnd - final integration time

Returns:

state at end of step

Throws:

MathIllegalStateException - if the interpolator throws one because the number of functions evaluations is exceeded

MathIllegalArgumentException - if the location of an event cannot be bracketed

MathIllegalArgumentException - if arrays dimensions do not match equations settings

sanityChecks

protected void sanityChecks(FieldODEState<T> initialState,
                            T t)
                     throws MathIllegalArgumentException

Check the integration span.

Parameters:

initialState - initial state

t - target time for the integration

Throws:

MathIllegalArgumentException - if integration span is too small

MathIllegalArgumentException - if adaptive step size integrators tolerance arrays dimensions are not compatible with equations settings

resetOccurred

protected boolean resetOccurred()

Check if a reset occurred while last step was accepted.

Returns:

true if a reset occurred while last step was accepted

setStepSize

protected void setStepSize(T stepSize)

Set the current step size.

Parameters:

stepSize - step size to set

getStepSize

protected T getStepSize()

Get the current step size.

Returns:

current step size

setStepStart

protected void setStepStart(FieldODEStateAndDerivative<T> stepStart)

Set current step start.

Parameters:

stepStart - step start

getStepStart

public FieldODEStateAndDerivative<T> getStepStart()

Get the state at step start time t_i.

This method can be called during integration (typically by the object implementing the differential equations problem) if the value of the current step that is attempted is needed.

The result is undefined if the method is called outside of calls to integrate.

Specified by:

getStepStart in interface FieldODEIntegrator<T extends CalculusFieldElement<T>>

Returns:

state at step start time t_i

setIsLastStep

protected void setIsLastStep(boolean isLastStep)

Set the last state flag.

Parameters:

isLastStep - if true, this step is the last one

isLastStep

protected boolean isLastStep()

Check if this step is the last one.

Returns:

true if this step is the last one