BSpline.h

/* -*- mode: c++; c-basic-offset: 4; -*- */
/************************************************************************
 * Copyright 2022 University Corporation for Atmospheric Research.
 * All rights reserved.
 *
 * Use of this code is subject to the standard BSD license:
 *
 *  http://www.opensource.org/licenses/bsd-license.html
 *
 * See the COPYRIGHT file in the source distribution for the license text,
 * or see this web page:
 *
 *  http://www.eol.ucar.edu/homes/granger/bspline/doc/
 *
 *************************************************************************/
 
#ifndef BSPLINE_H
#define BSPLINE_H
 
#include "BSplineBase.h"
#include <vector>
 
template <class T> struct BSplineP;
 
template <class T> class BSpline : public BSplineBase<T>
{
public:
    BSpline(const T* x, int nx, /* independent variable */
            const T* y,         /* dependent values @ ea X */
            double wl,          /* cutoff wavelength */
            int bc_type = BSplineBase<T>::BC_ZERO_SECOND, int num_nodes = 0);
 
    BSpline(BSplineBase<T>& base, const T* y);
 
    bool solve(const T* y);
 
    T evaluate(T x);
 
    T slope(T x);
 
    T coefficient(int n);
 
    virtual ~BSpline();
 
    using BSplineBase<T>::Debug;
    using BSplineBase<T>::Basis;
    using BSplineBase<T>::DBasis;
 
protected:
    using BSplineBase<T>::OK;
    using BSplineBase<T>::M;
    using BSplineBase<T>::NX;
    using BSplineBase<T>::DX;
    using BSplineBase<T>::base;
    using BSplineBase<T>::xmin;
    using BSplineBase<T>::xmax;
 
    // Our hidden state structure
    BSplineP<T>* s;
    T mean; // Fit without mean and add it in later
};
 
#endif