1 /* Copyright 2002-2013 CS Systèmes d'Information
2 * Licensed to CS Systèmes d'Information (CS) under one or more
3 * contributor license agreements. See the NOTICE file distributed with
4 * this work for additional information regarding copyright ownership.
5 * CS licenses this file to You under the Apache License, Version 2.0
6 * (the "License"); you may not use this file except in compliance with
7 * the License. You may obtain a copy of the License at
8 *
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Unless required by applicable law or agreed to in writing, software
12 * distributed under the License is distributed on an "AS IS" BASIS,
13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 * See the License for the specific language governing permissions and
15 * limitations under the License.
16 */
17 package org.orekit.propagation.events;
18
19 import java.io.Serializable;
20
21 import org.orekit.errors.OrekitException;
22 import org.orekit.propagation.SpacecraftState;
23 import org.orekit.time.AbsoluteDate;
24
25 /** This interface represents space-dynamics aware events detectors.
26 *
27 * <p>It mirrors the {@link org.apache.commons.math3.ode.events.EventHandler
28 * EventHandler} interface from <a href="http://commons.apache.org/math/">
29 * Apache Commons Math</a> but provides a space-dynamics interface to the
30 * methods.</p>
31 *
32 * <p>Events detectors are a useful solution to meet the requirements
33 * of propagators concerning discrete conditions. The state of each
34 * event detector is queried by the integrator at each step. When the
35 * sign of the underlying g switching function changes, the step is rejected
36 * and reduced, in order to make sure the sign changes occur only at steps
37 * boundaries.</p>
38 *
39 * <p>When step ends exactly at a switching function sign change, the corresponding
40 * event is triggered, by calling the {@link #eventOccurred(SpacecraftState, boolean)}
41 * method. The method can do whatever it needs with the event (logging it, performing
42 * some processing, ignore it ...). The return value of the method will be used by
43 * the propagator to stop or resume propagation, possibly changing the state vector.<p>
44 *
45 * @author Luc Maisonobe
46 * @author Véronique Pommier-Maurussane
47 */
48 public interface EventDetector extends Serializable {
49
50 /** Enumerate for actions to be performed when an event occurs.
51 * @deprecated as of 6.1, replaced with {@link org.orekit.propagation.events.handlers.EventHandler.Action}
52 */
53 @Deprecated
54 public enum Action {
55
56 /** Stop indicator.
57 * <p>This value should be used as the return value of the {@link
58 * #eventOccurred eventOccurred} method when the propagation should be
59 * stopped after the event ending the current step.</p>
60 */
61 STOP,
62
63 /** Reset state indicator.
64 * <p>This value should be used as the return value of the {@link
65 * #eventOccurred eventOccurred} method when the propagation should
66 * go on after the event ending the current step, with a new state
67 * (which will be retrieved thanks to the {@link #resetState
68 * resetState} method).</p>
69 */
70 RESET_STATE,
71
72 /** Reset derivatives indicator.
73 * <p>This value should be used as the return value of the {@link
74 * #eventOccurred eventOccurred} method when the propagation should
75 * go on after the event ending the current step, with recomputed
76 * derivatives vector.</p>
77 */
78 RESET_DERIVATIVES,
79
80 /** Continue indicator.
81 * <p>This value should be used as the return value of the {@link
82 * #eventOccurred eventOccurred} method when the propagation should go
83 * on after the event ending the current step.</p>
84 */
85 CONTINUE;
86
87 }
88
89 /** Initialize event handler at the start of a propagation.
90 * <p>
91 * This method is called once at the start of the propagation. It
92 * may be used by the event handler to initialize some internal data
93 * if needed.
94 * </p>
95 * @param s0 initial state
96 * @param t target time for the integration
97 */
98 void init(SpacecraftState s0, AbsoluteDate t);
99
100 /** Compute the value of the switching function.
101 * This function must be continuous (at least in its roots neighborhood),
102 * as the integrator will need to find its roots to locate the events.
103 * @param s the current state information: date, kinematics, attitude
104 * @return value of the switching function
105 * @exception OrekitException if some specific error occurs
106 */
107 double g(SpacecraftState s) throws OrekitException;
108
109 /** Handle an event and choose what to do next.
110
111 * <p>The scheduling between this method and the {@link
112 * org.orekit.propagation.sampling.OrekitStepHandler OrekitStepHandler} method {@link
113 * org.orekit.propagation.sampling.OrekitStepHandler#handleStep(
114 * org.orekit.propagation.sampling.OrekitStepInterpolator, boolean)
115 * handleStep(interpolator, isLast)} is to call this method first and
116 * <code>handleStep</code> afterwards. This scheduling allows the propagator to
117 * pass <code>true</code> as the <code>isLast</code> parameter to the step
118 * handler to make it aware the step will be the last one if this method
119 * returns {@link EventDetector.Action#STOP}. As the interpolator may be used to navigate back
120 * throughout the last step (as {@link
121 * org.orekit.propagation.sampling.OrekitStepNormalizer OrekitStepNormalizer}
122 * does for example), user code called by this method and user
123 * code called by step handlers may experience apparently out of order values
124 * of the independent time variable. As an example, if the same user object
125 * implements both this {@link EventDetector EventDetector} interface and the
126 * {@link org.orekit.propagation.sampling.OrekitFixedStepHandler OrekitFixedStepHandler}
127 * interface, a <em>forward</em> integration may call its
128 * <code>eventOccurred</code> method with a state at 2000-01-01T00:00:10 first
129 * and call its <code>handleStep</code> method with a state at 2000-01-01T00:00:09
130 * afterwards. Such out of order calls are limited to the size of the
131 * integration step for {@link
132 * org.orekit.propagation.sampling.OrekitStepHandler variable step handlers} and
133 * to the size of the fixed step for {@link
134 * org.orekit.propagation.sampling.OrekitFixedStepHandler fixed step handlers}.</p>
135
136 * @param s the current state information : date, kinematics, attitude
137 * @param increasing if true, the value of the switching function increases
138 * when times increases around event (note that increase is measured with respect
139 * to physical time, not with respect to propagation which may go backward in time)
140 * @return one of {@link EventDetector.Action#STOP}, {@link EventDetector.Action#RESET_STATE},
141 * {@link EventDetector.Action#RESET_DERIVATIVES} or {@link EventDetector.Action#CONTINUE}
142 * @exception OrekitException if some specific error occurs
143 * @deprecated as of 6.1 replaced by {@link
144 * org.orekit.propagation.events.handlers.EventHandler#eventOccurred(SpacecraftState,
145 * EventDetector, boolean)}
146 */
147 @Deprecated
148 Action eventOccurred(SpacecraftState s, boolean increasing) throws OrekitException;
149
150 /** Reset the state prior to continue propagation.
151 * <p>This method is called after the step handler has returned and
152 * before the next step is started, but only when {@link
153 * #eventOccurred} has itself returned the {@link EventDetector.Action#RESET_STATE}
154 * indicator. It allows the user to reset the state for the next step,
155 * without perturbing the step handler of the finishing step. If the
156 * {@link #eventOccurred} never returns the {@link EventDetector.Action#RESET_STATE}
157 * indicator, this function will never be called, and it is safe to simply return null.</p>
158 * @param oldState old state
159 * @return new state
160 * @exception OrekitException if the state cannot be reseted
161 * @deprecated as of 6.1 replaced by {@link
162 * org.orekit.propagation.events.handlers.EventHandler#resetState(EventDetector, SpacecraftState)}
163 */
164 @Deprecated
165 SpacecraftState resetState(SpacecraftState oldState) throws OrekitException;
166
167 /** Get the convergence threshold in the event time search.
168 * @return convergence threshold (s)
169 */
170 double getThreshold();
171
172 /** Get maximal time interval between switching function checks.
173 * @return maximal time interval (s) between switching function checks
174 */
175 double getMaxCheckInterval();
176
177 /** Get maximal number of iterations in the event time search.
178 * @return maximal number of iterations in the event time search
179 */
180 int getMaxIterationCount();
181
182 }