Lau1DCubicSpline.hh

Go to the documentation of this file.

// Copyright University of Warwick 2006 - 2015.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

// Authors:
// Thomas Latham
// John Back
// Paul Harrison

#ifndef LAU_1DCUBICSPLINE
#define LAU_1DCUBICSPLINE

#include <vector>

#include "Rtypes.h"

class Lau1DCubicSpline {

        public:
                enum LauSplineType {
                        StandardSpline, 
                        AkimaSpline, 
                        LinearInterpolation 
                };


                enum LauSplineBoundaryType {
                        Clamped, 
                        Natural, 
                        NotAKnot 
                };


                Lau1DCubicSpline(const std::vector<Double_t>& xs, const std::vector<Double_t>& ys,
                                 LauSplineType type = Lau1DCubicSpline::StandardSpline,
                                 LauSplineBoundaryType leftBound = Lau1DCubicSpline::NotAKnot, 
                                 LauSplineBoundaryType rightBound = Lau1DCubicSpline::NotAKnot,
                                 Double_t dydx0 = 0.0, Double_t dydxn = 0.0);

                virtual ~Lau1DCubicSpline();


                Double_t evaluate(Double_t x) const;


                void updateYValues(const std::vector<Double_t>& ys);


                void updateType(LauSplineType type);


                void updateBoundaryConditions(LauSplineBoundaryType leftBound, 
                                              LauSplineBoundaryType rightBound, 
                                              Double_t dydx0 = 0.0,
                                              Double_t dydxn = 0.0);

        private:
                Lau1DCubicSpline( const Lau1DCubicSpline& rhs );

                Lau1DCubicSpline& operator=(const Lau1DCubicSpline& rhs);

                void init();

                void calcDerivatives();

                void calcDerivativesStandard();
                void calcDerivativesAkima();

                const UInt_t nKnots_;

                std::vector<Double_t> x_;
                std::vector<Double_t> y_;
                std::vector<Double_t> dydx_;

                std::vector<Double_t> a_;
                std::vector<Double_t> b_;
                std::vector<Double_t> c_;
                std::vector<Double_t> d_;

                LauSplineType type_;

                LauSplineBoundaryType leftBound_;
                LauSplineBoundaryType rightBound_;

                Double_t dydx0_;
                Double_t dydxn_;

                ClassDef(Lau1DCubicSpline, 0); // Class for defining a 1D cubic spline
};

#endif