special_generic.h

Go to the documentation of this file.

#ifndef SPECIAL_GENERIC_H
#define SPECIAL_GENERIC_H
 
#include <cmath>
#include "scalar.h"
 
double Ei(double x);
 
double Li2(double x);
 
scalar_t CLi2(scalar_t x);
 
scalar_t CLi3(scalar_t x);
 
scalar_t CLi4(scalar_t x);
 
double Li3(double x);
 
double Cl2(double x);
 
double Cl3(double x);
 
// Bessel
 
double I0(double x);
double I1(double x);
double K0(double x);
double K1(double x);
double K2(double x);
double K3(double x);
double K4(double x);
double Lbessel(double x);
double Mbessel(double x);
double Nbessel(double x);
double K0exp(double x, double z);
double K1exp(double x, double z);
double K2exp(double x, double z);
 
double kron(int x, int y);
 
template <typename T> int sgn(T val) {
    return (T(0) < val) -(val<T(0));
}
 
template <class T>
std::enable_if_t<not std::numeric_limits<T>::is_integer, bool>
fpeq(T x, T y, std::size_t n) {
    const T m = std::min(std::fabs(x), std::fabs(y));
    const int exp = m < std::numeric_limits<T>::min()
                  ? std::numeric_limits<T>::min_exponent - 1
                  : std::ilogb(m);
    return std::fabs(x - y) <= n * std::ldexp(std::numeric_limits<T>::epsilon(), exp);
}
 
#endif // SPECIAL_GENERIC_H