1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17 package org.orekit.propagation;
18
19 import org.hipparchus.linear.RealMatrix;
20 import org.orekit.attitudes.AttitudeProvider;
21 import org.orekit.errors.OrekitException;
22 import org.orekit.errors.OrekitMessages;
23 import org.orekit.frames.Frame;
24 import org.orekit.propagation.sampling.StepHandlerMultiplexer;
25 import org.orekit.time.AbsoluteDate;
26 import org.orekit.utils.DataDictionary;
27 import org.orekit.utils.DoubleArrayDictionary;
28 import org.orekit.utils.TimeSpanMap;
29
30 import java.util.ArrayList;
31 import java.util.Collections;
32 import java.util.HashMap;
33 import java.util.LinkedList;
34 import java.util.List;
35 import java.util.Map;
36 import java.util.Queue;
37
38
39
40
41
42
43
44
45
46 public abstract class AbstractPropagator implements Propagator {
47
48
49 private final StepHandlerMultiplexer multiplexer;
50
51
52 private AbsoluteDate startDate;
53
54
55 private AttitudeProvider attitudeProvider;
56
57
58 private final List<AdditionalDataProvider<?>> additionalDataProviders;
59
60
61 private final Map<String, TimeSpanMap<Object>> unmanagedStates;
62
63
64 private SpacecraftState initialState;
65
66
67 private AbstractMatricesHarvester harvester;
68
69
70
71 protected AbstractPropagator() {
72 multiplexer = new StepHandlerMultiplexer();
73 additionalDataProviders = new ArrayList<>();
74 unmanagedStates = new HashMap<>();
75 harvester = null;
76 }
77
78
79
80
81 protected void setStartDate(final AbsoluteDate startDate) {
82 this.startDate = startDate;
83 }
84
85
86
87
88 protected AbsoluteDate getStartDate() {
89 return startDate;
90 }
91
92
93 public AttitudeProvider getAttitudeProvider() {
94 return attitudeProvider;
95 }
96
97
98 public void setAttitudeProvider(final AttitudeProvider attitudeProvider) {
99 this.attitudeProvider = attitudeProvider;
100 }
101
102
103 public SpacecraftState getInitialState() {
104 return initialState;
105 }
106
107
108 public Frame getFrame() {
109 return initialState.getFrame();
110 }
111
112
113 public void resetInitialState(final SpacecraftState state) {
114 initialState = state;
115 setStartDate(state.getDate());
116 }
117
118
119 public StepHandlerMultiplexer getMultiplexer() {
120 return multiplexer;
121 }
122
123
124 @Override
125 public void addAdditionalDataProvider(final AdditionalDataProvider<?> provider) {
126
127
128 if (isAdditionalDataManaged(provider.getName())) {
129
130 throw new OrekitException(OrekitMessages.ADDITIONAL_STATE_NAME_ALREADY_IN_USE,
131 provider.getName());
132 }
133
134
135 additionalDataProviders.add(provider);
136
137 }
138
139
140 @Override
141 public List<AdditionalDataProvider<?>> getAdditionalDataProviders() {
142 return Collections.unmodifiableList(additionalDataProviders);
143 }
144
145
146
147
148
149
150 public void removeAdditionalDataProvider(final String name) {
151 additionalDataProviders.removeIf(provider -> provider.getName().equals(name));
152 }
153
154
155 @Override
156 public MatricesHarvester setupMatricesComputation(final String stmName, final RealMatrix initialStm,
157 final DoubleArrayDictionary initialJacobianColumns) {
158 if (stmName == null) {
159 throw new OrekitException(OrekitMessages.NULL_ARGUMENT, "stmName");
160 }
161 harvester = createHarvester(stmName, initialStm, initialJacobianColumns);
162 return harvester;
163 }
164
165
166
167
168
169 public void clearMatricesComputation() {
170 harvester = null;
171 }
172
173
174
175
176
177
178
179
180
181
182
183 protected AbstractMatricesHarvester createHarvester(final String stmName, final RealMatrix initialStm,
184 final DoubleArrayDictionary initialJacobianColumns) {
185
186 throw new UnsupportedOperationException();
187 }
188
189
190
191
192
193 protected AbstractMatricesHarvester getHarvester() {
194 return harvester;
195 }
196
197
198
199
200
201
202 protected SpacecraftState updateUnmanagedData(final SpacecraftState original) {
203
204
205
206 SpacecraftState updated = original;
207
208
209 for (final Map.Entry<String, TimeSpanMap<Object>> entry : unmanagedStates.entrySet()) {
210 updated = updated.addAdditionalData(entry.getKey(),
211 entry.getValue().get(original.getDate()));
212 }
213
214 return updated;
215
216 }
217
218
219
220
221
222
223
224
225 public SpacecraftState updateAdditionalData(final SpacecraftState original) {
226
227
228 SpacecraftState updated = updateUnmanagedData(original);
229
230
231 final Queue<AdditionalDataProvider<?>> pending = new LinkedList<>(getAdditionalDataProviders());
232
233
234 int yieldCount = 0;
235 while (!pending.isEmpty()) {
236 final AdditionalDataProvider<?> provider = pending.remove();
237 if (provider.yields(updated)) {
238
239
240 pending.add(provider);
241 if (++yieldCount >= pending.size()) {
242
243
244
245 break;
246 }
247 } else {
248
249 updated = provider.update(updated);
250 yieldCount = 0;
251 }
252 }
253
254 return updated;
255
256 }
257
258
259
260
261
262
263 protected void initializeAdditionalData(final AbsoluteDate target) {
264 for (final AdditionalDataProvider<?> provider : additionalDataProviders) {
265 provider.init(initialState, target);
266 }
267 }
268
269
270 public boolean isAdditionalDataManaged(final String name) {
271 for (final AdditionalDataProvider<?> provider : additionalDataProviders) {
272 if (provider.getName().equals(name)) {
273 return true;
274 }
275 }
276 return false;
277 }
278
279
280 public String[] getManagedAdditionalData() {
281 final String[] managed = new String[additionalDataProviders.size()];
282 for (int i = 0; i < managed.length; ++i) {
283 managed[i] = additionalDataProviders.get(i).getName();
284 }
285 return managed;
286 }
287
288
289 public SpacecraftState propagate(final AbsoluteDate target) {
290 if (startDate == null) {
291 startDate = getInitialState().getDate();
292 }
293 return propagate(startDate, target);
294 }
295
296
297
298
299 protected void initializePropagation() {
300
301 unmanagedStates.clear();
302
303 if (initialState != null) {
304
305
306
307 for (final DataDictionary.Entry initial : initialState.getAdditionalDataValues().getData()) {
308 if (!isAdditionalDataManaged(initial.getKey())) {
309
310
311 unmanagedStates.put(initial.getKey(), new TimeSpanMap<>(initial.getValue()));
312 }
313 }
314 }
315 }
316
317
318
319
320 protected void stateChanged(final SpacecraftState state) {
321 final AbsoluteDate date = state.getDate();
322 final boolean forward = date.durationFrom(getStartDate()) >= 0.0;
323 for (final DataDictionary.Entry changed : state.getAdditionalDataValues().getData()) {
324 final TimeSpanMap<Object> tsm = unmanagedStates.get(changed.getKey());
325 if (tsm != null) {
326
327 if (forward) {
328 tsm.addValidAfter(changed.getValue(), date, false);
329 } else {
330 tsm.addValidBefore(changed.getValue(), date, false);
331 }
332 }
333 }
334 }
335
336 }