PCHIPInterpolator.h

#ifndef included_PCHIPInterpolator_h
#define included_PCHIPInterpolator_h
 
#include <vector>
#include <string>
#include <memory>
 
namespace CAROM {
 
class Matrix;
class Vector;
 
class PCHIPInterpolator
{
public:
 
 
    PCHIPInterpolator()
    {}
    ~PCHIPInterpolator()
    {}
 
    void interpolate(std::vector<Vector>& snapshot_ts,
                     std::vector<std::shared_ptr<Vector>>& snapshots,
                     std::vector<Vector>& output_ts,
                     std::vector<std::shared_ptr<Vector>>& output_snapshots);
 
    void interpolate(std::vector<Vector>& snapshot_ts,
                     std::vector<std::shared_ptr<Vector>>& snapshots,
                     int n_out,
                     std::vector<Vector>& output_ts,
                     std::vector<std::shared_ptr<Vector>>& output_snapshots);
 
private:
 
    double computeDerivative(double S1, double S2, double h1, double h2) const;
    double computeH1(double x, double xl, double xr) const;
    double computeH2(double x, double xl, double xr) const;
    double computeH3(double x, double xl, double xr) const;
    double computeH4(double x, double xl, double xr) const;
    double computePhi(double t) const;
    double computePsi(double t) const;
    int sign(double a) const;
};
 
}
#endif