DeBoor.H

Go to the documentation of this file.

/* Copyright 2000-2018 Matt Flax <flatmax@flatmax.org>
   This file is part of GTK+ IOStream class set

   GTK+ IOStream 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 2 of the License, or
   (at your option) any later version.

   GTK+ IOStream is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You have received a copy of the GNU General Public License
   along with GTK+ IOStream
 */
#ifndef DEBOOR_H
#define DEBOOR_H

#include <vector>
//#include <array>
using namespace std;

class DeBoor {
    float breakp[5]; 
    float coef[4][4]; 
    int l; 

    float **v; 
    float **a; 
    int N; 

//    vector<array<float, 7>> v; ///< work array
//    vector<array<float, 4>> a; ///< results array

    void createMatrices(int n);

    float **&createMatrix(int k);
    void deleteMatrices(); 

    void deleteMatrix(float **m);

    vector<float> resultCSAPS;
public:
    DeBoor();
    virtual ~DeBoor();

    void csaps(float *x, float *y, float *dy, int n, float s);

    float operator[](int i);
};

#endif // DEBOOR_H