ell_calc_core.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <cstddef>  // For size_t
#include <tuple>
 
class EllCalcCore {
  private:
    double _n_f;
    double _n_plus_1;
    double _half_n;
    double _inv_n;
    double _n_sq;
    double _cst1;
    double _cst2;
 
  public:
    explicit EllCalcCore(const size_t ndim) noexcept
        : _n_f{static_cast<double>(ndim)},
          _n_plus_1{_n_f + 1.0},
          _half_n{_n_f / 2.0},
          _inv_n{1.0 / _n_f},
          _n_sq{_n_f * _n_f},
          _cst1{_n_sq / (_n_sq - 1.0)},
          _cst2{2.0 / _n_plus_1} {}
 
    // Rule of five: explicitly default destructor and assignment operators
    ~EllCalcCore() noexcept = default;
    EllCalcCore& operator=(const EllCalcCore&) noexcept = default;
    EllCalcCore& operator=(EllCalcCore&&) noexcept = default;
 
    EllCalcCore(EllCalcCore&& E) noexcept = default;
 
    EllCalcCore(const EllCalcCore& E) noexcept = default;
 
    auto calc_parallel_cut(const double beta0, const double beta1, const double tsq) const noexcept
        -> std::tuple<double, double, double> {
        const auto b0b1 = beta0 * beta1;
        const auto eta = tsq + this->_n_f * b0b1;
        return this->calc_parallel_cut_fast(beta0, beta1, tsq, b0b1, eta);
    }
 
    auto calc_parallel_cut_fast_old(const double beta0, const double beta1, const double tsq,
                                    const double b0b1, const double eta) const noexcept
        -> std::tuple<double, double, double>;
 
    auto calc_parallel_cut_fast(const double beta0, const double beta1, const double tsq,
                                const double b0b1, const double eta) const noexcept
        -> std::tuple<double, double, double>;
 
    auto calc_parallel_central_cut(const double beta1, const double tsq) const noexcept
        -> std::tuple<double, double, double>;
 
    auto calc_bias_cut(const double beta, const double tau) const noexcept
        -> std::tuple<double, double, double> {
        return this->calc_bias_cut_fast(beta, tau, tau + this->_n_f * beta);
    }
 
    auto calc_bias_cut_fast(const double beta, const double tau, const double eta) const noexcept
        -> std::tuple<double, double, double>;
 
    auto calc_central_cut(const double tau) const noexcept -> std::tuple<double, double, double>;
};  // } EllCalcCore