/* Copyright 2022-2025 Romain Serra
    * Licensed to CS GROUP (CS) under one or more
    * contributor license agreements.  See the NOTICE file distributed with
    * this work for additional information regarding copyright ownership.
    * CS licenses this file to You under the Apache License, Version 2.0
    * (the "License"); you may not use this file except in compliance with
    * the License.  You may obtain a copy of the License at
    *
    *   http://www.apache.org/licenses/LICENSE-2.0
   *
   * Unless required by applicable law or agreed to in writing, software
   * distributed under the License is distributed on an "AS IS" BASIS,
   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
   * See the License for the specific language governing permissions and
   * limitations under the License.
   */
  package org.orekit.control.indirect.adjoint.cost;
  
  import org.hipparchus.geometry.euclidean.threed.Vector3D;
  import org.junit.jupiter.api.Assertions;
  import org.junit.jupiter.api.Test;
  import org.junit.jupiter.params.ParameterizedTest;
  import org.junit.jupiter.params.provider.ValueSource;
  
  class CartesianFlightDurationCostTest {
  
      private static final String ADJOINT_NAME = "adjoint";
  
      @Test
      void testGetHamiltonianContribution() {
          // GIVEN
          final CartesianFlightDurationCost cost = new CartesianFlightDurationCost(ADJOINT_NAME, 1, 2);
          // WHEN
          final double contribution = cost.getHamiltonianContribution(new double[6], 1);
          // THEN
          Assertions.assertEquals(-1, contribution);
      }
  
      @ParameterizedTest
      @ValueSource(booleans = {false, true})
      void testUpdateFieldAdjointDerivatives(final boolean withMass) {
          // GIVEN
          final double massFlowRateFactor = withMass ? 1 : 0;
          final CartesianFlightDurationCost cost = new CartesianFlightDurationCost(ADJOINT_NAME, massFlowRateFactor, 2);
          final double[] adjoint = new double[withMass ? 7 : 6];
          adjoint[3] = 1;
          final double[] derivatives = new double[adjoint.length];
          // WHEN
          cost.updateAdjointDerivatives(adjoint, 1, derivatives);
          // THEN
          for (int i = 0; i < 6; ++i) {
              Assertions.assertEquals(0., derivatives[i]);
          }
          if (withMass) {
              Assertions.assertNotEquals(0., derivatives[derivatives.length - 1]);
          } else {
              Assertions.assertEquals(0., derivatives[derivatives.length - 1]);
          }
      }
  
      @Test
      void testGetThrustAccelerationVector() {
          // GIVEN
          final CartesianFlightDurationCost cost = new CartesianFlightDurationCost(ADJOINT_NAME, 1, 2);
          // WHEN
          final Vector3D contribution = cost.getThrustAccelerationVector(new double[] {0, 0, 0, 1, 2, 3}, 1);
          // THEN
          Assertions.assertEquals(new Vector3D(1, 2, 3).normalize().scalarMultiply(cost.getMaximumThrustMagnitude()),
                  contribution);
      }
  }