dirac_to_dirac_approx_i.h

#ifndef DIRAC_TO_DIRAC_APPROX_I_H
#define DIRAC_TO_DIRAC_APPROX_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_i {
 public:
  using GSLVectorType = typename GSLTemplateTypeAlias<T>::VectorType;
  using GSLVectorViewType = typename GSLTemplateTypeAlias<T>::VectorViewType;
  using GSLMatrixType = typename GSLTemplateTypeAlias<T>::MatrixType;
  using GSLMatrixViewType = typename GSLTemplateTypeAlias<T>::MatrixViewType;
 
  virtual ~dirac_to_dirac_approx_i() = default;
 
  virtual bool approximate(const T* y, size_t M, size_t L, size_t N,
                           size_t bMax, T* x, const T* wX, const T* wY,
                           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, const T* wX,
                                              const T* wY) = 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,
      const T* wX, const T* wY) = 0;
 
  virtual bool approximate(const GSLVectorType* y, size_t L, size_t N,
                           size_t bMax, GSLVectorType* x,
                           const GSLVectorType* wX, const GSLVectorType* wY,
                           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, const GSLVectorType* wX, const GSLVectorType* wY) = 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, const GSLVectorType* wX,
      const GSLVectorType* wY) = 0;
 
  virtual bool approximate(GSLMatrixType* y, size_t L, size_t bMax,
                           GSLMatrixType* x, const GSLVectorType* wX,
                           const GSLVectorType* wY, 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,
                                              const GSLVectorType* wX,
                                              const GSLVectorType* wY) = 0;
 
  virtual void modified_van_mises_distance_sq_derivative(
      GSLMatrixType* gradient, GSLMatrixType* y, size_t L, size_t bMax,
      GSLMatrixType* x, const GSLVectorType* wX, const GSLVectorType* wY) = 0;
};
 
#endif  // DIRAC_TO_DIRAC_APPROX_I_H