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.propagation.semianalytical.dsst.utilities;
18  
19  import org.hipparchus.Field;
20  import org.hipparchus.CalculusFieldElement;
21  import org.hipparchus.util.FastMath;
22  import org.hipparchus.util.MathArrays;
23  
24  /** Interpolation grid where points obey a maximum time gap.
25   * <p>
26   * The grid is adapted to the step considered,
27   * meaning that for short steps, the grid will have numerous points.
28   * </p>
29   *
30   * @author Luc Maisonobe
31   * @since 7.1
32   * @param <T> type of the field elements
33   */
34  public class FieldMaxGapInterpolationGrid <T extends CalculusFieldElement<T>> implements FieldInterpolationGrid<T> {
35  
36      /** Maximum time gap. */
37      private final T maxGap;
38  
39      /** Field used by default. */
40      private final Field<T> field;
41  
42      /** Constructor.
43       * @param field field used by default
44       * @param maxGap maximum time gap between interpolation points
45       */
46      public FieldMaxGapInterpolationGrid(final Field<T> field, final T maxGap) {
47          this.field  = field;
48          this.maxGap = maxGap;
49      }
50  
51      /** {@inheritDoc} */
52      @Override
53      public T[] getGridPoints(final T stepStart, final T stepEnd) {
54          final int pointsPerStep = FastMath.max(2, (int) FastMath.ceil(FastMath.abs(stepEnd.getReal() - stepStart.getReal()) / maxGap.getReal()));;
55          final T[] grid = MathArrays.buildArray(field, pointsPerStep);
56  
57          final T stepSize = stepEnd.subtract(stepStart).divide(pointsPerStep - 1.);
58          for (int i = 0; i < pointsPerStep; i++) {
59              grid[i] = stepSize.multiply(i).add(stepStart);
60          }
61  
62          return grid;
63      }
64  }