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