Class EventState
java.lang.Object
org.apache.commons.math3.ode.events.EventState
This class handles the state for one
event handler during integration steps.
Each time the integrator proposes a step, the event handler switching function should be checked. This class handles the state of one handler during one integration step, with references to the state at the end of the preceding step. This information is used to decide if the handler should trigger an event or not during the proposed step.
- Since:
- 1.2
-
Nested Class Summary
Nested ClassesModifier and TypeClassDescriptionprivate static classLocal wrapper to propagate exceptions. -
Field Summary
FieldsModifier and TypeFieldDescriptionprivate final doubleConvergence threshold for event localization.private ExpandableStatefulODEEquation being integrated.private booleanIntegration direction.private doubleValue of the events handler at the beginning of the step.private booleanSimulated sign of g0 (we cheat when crossing events).private final EventHandlerEvent handler.private booleanVariation direction around pending event.private final doubleMaximal time interval between events handler checks.private final intUpper limit in the iteration count for event localization.private EventHandler.ActionNext action indicator.private booleanIndicator of event expected during the step.private doubleOccurrence time of the pending event.private doubleOccurrence time of the previous event.private final UnivariateSolverRoot-finding algorithm to use to detect state events.private doubleTime at the beginning of the step. -
Constructor Summary
ConstructorsConstructorDescriptionEventState(EventHandler handler, double maxCheckInterval, double convergence, int maxIterationCount, UnivariateSolver solver) Simple constructor. -
Method Summary
Modifier and TypeMethodDescriptionbooleanevaluateStep(StepInterpolator interpolator) Evaluate the impact of the proposed step on the event handler.private double[]getCompleteState(StepInterpolator interpolator) Get the complete state (primary and secondary).doubleGet the convergence threshold for event localization.Get the underlying event handler.doubleGet the occurrence time of the event triggered in the current step.doubleGet the maximal time interval between events handler checks.intGet the upper limit in the iteration count for event localization.voidreinitializeBegin(StepInterpolator interpolator) Reinitialize the beginning of the step.booleanreset(double t, double[] y) Let the event handler reset the state if it wants.voidsetExpandable(ExpandableStatefulODE expandable) Set the equation.voidstepAccepted(double t, double[] y) Acknowledge the fact the step has been accepted by the integrator.booleanstop()Check if the integration should be stopped at the end of the current step.
-
Field Details
-
handler
Event handler. -
maxCheckInterval
private final double maxCheckIntervalMaximal time interval between events handler checks. -
convergence
private final double convergenceConvergence threshold for event localization. -
maxIterationCount
private final int maxIterationCountUpper limit in the iteration count for event localization. -
expandable
Equation being integrated. -
t0
private double t0Time at the beginning of the step. -
g0
private double g0Value of the events handler at the beginning of the step. -
g0Positive
private boolean g0PositiveSimulated sign of g0 (we cheat when crossing events). -
pendingEvent
private boolean pendingEventIndicator of event expected during the step. -
pendingEventTime
private double pendingEventTimeOccurrence time of the pending event. -
previousEventTime
private double previousEventTimeOccurrence time of the previous event. -
forward
private boolean forwardIntegration direction. -
increasing
private boolean increasingVariation direction around pending event. (this is considered with respect to the integration direction) -
nextAction
Next action indicator. -
solver
Root-finding algorithm to use to detect state events.
-
-
Constructor Details
-
EventState
public EventState(EventHandler handler, double maxCheckInterval, double convergence, int maxIterationCount, UnivariateSolver solver) Simple constructor.- Parameters:
handler- event handlermaxCheckInterval- maximal time interval between switching function checks (this interval prevents missing sign changes in case the integration steps becomes very large)convergence- convergence threshold in the event time searchmaxIterationCount- upper limit of the iteration count in the event time searchsolver- Root-finding algorithm to use to detect state events
-
-
Method Details
-
getEventHandler
Get the underlying event handler.- Returns:
- underlying event handler
-
setExpandable
Set the equation.- Parameters:
expandable- equation being integrated
-
getMaxCheckInterval
public double getMaxCheckInterval()Get the maximal time interval between events handler checks.- Returns:
- maximal time interval between events handler checks
-
getConvergence
public double getConvergence()Get the convergence threshold for event localization.- Returns:
- convergence threshold for event localization
-
getMaxIterationCount
public int getMaxIterationCount()Get the upper limit in the iteration count for event localization.- Returns:
- upper limit in the iteration count for event localization
-
reinitializeBegin
Reinitialize the beginning of the step.- Parameters:
interpolator- valid for the current step- Throws:
MaxCountExceededException- if the interpolator throws one because the number of functions evaluations is exceeded
-
getCompleteState
Get the complete state (primary and secondary).- Parameters:
interpolator- interpolator to use- Returns:
- complete state
-
evaluateStep
public boolean evaluateStep(StepInterpolator interpolator) throws MaxCountExceededException, NoBracketingException Evaluate the impact of the proposed step on the event handler.- Parameters:
interpolator- step interpolator for the proposed step- Returns:
- true if the event handler triggers an event before the end of the proposed step
- Throws:
MaxCountExceededException- if the interpolator throws one because the number of functions evaluations is exceededNoBracketingException- if the event cannot be bracketed
-
getEventTime
public double getEventTime()Get the occurrence time of the event triggered in the current step.- Returns:
- occurrence time of the event triggered in the current step or infinity if no events are triggered
-
stepAccepted
public void stepAccepted(double t, double[] y) Acknowledge the fact the step has been accepted by the integrator.- Parameters:
t- value of the independent time variable at the end of the stepy- array containing the current value of the state vector at the end of the step
-
stop
public boolean stop()Check if the integration should be stopped at the end of the current step.- Returns:
- true if the integration should be stopped
-
reset
public boolean reset(double t, double[] y) Let the event handler reset the state if it wants.- Parameters:
t- value of the independent time variable at the beginning of the next stepy- array were to put the desired state vector at the beginning of the next step- Returns:
- true if the integrator should reset the derivatives too
-