#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r20764 = 1; float r20765 = atan2(1.0, 0.0); float r20766 = 4; float r20767 = r20765 / r20766; float r20768 = r20764 / r20767; float r20769 = f; float r20770 = r20767 * r20769; float r20771 = exp(r20770); float r20772 = -r20770; float r20773 = exp(r20772); float r20774 = r20771 + r20773; float r20775 = r20771 - r20773; float r20776 = r20774 / r20775; float r20777 = log(r20776); float r20778 = r20768 * r20777; float r20779 = -r20778; return r20779; } double f_id(double f) { double r20780 = 1; double r20781 = atan2(1.0, 0.0); double r20782 = 4; double r20783 = r20781 / r20782; double r20784 = r20780 / r20783; double r20785 = f; double r20786 = r20783 * r20785; double r20787 = exp(r20786); double r20788 = -r20786; double r20789 = exp(r20788); double r20790 = r20787 + r20789; double r20791 = r20787 - r20789; double r20792 = r20790 / r20791; double r20793 = log(r20792); double r20794 = r20784 * r20793; double r20795 = -r20794; return r20795; } double f_of(float f) { float r20796 = 4; float r20797 = atan2(1.0, 0.0); float r20798 = r20796 / r20797; float r20799 = cbrt(r20798); float r20800 = log(r20799); float r20801 = r20796 + r20796; float r20802 = r20801 / r20797; float r20803 = r20800 * r20802; float r20804 = f; float r20805 = log(r20804); float r20806 = r20800 - r20805; float r20807 = r20798 * r20806; float r20808 = r20803 + r20807; float r20809 = r20804 * r20804; float r20810 = 1/12; float r20811 = r20810 * r20797; float r20812 = r20809 * r20811; float r20813 = 7/5760; float r20814 = r20813 * r20797; float r20815 = pow(r20804, r20796); float r20816 = r20814 * r20815; float r20817 = r20797 * r20797; float r20818 = r20816 * r20817; float r20819 = r20812 - r20818; float r20820 = r20808 + r20819; float r20821 = -r20820; return r20821; } double f_od(double f) { double r20822 = 4; double r20823 = atan2(1.0, 0.0); double r20824 = r20822 / r20823; double r20825 = cbrt(r20824); double r20826 = log(r20825); double r20827 = r20822 + r20822; double r20828 = r20827 / r20823; double r20829 = r20826 * r20828; double r20830 = f; double r20831 = log(r20830); double r20832 = r20826 - r20831; double r20833 = r20824 * r20832; double r20834 = r20829 + r20833; double r20835 = r20830 * r20830; double r20836 = 1/12; double r20837 = r20836 * r20823; double r20838 = r20835 * r20837; double r20839 = 7/5760; double r20840 = r20839 * r20823; double r20841 = pow(r20830, r20822); double r20842 = r20840 * r20841; double r20843 = r20823 * r20823; double r20844 = r20842 * r20843; double r20845 = r20838 - r20844; double r20846 = r20834 + r20845; double r20847 = -r20846; return r20847; } 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 r20848, r20849, r20850, r20851, r20852, r20853, r20854, r20855, r20856, r20857, r20858, r20859, r20860, r20861, r20862, r20863; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20848, "1", 10, MPFR_RNDN); mpfr_init(r20849); mpfr_init_set_str(r20850, "4", 10, MPFR_RNDN); mpfr_init(r20851); mpfr_init(r20852); mpfr_init(r20853); mpfr_init(r20854); mpfr_init(r20855); mpfr_init(r20856); mpfr_init(r20857); mpfr_init(r20858); mpfr_init(r20859); mpfr_init(r20860); mpfr_init(r20861); mpfr_init(r20862); mpfr_init(r20863); } double f_im(double f) { ; mpfr_const_pi(r20849, MPFR_RNDN); ; mpfr_div(r20851, r20849, r20850, MPFR_RNDN); mpfr_div(r20852, r20848, r20851, MPFR_RNDN); mpfr_set_d(r20853, f, MPFR_RNDN); mpfr_mul(r20854, r20851, r20853, MPFR_RNDN); mpfr_exp(r20855, r20854, MPFR_RNDN); mpfr_neg(r20856, r20854, MPFR_RNDN); mpfr_exp(r20857, r20856, MPFR_RNDN); mpfr_add(r20858, r20855, r20857, MPFR_RNDN); mpfr_sub(r20859, r20855, r20857, MPFR_RNDN); mpfr_div(r20860, r20858, r20859, MPFR_RNDN); mpfr_log(r20861, r20860, MPFR_RNDN); mpfr_mul(r20862, r20852, r20861, MPFR_RNDN); mpfr_neg(r20863, r20862, MPFR_RNDN); return mpfr_get_d(r20863, MPFR_RNDN); } static mpfr_t r20864, r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872, r20873, r20874, r20875, r20876, r20877, r20878, r20879, r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20864, "4", 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); mpfr_init(r20873); mpfr_init(r20874); mpfr_init(r20875); mpfr_init(r20876); mpfr_init(r20877); mpfr_init_set_str(r20878, "1/12", 10, MPFR_RNDN); mpfr_init(r20879); mpfr_init(r20880); mpfr_init_set_str(r20881, "7/5760", 10, MPFR_RNDN); mpfr_init(r20882); mpfr_init(r20883); mpfr_init(r20884); mpfr_init(r20885); mpfr_init(r20886); mpfr_init(r20887); mpfr_init(r20888); mpfr_init(r20889); } double f_fm(double f) { ; mpfr_const_pi(r20865, MPFR_RNDN); mpfr_div(r20866, r20864, r20865, MPFR_RNDN); mpfr_cbrt(r20867, r20866, MPFR_RNDN); mpfr_log(r20868, r20867, MPFR_RNDN); mpfr_add(r20869, r20864, r20864, MPFR_RNDN); mpfr_div(r20870, r20869, r20865, MPFR_RNDN); mpfr_mul(r20871, r20868, r20870, MPFR_RNDN); mpfr_set_d(r20872, f, MPFR_RNDN); mpfr_log(r20873, r20872, MPFR_RNDN); mpfr_sub(r20874, r20868, r20873, MPFR_RNDN); mpfr_mul(r20875, r20866, r20874, MPFR_RNDN); mpfr_add(r20876, r20871, r20875, MPFR_RNDN); mpfr_mul(r20877, r20872, r20872, MPFR_RNDN); ; mpfr_mul(r20879, r20878, r20865, MPFR_RNDN); mpfr_mul(r20880, r20877, r20879, MPFR_RNDN); ; mpfr_mul(r20882, r20881, r20865, MPFR_RNDN); mpfr_pow(r20883, r20872, r20864, MPFR_RNDN); mpfr_mul(r20884, r20882, r20883, MPFR_RNDN); mpfr_mul(r20885, r20865, r20865, MPFR_RNDN); mpfr_mul(r20886, r20884, r20885, MPFR_RNDN); mpfr_sub(r20887, r20880, r20886, MPFR_RNDN); mpfr_add(r20888, r20876, r20887, MPFR_RNDN); mpfr_neg(r20889, r20888, MPFR_RNDN); return mpfr_get_d(r20889, MPFR_RNDN); } static mpfr_t r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20890, "4", 10, MPFR_RNDN); mpfr_init(r20891); mpfr_init(r20892); mpfr_init(r20893); mpfr_init(r20894); mpfr_init(r20895); mpfr_init(r20896); mpfr_init(r20897); mpfr_init(r20898); mpfr_init(r20899); mpfr_init(r20900); mpfr_init(r20901); mpfr_init(r20902); mpfr_init(r20903); mpfr_init_set_str(r20904, "1/12", 10, MPFR_RNDN); mpfr_init(r20905); mpfr_init(r20906); mpfr_init_set_str(r20907, "7/5760", 10, MPFR_RNDN); mpfr_init(r20908); mpfr_init(r20909); mpfr_init(r20910); mpfr_init(r20911); mpfr_init(r20912); mpfr_init(r20913); mpfr_init(r20914); mpfr_init(r20915); } double f_dm(double f) { ; mpfr_const_pi(r20891, MPFR_RNDN); mpfr_div(r20892, r20890, r20891, MPFR_RNDN); mpfr_cbrt(r20893, r20892, MPFR_RNDN); mpfr_log(r20894, r20893, MPFR_RNDN); mpfr_add(r20895, r20890, r20890, MPFR_RNDN); mpfr_div(r20896, r20895, r20891, MPFR_RNDN); mpfr_mul(r20897, r20894, r20896, MPFR_RNDN); mpfr_set_d(r20898, f, MPFR_RNDN); mpfr_log(r20899, r20898, MPFR_RNDN); mpfr_sub(r20900, r20894, r20899, MPFR_RNDN); mpfr_mul(r20901, r20892, r20900, MPFR_RNDN); mpfr_add(r20902, r20897, r20901, MPFR_RNDN); mpfr_mul(r20903, r20898, r20898, MPFR_RNDN); ; mpfr_mul(r20905, r20904, r20891, MPFR_RNDN); mpfr_mul(r20906, r20903, r20905, MPFR_RNDN); ; mpfr_mul(r20908, r20907, r20891, MPFR_RNDN); mpfr_pow(r20909, r20898, r20890, MPFR_RNDN); mpfr_mul(r20910, r20908, r20909, MPFR_RNDN); mpfr_mul(r20911, r20891, r20891, MPFR_RNDN); mpfr_mul(r20912, r20910, r20911, MPFR_RNDN); mpfr_sub(r20913, r20906, r20912, MPFR_RNDN); mpfr_add(r20914, r20902, r20913, MPFR_RNDN); mpfr_neg(r20915, r20914, MPFR_RNDN); return mpfr_get_d(r20915, MPFR_RNDN); }