dirac_to_dirac_approx_function_i.h

#ifndef DIRAC_TO_DIRAC_APPROX_FUNCTION_I_H
#define DIRAC_TO_DIRAC_APPROX_FUNCTION_I_H
 
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_matrix_float.h>
#include <gsl/gsl_matrix_long_double.h>
 
#include <type_traits>
 
#include "approximate_options.h"
#include "gsl_minimizer.h"
#include "gsl_vector_matrix_types.h"
 
template <typename T>
class dirac_to_dirac_approx_function_i {
 public:
  using GSLVectorType = typename GSLTemplateTypeAlias<T>::VectorType;
  using GSLVectorViewType = typename GSLTemplateTypeAlias<T>::VectorViewType;
  using GSLMatrixType = typename GSLTemplateTypeAlias<T>::MatrixType;
  using wXf = void (*)(const double* x, double* res, size_t L, size_t N);
  using wXd = wXf;  // (*)(const double* x, double* res, size_t L, size_t N);
 
  virtual ~dirac_to_dirac_approx_function_i() = default;
 
  virtual bool approximate(const T* y, size_t M, size_t L, size_t N, T* x,
                           wXf wXcallback, wXd wXDcallback,
                           GslminimizerResult* result,
                           const ApproximateOptions& options) = 0;
 
  virtual void modified_van_mises_distance_sq(T* distance, const T* y, size_t M,
                                              size_t L, size_t N, size_t bMax,
                                              T* x, wXf wXcallback,
                                              wXd wXDcallback) = 0;
 
  virtual void modified_van_mises_distance_sq_derivative(
      T* gradient, const T* y, size_t M, size_t L, size_t N, size_t bMax, T* x,
      wXf wXcallback, wXd wXDcallback) = 0;
 
  virtual bool approximate(const GSLVectorType* y, size_t L, size_t N,
                           GSLVectorType* x, wXf wXcallback, wXd wXDcallback,
                           GslminimizerResult* result,
                           const ApproximateOptions& options) = 0;
 
  virtual void modified_van_mises_distance_sq(T* distance,
                                              const GSLVectorType* y, size_t L,
                                              size_t N, size_t bMax,
                                              GSLVectorType* x, wXf wXcallback,
                                              wXd wXDcallback) = 0;
 
  virtual void modified_van_mises_distance_sq_derivative(
      GSLMatrixType* gradient, const GSLVectorType* y, size_t L, size_t N,
      size_t bMax, GSLVectorType* x, wXf wXcallback, wXd wXDcallback) = 0;
 
  virtual bool approximate(GSLMatrixType* y, size_t L, GSLMatrixType* x,
                           wXf wXcallback, wXd wXDcallback,
                           GslminimizerResult* result,
                           const ApproximateOptions& options) = 0;
 
  virtual void modified_van_mises_distance_sq(T* distance, GSLMatrixType* y,
                                              size_t L, size_t bMax,
                                              GSLMatrixType* x, wXf wXcallback,
                                              wXd wXDcallback) = 0;
 
  virtual void modified_van_mises_distance_sq_derivative(
      GSLMatrixType* gradient, GSLMatrixType* y, size_t L, size_t bMax,
      GSLMatrixType* x, wXf wXcallback, wXd wXDcallback) = 0;
};
 
#endif  // DIRAC_TO_DIRAC_APPROX_FUNCTION_I_H