ell_calc.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <cassert>
 
#include "ell_calc_core.hpp"
#include "ell_config.hpp"
 
class EllCalc {
  public:
    bool use_parallel_cut = true;
 
  protected:
    double _n_f;
    EllCalcCore _helper;
 
  public:
    explicit EllCalc(const size_t ndim) : _n_f{static_cast<double>(ndim)}, _helper{ndim} {
        assert(ndim >= 2U);  // do not accept one-dimensional
    }
 
    EllCalc(EllCalc&& E) = default;
 
    EllCalc& operator=(const EllCalc& other) = default;
 
    EllCalc& operator=(EllCalc&& other) = default;
 
    ~EllCalc() = default;
 
    EllCalc(const EllCalc& E) = default;
 
    auto calc_parallel_bias_cut(double beta0, double beta1, double tsq) const -> CutResult;
 
    auto calc_parallel_central_cut(double beta1, double tsq) const -> CutResult;
 
    auto calc_bias_cut(double beta, double tsq) const -> CutResult;
 
    auto calc_central_cut(double tsq) const -> CutResult;
 
    auto calc_parallel_bias_cut_q(double beta0, double beta1, double tsq) const -> CutResult;
 
    auto calc_bias_cut_q(double beta, double tsq) const -> CutResult;
};  // } EllCalc