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 @Override
147 public MatricesHarvester setupMatricesComputation(final String stmName, final RealMatrix initialStm,
148 final DoubleArrayDictionary initialJacobianColumns) {
149 if (stmName == null) {
150 throw new OrekitException(OrekitMessages.NULL_ARGUMENT, "stmName");
151 }
152 harvester = createHarvester(stmName, initialStm, initialJacobianColumns);
153 return harvester;
154 }
155
156
157
158
159
160
161
162
163
164
165
166 protected AbstractMatricesHarvester createHarvester(final String stmName, final RealMatrix initialStm,
167 final DoubleArrayDictionary initialJacobianColumns) {
168
169 throw new UnsupportedOperationException();
170 }
171
172
173
174
175
176 protected AbstractMatricesHarvester getHarvester() {
177 return harvester;
178 }
179
180
181
182
183
184
185 protected SpacecraftState updateUnmanagedData(final SpacecraftState original) {
186
187
188
189 SpacecraftState updated = original;
190
191
192 for (final Map.Entry<String, TimeSpanMap<Object>> entry : unmanagedStates.entrySet()) {
193 updated = updated.addAdditionalData(entry.getKey(),
194 entry.getValue().get(original.getDate()));
195 }
196
197 return updated;
198
199 }
200
201
202
203
204
205
206
207
208 public SpacecraftState updateAdditionalData(final SpacecraftState original) {
209
210
211 SpacecraftState updated = updateUnmanagedData(original);
212
213
214 final Queue<AdditionalDataProvider<?>> pending = new LinkedList<>(getAdditionalDataProviders());
215
216
217 int yieldCount = 0;
218 while (!pending.isEmpty()) {
219 final AdditionalDataProvider<?> provider = pending.remove();
220 if (provider.yields(updated)) {
221
222
223 pending.add(provider);
224 if (++yieldCount >= pending.size()) {
225
226
227
228 break;
229 }
230 } else {
231
232 updated = provider.update(updated);
233 yieldCount = 0;
234 }
235 }
236
237 return updated;
238
239 }
240
241
242
243
244
245
246 protected void initializeAdditionalData(final AbsoluteDate target) {
247 for (final AdditionalDataProvider<?> provider : additionalDataProviders) {
248 provider.init(initialState, target);
249 }
250 }
251
252
253 public boolean isAdditionalDataManaged(final String name) {
254 for (final AdditionalDataProvider<?> provider : additionalDataProviders) {
255 if (provider.getName().equals(name)) {
256 return true;
257 }
258 }
259 return false;
260 }
261
262
263 public String[] getManagedAdditionalData() {
264 final String[] managed = new String[additionalDataProviders.size()];
265 for (int i = 0; i < managed.length; ++i) {
266 managed[i] = additionalDataProviders.get(i).getName();
267 }
268 return managed;
269 }
270
271
272 public SpacecraftState propagate(final AbsoluteDate target) {
273 if (startDate == null) {
274 startDate = getInitialState().getDate();
275 }
276 return propagate(startDate, target);
277 }
278
279
280
281
282 protected void initializePropagation() {
283
284 unmanagedStates.clear();
285
286 if (initialState != null) {
287
288
289
290 for (final DataDictionary.Entry initial : initialState.getAdditionalDataValues().getData()) {
291 if (!isAdditionalDataManaged(initial.getKey())) {
292
293
294 unmanagedStates.put(initial.getKey(), new TimeSpanMap<>(initial.getValue()));
295 }
296 }
297 }
298 }
299
300
301
302
303 protected void stateChanged(final SpacecraftState state) {
304 final AbsoluteDate date = state.getDate();
305 final boolean forward = date.durationFrom(getStartDate()) >= 0.0;
306 for (final DataDictionary.Entry changed : state.getAdditionalDataValues().getData()) {
307 final TimeSpanMap<Object> tsm = unmanagedStates.get(changed.getKey());
308 if (tsm != null) {
309
310 if (forward) {
311 tsm.addValidAfter(changed.getValue(), date, false);
312 } else {
313 tsm.addValidBefore(changed.getValue(), date, false);
314 }
315 }
316 }
317 }
318
319 }