1 /* Copyright 2002-2025 CS GROUP
2 * Licensed to CS GROUP (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.estimation.leastsquares;
18
19 import org.orekit.estimation.measurements.ObservedMeasurement;
20 import org.orekit.orbits.Orbit;
21 import org.orekit.propagation.MatricesHarvester;
22 import org.orekit.propagation.PropagationType;
23 import org.orekit.propagation.Propagator;
24 import org.orekit.propagation.SpacecraftState;
25 import org.orekit.propagation.conversion.PropagatorBuilder;
26 import org.orekit.propagation.semianalytical.dsst.DSSTPropagator;
27 import org.orekit.utils.ParameterDriversList;
28
29 import java.util.List;
30
31 /** Bridge between {@link ObservedMeasurement measurements} and {@link
32 * org.hipparchus.optim.nonlinear.vector.leastsquares.LeastSquaresProblem
33 * least squares problems}.
34 * <p>
35 * This class is an adaption of the {@link BatchLSModel} class
36 * for the {@link DSSTPropagator DSST propagator}.
37 * </p>
38 * @author Luc Maisonobe
39 * @author Bryan Cazabonne
40 * @since 10.0
41 *
42 */
43 public class DSSTBatchLSModel extends AbstractBatchLSModel {
44
45 /** Name of the State Transition Matrix state. */
46 private static final String STM_NAME = DSSTBatchLSModel.class.getName() + "-derivatives";
47
48 /** Type of the orbit used for the propagation.*/
49 private PropagationType propagationType;
50
51 /** Simple constructor.
52 * @param propagatorBuilders builders to use for propagation
53 * @param measurements measurements
54 * @param estimatedMeasurementsParameters estimated measurements parameters
55 * @param observer observer to be notified at model calls
56 * @param propagationType type of the orbit used for the propagation (mean or osculating)
57 */
58 public DSSTBatchLSModel(final PropagatorBuilder[] propagatorBuilders,
59 final List<ObservedMeasurement<?>> measurements,
60 final ParameterDriversList estimatedMeasurementsParameters,
61 final ModelObserver observer,
62 final PropagationType propagationType) {
63 // call super constructor
64 super(propagatorBuilders, measurements, estimatedMeasurementsParameters, observer);
65 this.propagationType = propagationType;
66 }
67
68 /** {@inheritDoc} */
69 @Override
70 protected MatricesHarvester configureHarvester(final Propagator propagator) {
71 return propagator.setupMatricesComputation(STM_NAME, null, null);
72 }
73
74 /** {@inheritDoc} */
75 @Override
76 protected Orbit configureOrbits(final MatricesHarvester harvester, final Propagator propagator) {
77 // Cast
78 final DSSTPropagator dsstPropagator = (DSSTPropagator) propagator;
79 // Mean orbit
80 final SpacecraftState initial = dsstPropagator.initialIsOsculating() ?
81 DSSTPropagator.computeMeanState(dsstPropagator.getInitialState(), dsstPropagator.getAttitudeProvider(), dsstPropagator.getAllForceModels()) :
82 dsstPropagator.getInitialState();
83 // Compute short period derivatives at the beginning of the iteration
84 harvester.setReferenceState(initial);
85 return initial.getOrbit();
86 }
87
88 }