external_data_interpolator.hh

Go to the documentation of this file.

#pragma once
 
// standard library
#include <functional>
#include <string>
#include <vector>
 
// external libraries
#include <boost/math/interpolators/barycentric_rational.hpp>
 
namespace DiFfRG
{
  class ExternalDataInterpolator
  {
  public:
    ExternalDataInterpolator(std::vector<std::string> input_files,
                             std::vector<std::function<double(double)>> post_processors, char separator = ',',
                             bool has_header = false, uint x_column = 0, uint order = 1);
 
    ExternalDataInterpolator(std::vector<std::string> input_files, char separator = ',', bool has_header = false,
                             uint x_column = 0, uint order = 1);
 
    ExternalDataInterpolator();
 
  private:
    void setup(const std::vector<std::string> &input_files,
               const std::vector<std::function<double(double)>> &post_processors, char separator = ',',
               bool has_header = false, uint x_column = 0, uint order = 1);
 
    bool check_consistency(double tolerance) const;
 
  public:
    double value(double x, uint i) const;
 
    double derivative(double x, uint i) const;
 
  private:
    std::string input_file;
    std::vector<std::vector<double>> file_data;
    std::vector<boost::math::interpolators::barycentric_rational<double>> interpolators;
 
    double max_x, min_x;
  };
} // namespace DiFfRG