kernels.h

Go to the documentation of this file.

#ifndef PURIFY_KERNELS_H
#define PURIFY_KERNELS_H
 
#include "purify/config.h"
#include "purify/types.h"
#include <array>
#include <map>
#include <tuple>
#include <boost/math/special_functions/bessel.hpp>
#include <boost/math/special_functions/sinc.hpp>
 
namespace purify {
 
namespace kernels {
enum class kernel { kb, gauss, box, pswf, kbmin, gauss_alt, kb_presample };
const std::map<std::string, kernel> kernel_from_string = {{"kb", kernel::kb},
                                                          {"gauss", kernel::gauss},
                                                          {"box", kernel::box},
                                                          {"pswf", kernel::pswf},
                                                          {"kbmin", kernel::kbmin},
                                                          {"kb_presample", kernel::kb_presample},
                                                          {"gauss_alt", kernel::gauss_alt}};
 
t_real kaiser_bessel(const t_real x, const t_real J);
t_real kaiser_bessel_general(const t_real x, const t_real J, const t_real alpha);
t_real ft_kaiser_bessel_general(const t_real x, const t_real J, const t_real alpha);
t_real ft_kaiser_bessel(const t_real x, const t_real J);
t_real gaussian(const t_real x, const t_real J);
t_real ft_gaussian(const t_real x, const t_real J);
t_real calc_for_pswf(const t_real x, const t_real J, const t_real alpha);
t_real pswf(const t_real x, const t_real J);
t_real ft_pswf(const t_real x, const t_real J);
std::vector<t_real> kernel_samples(const t_int total_samples,
                                   const std::function<t_real(t_real)> kernelu);
t_real kernel_zero_interp(const std::vector<t_real> &samples, const t_real x, const t_real J);
t_real kernel_linear_interp(const Vector<t_real> &samples, const t_real x, const t_real J);
t_real pill_box(const t_real x, const t_real J);
t_real ft_pill_box(const t_real x, const t_real J);
t_real gaussian_general(const t_real x, const t_real J, const t_real sigma);
t_real ft_gaussian_general(const t_real x, const t_real J, const t_real sigma);
}  // namespace kernels
std::tuple<std::function<t_real(t_real)>, std::function<t_real(t_real)>,
           std::function<t_real(t_real)>, std::function<t_real(t_real)>>
create_kernels(const kernels::kernel kernel_name, const t_uint Ju_, const t_uint Jv_,
               const t_real ftsizeu_, const t_real ftsizev_, const t_real oversample_ratio);
std::tuple<std::function<t_real(t_real)>, std::function<t_real(t_real)>> create_radial_ftkernel(
    const kernels::kernel kernel_name_, const t_uint Ju_, const t_real oversample_ratio);
}  // namespace purify
 
#endif