GSLInterpolation.h

#ifndef _theplu_yat_regression_gslinterpolation_
#define _theplu_yat_regression_gslinterpolation_

// $Id: GSLInterpolation.h 4207 2022-08-26 04:36:28Z peter $

/*
  Copyright (C) 2004, 2005, 2006, 2007, 2008 Jari Häkkinen, Peter Johansson
  Copyright (C) 2009 Jari Häkkinen
  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 <http://www.gnu.org/licenses/>.
*/

#include <memory>

#include <gsl/gsl_interp.h>

namespace theplu {
namespace yat {
namespace utility {
  class VectorBase;
}
namespace regression {

  class GSLInterpolation
  {

  public:
    GSLInterpolation(const gsl_interp_type*, const utility::VectorBase& x,
                     const utility::VectorBase& y);

    virtual ~GSLInterpolation(void)=0;

    double evaluate(double x);

    double evaluate_derivative(double x);

    double evaluate_derivative2(double x);

    double evaluate_integral(double a, double b);

    double evaluation(void) const;

    unsigned int min_size(void) const;

  private:
    GSLInterpolation(const GSLInterpolation&);
    // no assignment
    GSLInterpolation& operator=(const GSLInterpolation&);

    gsl_interp_accel* accelerator_;
    double evaluation_;
    gsl_interp* interpolator_;
    std::unique_ptr<double[]> x_;
    std::unique_ptr<double[]> y_;
  };

}}} // of namespaces regression, yat, and theplu

#endif