interpolation.hh

Go to the documentation of this file.

#pragma once
 
#include <DiFfRG/common/kokkos.hh>
#include <DiFfRG/common/types.hh>
 
namespace DiFfRG
{
  template <typename T>
  concept has_interpolator_types = requires {
    typename T::value_type;
    typename T::ctype;
  };
 
  template <typename T>
  concept has_interpolator_methods = has_interpolator_types<T> && requires(T t) {
    t.get_coordinates();
    t.template get_on<CPU_memory>();
    t.template get_on<GPU_memory>();
    requires has_n_call_operator<T, typename T::ctype, T::dim>;
  };
 
  template <typename T>
  concept is_interpolator = has_interpolator_types<T> && has_interpolator_methods<T>;
} // namespace DiFfRG
 
#include <DiFfRG/discretization/coordinates/combined_coordinates.hh>
#include <DiFfRG/discretization/coordinates/coordinates.hh>
#include <DiFfRG/discretization/coordinates/stack_coordinates.hh>
 
#include <DiFfRG/physics/interpolation/linear_interpolator.hh>
 
#include <DiFfRG/physics/interpolation/spline_interpolator_1d.hh>
#include <DiFfRG/physics/interpolation/spline_interpolator_1d_stack.hh>