Interpolation.h

Go to the documentation of this file.

#pragma once

#include <array>
#include <vector>
#include <functional>

namespace NuTo
{
namespace Math
{
class Interpolation
{
public:
    Interpolation(std::vector<std::array<double, 2>> data, unsigned numNeighborPoints);

    virtual ~Interpolation() = default;
    virtual double operator()(double x) const = 0;

    virtual double derivative(double x) const = 0;

    std::function<double(double)> f = [this](double x) { return operator()(x); };

    std::function<double(double)> df = [this](double x) { return derivative(x); };

protected:
    std::vector<std::array<double, 2>> mData;

    unsigned mNumNeighborPoints;

    unsigned bisection(double x) const;
};
} // namespace Math
} // namespace NuTo