df.h

#ifndef theplu_yat_utility_multivariable_df
#define theplu_yat_utility_multivariable_df

// $Id: df.h 4252 2022-11-18 02:54:04Z peter $

/*
  Copyright (C) 2022 Peter Johansson

  This file is part of the yat library, http://dev.thep.lu.se/yat

  The yat library is free software; you can redistribute it and/or
  modify it under the terms of the GNU General Public License as
  published by the Free Software Foundation; either version 3 of the
  License, or (at your option) any later version.

  The yat library is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with yat. If not, see <https://www.gnu.org/licenses/>.
*/

#include "f.h"

#include <yat/utility/VectorConstView.h>
#include <yat/utility/VectorView.h>

#include <gsl/gsl_vector.h>

namespace theplu {
namespace yat {
namespace utility {

  namespace multivariable
  {
    template<class FUNC>
    void df(const gsl_vector* gsl_x, void* params, gsl_vector* gradient)
    {
      FUNC* func = static_cast<FUNC*>(params);
      yat::utility::VectorConstView x(gsl_x);
      yat::utility::VectorView g(gradient);
      (*func)(x, g);
    }


    template<class FUNC>
    void fdf(const gsl_vector* x, void* params,
             double* y, gsl_vector* gradient)
    {
      *y = f<FUNC>(x, params);
      df<FUNC>(x, params, gradient);
    }
  } // end namespace multivariable

}}} // end namespaces utility, yat and theplu
#endif