GSLSpline.h

Go to the documentation of this file.

//
// Copyright (c) 2023, Paul Scherrer Institute, Villigen PSI, Switzerland
// All rights reserved
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
 
#ifndef OPALX_GSL_SPLINE_H
#define OPALX_GSL_SPLINE_H
 
#include "CubicSpline.h"
#include "LinearSpline.h"
 
// Compatibility layer for GSL-like interface
 
constexpr int gsl_interp_linear  = 1;
constexpr int gsl_interp_cspline = 0;
 
using gsl_spline = AbstractSpline;
 
using gsl_interp_accel = AbstractSpline::Accelerator;
 
inline gsl_spline* gsl_spline_alloc(const int type, size_t /*size*/) {
    gsl_spline* result{};
    if (type == gsl_interp_cspline) {
        result = new CubicSpline();
    } else if (type == gsl_interp_linear) {
        result = new LinearSpline();
    }
    return result;
}
 
inline gsl_interp_accel* gsl_interp_accel_alloc() { return new AbstractSpline::Accelerator(); }
 
inline void gsl_spline_init(gsl_spline* spline, const double* x, const double* y, const size_t n) {
    spline->init(std::vector<double>{x, x + n}, std::vector<double>{y, y + n});
}
 
inline double gsl_spline_eval(const gsl_spline* spline, const double x, gsl_interp_accel* accel) {
    return spline->eval(x, *accel);
}
 
inline double gsl_spline_eval_integ(
        const gsl_spline* spline, const double xa, const double xb, gsl_interp_accel* accel) {
    return spline->evalIntegral(xa, xb, *accel);
}
 
inline void gsl_spline_free(const gsl_spline* spline) { delete spline; }
 
inline void gsl_interp_accel_free(const gsl_interp_accel* accel) { delete accel; }
 
inline void gsl_interp_accel_reset(gsl_interp_accel* accel) {
    if (accel) {
        accel->reset();
    }
}
 
#endif  // OPALX_GSL_SPLINE_H