#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r21959 = 1; float r21960 = atan2(1.0, 0.0); float r21961 = 4; float r21962 = r21960 / r21961; float r21963 = r21959 / r21962; float r21964 = f; float r21965 = r21962 * r21964; float r21966 = exp(r21965); float r21967 = -r21965; float r21968 = exp(r21967); float r21969 = r21966 + r21968; float r21970 = r21966 - r21968; float r21971 = r21969 / r21970; float r21972 = log(r21971); float r21973 = r21963 * r21972; float r21974 = -r21973; return r21974; } double f_id(double f) { double r21975 = 1; double r21976 = atan2(1.0, 0.0); double r21977 = 4; double r21978 = r21976 / r21977; double r21979 = r21975 / r21978; double r21980 = f; double r21981 = r21978 * r21980; double r21982 = exp(r21981); double r21983 = -r21981; double r21984 = exp(r21983); double r21985 = r21982 + r21984; double r21986 = r21982 - r21984; double r21987 = r21985 / r21986; double r21988 = log(r21987); double r21989 = r21979 * r21988; double r21990 = -r21989; return r21990; } double f_of(float f) { float r21991 = 4; float r21992 = atan2(1.0, 0.0); float r21993 = r21991 / r21992; float r21994 = log(r21993); float r21995 = r21994 / r21992; float r21996 = r21991 * r21995; float r21997 = 1/12; float r21998 = f; float r21999 = 2; float r22000 = pow(r21998, r21999); float r22001 = r21992 * r22000; float r22002 = r21997 * r22001; float r22003 = r21996 + r22002; float r22004 = 7/5760; float r22005 = 3; float r22006 = pow(r21992, r22005); float r22007 = log1p(r22006); float r22008 = expm1(r22007); float r22009 = pow(r21998, r21991); float r22010 = r22008 * r22009; float r22011 = r22004 * r22010; float r22012 = log(r21998); float r22013 = r22012 / r21992; float r22014 = r21991 * r22013; float r22015 = r22011 + r22014; float r22016 = r22003 - r22015; float r22017 = -r22016; return r22017; } double f_od(double f) { double r22018 = 4; double r22019 = atan2(1.0, 0.0); double r22020 = r22018 / r22019; double r22021 = log(r22020); double r22022 = r22021 / r22019; double r22023 = r22018 * r22022; double r22024 = 1/12; double r22025 = f; double r22026 = 2; double r22027 = pow(r22025, r22026); double r22028 = r22019 * r22027; double r22029 = r22024 * r22028; double r22030 = r22023 + r22029; double r22031 = 7/5760; double r22032 = 3; double r22033 = pow(r22019, r22032); double r22034 = log1p(r22033); double r22035 = expm1(r22034); double r22036 = pow(r22025, r22018); double r22037 = r22035 * r22036; double r22038 = r22031 * r22037; double r22039 = log(r22025); double r22040 = r22039 / r22019; double r22041 = r22018 * r22040; double r22042 = r22038 + r22041; double r22043 = r22030 - r22042; double r22044 = -r22043; return r22044; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r22045, "1", 10, MPFR_RNDN); mpfr_init(r22046); mpfr_init_set_str(r22047, "4", 10, MPFR_RNDN); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); } double f_im(double f) { ; mpfr_const_pi(r22046, MPFR_RNDN); ; mpfr_div(r22048, r22046, r22047, MPFR_RNDN); mpfr_div(r22049, r22045, r22048, MPFR_RNDN); mpfr_set_d(r22050, f, MPFR_RNDN); mpfr_mul(r22051, r22048, r22050, MPFR_RNDN); mpfr_exp(r22052, r22051, MPFR_RNDN); mpfr_neg(r22053, r22051, MPFR_RNDN); mpfr_exp(r22054, r22053, MPFR_RNDN); mpfr_add(r22055, r22052, r22054, MPFR_RNDN); mpfr_sub(r22056, r22052, r22054, MPFR_RNDN); mpfr_div(r22057, r22055, r22056, MPFR_RNDN); mpfr_log(r22058, r22057, MPFR_RNDN); mpfr_mul(r22059, r22049, r22058, MPFR_RNDN); mpfr_neg(r22060, r22059, MPFR_RNDN); return mpfr_get_d(r22060, MPFR_RNDN); } static mpfr_t r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068, r22069, r22070, r22071, r22072, r22073, r22074, r22075, r22076, r22077, r22078, r22079, r22080, r22081, r22082, r22083, r22084, r22085, r22086, r22087; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r22061, "4", 10, MPFR_RNDN); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init_set_str(r22067, "1/12", 10, MPFR_RNDN); mpfr_init(r22068); mpfr_init_set_str(r22069, "2", 10, MPFR_RNDN); mpfr_init(r22070); mpfr_init(r22071); mpfr_init(r22072); mpfr_init(r22073); mpfr_init_set_str(r22074, "7/5760", 10, MPFR_RNDN); mpfr_init_set_str(r22075, "3", 10, MPFR_RNDN); mpfr_init(r22076); mpfr_init(r22077); mpfr_init(r22078); mpfr_init(r22079); mpfr_init(r22080); mpfr_init(r22081); mpfr_init(r22082); mpfr_init(r22083); mpfr_init(r22084); mpfr_init(r22085); mpfr_init(r22086); mpfr_init(r22087); } double f_fm(double f) { ; mpfr_const_pi(r22062, MPFR_RNDN); mpfr_div(r22063, r22061, r22062, MPFR_RNDN); mpfr_log(r22064, r22063, MPFR_RNDN); mpfr_div(r22065, r22064, r22062, MPFR_RNDN); mpfr_mul(r22066, r22061, r22065, MPFR_RNDN); ; mpfr_set_d(r22068, f, MPFR_RNDN); ; mpfr_pow(r22070, r22068, r22069, MPFR_RNDN); mpfr_mul(r22071, r22062, r22070, MPFR_RNDN); mpfr_mul(r22072, r22067, r22071, MPFR_RNDN); mpfr_add(r22073, r22066, r22072, MPFR_RNDN); ; ; mpfr_pow(r22076, r22062, r22075, MPFR_RNDN); mpfr_log1p(r22077, r22076, MPFR_RNDN); mpfr_expm1(r22078, r22077, MPFR_RNDN); mpfr_pow(r22079, r22068, r22061, MPFR_RNDN); mpfr_mul(r22080, r22078, r22079, MPFR_RNDN); mpfr_mul(r22081, r22074, r22080, MPFR_RNDN); mpfr_log(r22082, r22068, MPFR_RNDN); mpfr_div(r22083, r22082, r22062, MPFR_RNDN); mpfr_mul(r22084, r22061, r22083, MPFR_RNDN); mpfr_add(r22085, r22081, r22084, MPFR_RNDN); mpfr_sub(r22086, r22073, r22085, MPFR_RNDN); mpfr_neg(r22087, r22086, MPFR_RNDN); return mpfr_get_d(r22087, MPFR_RNDN); } static mpfr_t r22088, r22089, r22090, r22091, r22092, r22093, r22094, r22095, r22096, r22097, r22098, r22099, r22100, r22101, r22102, r22103, r22104, r22105, r22106, r22107, r22108, r22109, r22110, r22111, r22112, r22113, r22114; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init_set_str(r22088, "4", 10, MPFR_RNDN); mpfr_init(r22089); mpfr_init(r22090); mpfr_init(r22091); mpfr_init(r22092); mpfr_init(r22093); mpfr_init_set_str(r22094, "1/12", 10, MPFR_RNDN); mpfr_init(r22095); mpfr_init_set_str(r22096, "2", 10, MPFR_RNDN); mpfr_init(r22097); mpfr_init(r22098); mpfr_init(r22099); mpfr_init(r22100); mpfr_init_set_str(r22101, "7/5760", 10, MPFR_RNDN); mpfr_init_set_str(r22102, "3", 10, MPFR_RNDN); mpfr_init(r22103); mpfr_init(r22104); mpfr_init(r22105); mpfr_init(r22106); mpfr_init(r22107); mpfr_init(r22108); mpfr_init(r22109); mpfr_init(r22110); mpfr_init(r22111); mpfr_init(r22112); mpfr_init(r22113); mpfr_init(r22114); } double f_dm(double f) { ; mpfr_const_pi(r22089, MPFR_RNDN); mpfr_div(r22090, r22088, r22089, MPFR_RNDN); mpfr_log(r22091, r22090, MPFR_RNDN); mpfr_div(r22092, r22091, r22089, MPFR_RNDN); mpfr_mul(r22093, r22088, r22092, MPFR_RNDN); ; mpfr_set_d(r22095, f, MPFR_RNDN); ; mpfr_pow(r22097, r22095, r22096, MPFR_RNDN); mpfr_mul(r22098, r22089, r22097, MPFR_RNDN); mpfr_mul(r22099, r22094, r22098, MPFR_RNDN); mpfr_add(r22100, r22093, r22099, MPFR_RNDN); ; ; mpfr_pow(r22103, r22089, r22102, MPFR_RNDN); mpfr_log1p(r22104, r22103, MPFR_RNDN); mpfr_expm1(r22105, r22104, MPFR_RNDN); mpfr_pow(r22106, r22095, r22088, MPFR_RNDN); mpfr_mul(r22107, r22105, r22106, MPFR_RNDN); mpfr_mul(r22108, r22101, r22107, MPFR_RNDN); mpfr_log(r22109, r22095, MPFR_RNDN); mpfr_div(r22110, r22109, r22089, MPFR_RNDN); mpfr_mul(r22111, r22088, r22110, MPFR_RNDN); mpfr_add(r22112, r22108, r22111, MPFR_RNDN); mpfr_sub(r22113, r22100, r22112, MPFR_RNDN); mpfr_neg(r22114, r22113, MPFR_RNDN); return mpfr_get_d(r22114, MPFR_RNDN); }