LinearSpline.h

Go to the documentation of this file.

//
// Linear Spline Interpolation to replace GSL spline
//
// 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_LINEAR_SPLINE_H
#define OPALX_LINEAR_SPLINE_H
 
#include <vector>
#include "AbstractSpline.h"
 
// \brief Linear spline interpolation class that implements an equivelent of
// the gsl_interp_linear mode of the GSL spline class
class LinearSpline final : public AbstractSpline {
public:
    LinearSpline() = default;
 
    LinearSpline(const std::vector<double>& x, const std::vector<double>& y);
 
    void init(const std::vector<double>& x, const std::vector<double>& y) override;
 
    double eval(double x, Accelerator& accel) const override;
 
    double evalIntegral(double xa, double xb, Accelerator& accel) const override;
 
private:
    void computeCoefficients();
 
    void computeIntegrals();
 
    [[nodiscard]] double integral(size_t i, double dx) const;
 
    [[nodiscard]] double extrapolateLeft(double x) const;
 
    [[nodiscard]] double extrapolateRight(double x) const;
 
    // Members
    std::vector<double> m_{};  // Gradients
    std::vector<double> c_{};  // Constants
    std::vector<double> integrals_{};
};
 
#endif  // OPALX_LINEAR_SPLINE_H