BSpline.h

/* -*- mode: c++; c-basic-offset: 4; -*- */
/************************************************************************
 * Copyright 2009 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