gslinterpolator.h

#ifndef LSST_GAUSS2D_FIT_GSLINTERPOLATOR_H
#define LSST_GAUSS2D_FIT_GSLINTERPOLATOR_H
 
#ifdef LSST_GAUSS2D_FIT_HAS_GSL
 
#include "gsl.h"
#include "interpolation.h"
 
#include "lsst/gauss2d/object.h"
 
#include <gsl/gsl_errno.h>
#include <gsl/gsl_spline.h>
 
namespace lsst::gauss2d::fit {

class GSLInterpolator : public Object {
public:
    explicit GSLInterpolator(std::vector<double> x, std::vector<double> y,
                             InterpType interp_type = INTERPTYPE_DEFAULT);
    ~GSLInterpolator();
 
    const InterpType get_interp_type() const;
    static constexpr InterpType INTERPTYPE_DEFAULT = InterpType::cspline;

    double get_knot_x(size_t idx) const;
    double get_knot_y(size_t idx) const;

    double eval(double x) const;
    double eval_deriv(double x) const;
    size_t size() const;
 
    std::string repr(bool name_keywords = false,
                     std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
    std::string str() const override;
 
private:
    gsl_interp_accel* _acc;
    gsl_spline* _spline;
 
    const InterpType _interp_type;
    size_t _n_knots;
    std::vector<double> _x;
    std::vector<double> _y;
};
 
}  // namespace lsst::gauss2d::fit
 
#endif  // LSST_GAUSS2D_FIT_HAS_GSL
#endif  // LSST_GAUSS2D_FIT_GSLINTERPOLATOR_H