rbf_interpolation.h

#ifndef __RBF_INTERPOLATION__
#define __RBF_INTERPOLATION__

#include <eigen/Eigen/Dense>
#include "ROL_OptimizationProblem.hpp"
#include "ROL_StdVector.hpp"

namespace PHiLiP {
namespace ProperOrthogonalDecomposition {
using Eigen::MatrixXd;
using Eigen::VectorXd;
using Eigen::RowVectorXd;

class RBFInterpolation: public ROL::Objective<double>
{
public:
    RBFInterpolation(const MatrixXd& data_coordinates, const VectorXd& data_values, std::string kernel);

    void computeWeights();

    double radialBasisFunction(double r) const;

    double evaluate(const RowVectorXd& evaluate_coordinate) const;

    VectorXd weights;

    const MatrixXd data_coordinates;

    const VectorXd data_values;

    const std::string kernel;

    typedef std::vector<double> vector;
    typedef ROL::Vector<double>      V;

    template<class VectorType>
    ROL::Ptr<const vector> getVector( const V& x ) {
        return dynamic_cast<const VectorType&>((x)).getVector();
    }

    double value(const ROL::Vector<double> &x, double &/*tol*/ );

};

}
}


#endif