interpolate1d.hpp

Go to the documentation of this file.

/* Copyright (c) 2023 Otto Link. Distributed under the terms of the GNU General
   Public License. The full license is in the file LICENSE, distributed with
   this software. */
 
#pragma once
#include <map>
#include <vector>
 
#include <gsl/gsl_interp.h>
#include <gsl/gsl_spline.h>
 
namespace hmap
{
 
enum InterpolationMethod1D : int
{
  AKIMA,          
  AKIMA_PERIODIC, 
  CUBIC,          
  CUBIC_PERIODIC, 
  LINEAR,         
  POLYNOMIAL,     
  STEFFEN,        
};
 
class Interpolator1D
{
public:
  Interpolator1D(const std::vector<float> &x,
                 const std::vector<float> &y,
                 InterpolationMethod1D method = InterpolationMethod1D::LINEAR);
 
  ~Interpolator1D();
 
  float operator()(float x) const;
 
  float interpolate(float x) const;
 
private:
  gsl_spline         *interp; 
  gsl_interp_accel   *accel_; 
  std::vector<double> x_data; 
  std::vector<double> y_data; 
  double              xmin;   
  double              xmax;   
};
 
} // namespace hmap