#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r20745 = 1; float r20746 = atan2(1.0, 0.0); float r20747 = 4; float r20748 = r20746 / r20747; float r20749 = r20745 / r20748; float r20750 = f; float r20751 = r20748 * r20750; float r20752 = exp(r20751); float r20753 = -r20751; float r20754 = exp(r20753); float r20755 = r20752 + r20754; float r20756 = r20752 - r20754; float r20757 = r20755 / r20756; float r20758 = log(r20757); float r20759 = r20749 * r20758; float r20760 = -r20759; return r20760; } double f_id(double f) { double r20761 = 1; double r20762 = atan2(1.0, 0.0); double r20763 = 4; double r20764 = r20762 / r20763; double r20765 = r20761 / r20764; double r20766 = f; double r20767 = r20764 * r20766; double r20768 = exp(r20767); double r20769 = -r20767; double r20770 = exp(r20769); double r20771 = r20768 + r20770; double r20772 = r20768 - r20770; double r20773 = r20771 / r20772; double r20774 = log(r20773); double r20775 = r20765 * r20774; double r20776 = -r20775; return r20776; } double f_of(float f) { float r20777 = f; float r20778 = 1/12; float r20779 = r20777 * r20778; float r20780 = atan2(1.0, 0.0); float r20781 = r20779 * r20780; float r20782 = 4; float r20783 = r20782 / r20780; float r20784 = r20783 / r20777; float r20785 = r20781 + r20784; float r20786 = 1/2880; float r20787 = r20786 * r20780; float r20788 = 3; float r20789 = pow(r20777, r20788); float r20790 = r20787 * r20789; float r20791 = r20780 * r20780; float r20792 = r20790 * r20791; float r20793 = r20785 - r20792; float r20794 = log(r20793); float r20795 = -r20783; float r20796 = r20794 * r20795; return r20796; } double f_od(double f) { double r20797 = f; double r20798 = 1/12; double r20799 = r20797 * r20798; double r20800 = atan2(1.0, 0.0); double r20801 = r20799 * r20800; double r20802 = 4; double r20803 = r20802 / r20800; double r20804 = r20803 / r20797; double r20805 = r20801 + r20804; double r20806 = 1/2880; double r20807 = r20806 * r20800; double r20808 = 3; double r20809 = pow(r20797, r20808); double r20810 = r20807 * r20809; double r20811 = r20800 * r20800; double r20812 = r20810 * r20811; double r20813 = r20805 - r20812; double r20814 = log(r20813); double r20815 = -r20803; double r20816 = r20814 * r20815; return r20816; } 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 r20817, r20818, r20819, r20820, r20821, r20822, r20823, r20824, r20825, r20826, r20827, r20828, r20829, r20830, r20831, r20832; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20817, "1", 10, MPFR_RNDN); mpfr_init(r20818); mpfr_init_set_str(r20819, "4", 10, MPFR_RNDN); mpfr_init(r20820); mpfr_init(r20821); mpfr_init(r20822); mpfr_init(r20823); mpfr_init(r20824); mpfr_init(r20825); mpfr_init(r20826); mpfr_init(r20827); mpfr_init(r20828); mpfr_init(r20829); mpfr_init(r20830); mpfr_init(r20831); mpfr_init(r20832); } double f_im(double f) { ; mpfr_const_pi(r20818, MPFR_RNDN); ; mpfr_div(r20820, r20818, r20819, MPFR_RNDN); mpfr_div(r20821, r20817, r20820, MPFR_RNDN); mpfr_set_d(r20822, f, MPFR_RNDN); mpfr_mul(r20823, r20820, r20822, MPFR_RNDN); mpfr_exp(r20824, r20823, MPFR_RNDN); mpfr_neg(r20825, r20823, MPFR_RNDN); mpfr_exp(r20826, r20825, MPFR_RNDN); mpfr_add(r20827, r20824, r20826, MPFR_RNDN); mpfr_sub(r20828, r20824, r20826, MPFR_RNDN); mpfr_div(r20829, r20827, r20828, MPFR_RNDN); mpfr_log(r20830, r20829, MPFR_RNDN); mpfr_mul(r20831, r20821, r20830, MPFR_RNDN); mpfr_neg(r20832, r20831, MPFR_RNDN); return mpfr_get_d(r20832, MPFR_RNDN); } static mpfr_t r20833, r20834, r20835, r20836, r20837, r20838, r20839, r20840, r20841, r20842, r20843, r20844, r20845, r20846, r20847, r20848, r20849, r20850, r20851, r20852; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20833); mpfr_init_set_str(r20834, "1/12", 10, MPFR_RNDN); mpfr_init(r20835); mpfr_init(r20836); mpfr_init(r20837); mpfr_init_set_str(r20838, "4", 10, MPFR_RNDN); mpfr_init(r20839); mpfr_init(r20840); mpfr_init(r20841); mpfr_init_set_str(r20842, "1/2880", 10, MPFR_RNDN); mpfr_init(r20843); mpfr_init_set_str(r20844, "3", 10, MPFR_RNDN); mpfr_init(r20845); mpfr_init(r20846); mpfr_init(r20847); mpfr_init(r20848); mpfr_init(r20849); mpfr_init(r20850); mpfr_init(r20851); mpfr_init(r20852); } double f_fm(double f) { mpfr_set_d(r20833, f, MPFR_RNDN); ; mpfr_mul(r20835, r20833, r20834, MPFR_RNDN); mpfr_const_pi(r20836, MPFR_RNDN); mpfr_mul(r20837, r20835, r20836, MPFR_RNDN); ; mpfr_div(r20839, r20838, r20836, MPFR_RNDN); mpfr_div(r20840, r20839, r20833, MPFR_RNDN); mpfr_add(r20841, r20837, r20840, MPFR_RNDN); ; mpfr_mul(r20843, r20842, r20836, MPFR_RNDN); ; mpfr_pow(r20845, r20833, r20844, MPFR_RNDN); mpfr_mul(r20846, r20843, r20845, MPFR_RNDN); mpfr_mul(r20847, r20836, r20836, MPFR_RNDN); mpfr_mul(r20848, r20846, r20847, MPFR_RNDN); mpfr_sub(r20849, r20841, r20848, MPFR_RNDN); mpfr_log(r20850, r20849, MPFR_RNDN); mpfr_neg(r20851, r20839, MPFR_RNDN); mpfr_mul(r20852, r20850, r20851, MPFR_RNDN); return mpfr_get_d(r20852, MPFR_RNDN); } static mpfr_t r20853, r20854, r20855, r20856, r20857, r20858, r20859, r20860, r20861, r20862, r20863, r20864, r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20853); mpfr_init_set_str(r20854, "1/12", 10, MPFR_RNDN); mpfr_init(r20855); mpfr_init(r20856); mpfr_init(r20857); mpfr_init_set_str(r20858, "4", 10, MPFR_RNDN); mpfr_init(r20859); mpfr_init(r20860); mpfr_init(r20861); mpfr_init_set_str(r20862, "1/2880", 10, MPFR_RNDN); mpfr_init(r20863); mpfr_init_set_str(r20864, "3", 10, MPFR_RNDN); mpfr_init(r20865); mpfr_init(r20866); mpfr_init(r20867); mpfr_init(r20868); mpfr_init(r20869); mpfr_init(r20870); mpfr_init(r20871); mpfr_init(r20872); } double f_dm(double f) { mpfr_set_d(r20853, f, MPFR_RNDN); ; mpfr_mul(r20855, r20853, r20854, MPFR_RNDN); mpfr_const_pi(r20856, MPFR_RNDN); mpfr_mul(r20857, r20855, r20856, MPFR_RNDN); ; mpfr_div(r20859, r20858, r20856, MPFR_RNDN); mpfr_div(r20860, r20859, r20853, MPFR_RNDN); mpfr_add(r20861, r20857, r20860, MPFR_RNDN); ; mpfr_mul(r20863, r20862, r20856, MPFR_RNDN); ; mpfr_pow(r20865, r20853, r20864, MPFR_RNDN); mpfr_mul(r20866, r20863, r20865, MPFR_RNDN); mpfr_mul(r20867, r20856, r20856, MPFR_RNDN); mpfr_mul(r20868, r20866, r20867, MPFR_RNDN); mpfr_sub(r20869, r20861, r20868, MPFR_RNDN); mpfr_log(r20870, r20869, MPFR_RNDN); mpfr_neg(r20871, r20859, MPFR_RNDN); mpfr_mul(r20872, r20870, r20871, MPFR_RNDN); return mpfr_get_d(r20872, MPFR_RNDN); }