/* 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.util.FastMath;
  import org.junit.jupiter.api.Assertions;
  import org.junit.jupiter.params.ParameterizedTest;
  import org.junit.jupiter.params.provider.ValueSource;
  
  class LogarithmicBarrierCartesianFuelTest {
  
      private static final String ADJOINT_NAME = "adjoint";
  
      @ParameterizedTest
      @ValueSource(doubles = {0., 0.1, 0.5, 0.9})
      void testEvaluatePenaltyFunction(final double norm) {
          // GIVEN
          final LogarithmicBarrierCartesianFuel penalizedCartesianFuel = new LogarithmicBarrierCartesianFuel(ADJOINT_NAME,
                  1., 2., 0.5);
          // WHEN
          final double actualPenalty = penalizedCartesianFuel.evaluatePenaltyFunction(norm);
          // THEN
          Assertions.assertEquals(FastMath.log(norm) + FastMath.log(1 - norm), actualPenalty);
      }
  
      @ParameterizedTest
      @ValueSource(booleans = {false, true})
      void testUpdateFieldAdjointDerivatives(final boolean withMass) {
          // GIVEN
          final double massFlowRateFactor = withMass ? 1 : 0;
          final LogarithmicBarrierCartesianFuel cost = new LogarithmicBarrierCartesianFuel("adjoint", massFlowRateFactor, 2, 0.5);
          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]);
          }
      }
  }