#ifndef MANDEL_TYPES_H #define MANDEL_TYPES_H #include #include #include "Fixed.h" #ifdef WITH_BOOST # include # if defined(__GNUC__) || defined(__INTEL_COMPILER) # include # endif # include #endif namespace mnd { #ifdef WITH_BOOST #if 0//defined(__GNUC__) || defined(__INTEL_COMPILER) using Float128 = boost::multiprecision::float128; #else //using Float128 = boost::multiprecision::cpp_bin_float_quad; using Float128 = boost::multiprecision::number< boost::multiprecision::backends::cpp_bin_float< 112, boost::multiprecision::backends::digit_base_2, void, boost::int16_t, -16382, 16383>, boost::multiprecision::et_off>; #endif inline Float128 abs(const Float128& x) { return boost::multiprecision::abs(x); } inline Float128 floor(const Float128& x) { return boost::multiprecision::floor(x); } inline Float128 log(const Float128& x) { return boost::multiprecision::log(x); } inline Float128 log2(const Float128& x) { return boost::multiprecision::log2(x); } inline Float128 pow(const Float128& x, const Float128& y) { return boost::multiprecision::pow(x, y); } /* using Float256 = boost::multiprecision::number< boost::multiprecision::backends::cpp_bin_float< 240, boost::multiprecision::backends::digit_base_2, void, boost::int16_t, -16382, 16383>, boost::multiprecision::et_off>; */ using Float256 = long double; /*inline Float256 abs(const Float256& x) { return boost::multiprecision::abs(x); } inline Float256 floor(const Float256& x) { return boost::multiprecision::floor(x); } inline Float256 log(const Float256& x) { return boost::multiprecision::log(x); } inline Float256 log2(const Float256& x) { return boost::multiprecision::log2(x); } inline Float256 pow(const Float256& x, const Float256& y) { return boost::multiprecision::pow(x, y); } */ using Real = Float128; using Integer = boost::multiprecision::int256_t; #else using Real = double; using Integer = int64_t; #endif inline double abs(double x) { return ::abs(x); } inline float abs(float x) { return ::abs(x); } inline double floor(double x) { return ::floor(x); } inline float floor(float x) { return ::floorf(x); } inline double log(double x) { return ::log(x); } inline float log(float x) { return ::logf(x); } inline double log2(double x) { return ::log2(x); } inline float log2(float x) { return ::log2f(x); } inline double pow(double x, double y) { return ::pow(x, y); } inline float pow(float x, float y) { return ::powf(x, y); } } #endif // MANDEL_TYPES_H