#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r25789 = x; float r25790 = exp(r25789); float r25791 = 1; float r25792 = r25790 - r25791; float r25793 = r25792 / r25789; return r25793; } double f_id(double x) { double r25794 = x; double r25795 = exp(r25794); double r25796 = 1; double r25797 = r25795 - r25796; double r25798 = r25797 / r25794; return r25798; } double f_of(float x) { float r25799 = x; float r25800 = -0.00010281873558823207; bool r25801 = r25799 <= r25800; float r25802 = r25799 + r25799; float r25803 = exp(r25802); float r25804 = 1; float r25805 = r25803 - r25804; float r25806 = exp(r25799); float r25807 = r25806 + r25804; float r25808 = r25805 / r25807; float r25809 = r25808 / r25799; float r25810 = 1/2; float r25811 = 2; float r25812 = pow(r25799, r25811); float r25813 = r25810 * r25812; float r25814 = 1/6; float r25815 = 3; float r25816 = pow(r25799, r25815); float r25817 = r25814 * r25816; float r25818 = r25817 + r25799; float r25819 = r25813 + r25818; float r25820 = r25819 / r25799; float r25821 = r25801 ? r25809 : r25820; return r25821; } double f_od(double x) { double r25822 = x; double r25823 = -0.00010281873558823207; bool r25824 = r25822 <= r25823; double r25825 = r25822 + r25822; double r25826 = exp(r25825); double r25827 = 1; double r25828 = r25826 - r25827; double r25829 = exp(r25822); double r25830 = r25829 + r25827; double r25831 = r25828 / r25830; double r25832 = r25831 / r25822; double r25833 = 1/2; double r25834 = 2; double r25835 = pow(r25822, r25834); double r25836 = r25833 * r25835; double r25837 = 1/6; double r25838 = 3; double r25839 = pow(r25822, r25838); double r25840 = r25837 * r25839; double r25841 = r25840 + r25822; double r25842 = r25836 + r25841; double r25843 = r25842 / r25822; double r25844 = r25824 ? r25832 : r25843; return r25844; } 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 r25845, r25846, r25847, r25848, r25849; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25845); mpfr_init(r25846); mpfr_init_set_str(r25847, "1", 10, MPFR_RNDN); mpfr_init(r25848); mpfr_init(r25849); } double f_im(double x) { mpfr_set_d(r25845, x, MPFR_RNDN); mpfr_exp(r25846, r25845, MPFR_RNDN); ; mpfr_sub(r25848, r25846, r25847, MPFR_RNDN); mpfr_div(r25849, r25848, r25845, MPFR_RNDN); return mpfr_get_d(r25849, MPFR_RNDN); } static mpfr_t r25850, r25851, r25852, r25853, r25854, r25855, r25856, r25857, r25858, r25859, r25860, r25861, r25862, r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871, r25872; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25850); mpfr_init_set_str(r25851, "-0.00010281873558823207", 10, MPFR_RNDN); mpfr_init(r25852); mpfr_init(r25853); mpfr_init(r25854); mpfr_init_set_str(r25855, "1", 10, MPFR_RNDN); mpfr_init(r25856); mpfr_init(r25857); mpfr_init(r25858); mpfr_init(r25859); mpfr_init(r25860); mpfr_init_set_str(r25861, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r25862, "2", 10, MPFR_RNDN); mpfr_init(r25863); mpfr_init(r25864); mpfr_init_set_str(r25865, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r25866, "3", 10, MPFR_RNDN); mpfr_init(r25867); mpfr_init(r25868); mpfr_init(r25869); mpfr_init(r25870); mpfr_init(r25871); mpfr_init(r25872); } double f_fm(double x) { mpfr_set_d(r25850, x, MPFR_RNDN); ; mpfr_set_si(r25852, mpfr_cmp(r25850, r25851) <= 0, MPFR_RNDN); mpfr_add(r25853, r25850, r25850, MPFR_RNDN); mpfr_exp(r25854, r25853, MPFR_RNDN); ; mpfr_sub(r25856, r25854, r25855, MPFR_RNDN); mpfr_exp(r25857, r25850, MPFR_RNDN); mpfr_add(r25858, r25857, r25855, MPFR_RNDN); mpfr_div(r25859, r25856, r25858, MPFR_RNDN); mpfr_div(r25860, r25859, r25850, MPFR_RNDN); ; ; mpfr_pow(r25863, r25850, r25862, MPFR_RNDN); mpfr_mul(r25864, r25861, r25863, MPFR_RNDN); ; ; mpfr_pow(r25867, r25850, r25866, MPFR_RNDN); mpfr_mul(r25868, r25865, r25867, MPFR_RNDN); mpfr_add(r25869, r25868, r25850, MPFR_RNDN); mpfr_add(r25870, r25864, r25869, MPFR_RNDN); mpfr_div(r25871, r25870, r25850, MPFR_RNDN); if (mpfr_get_si(r25852, MPFR_RNDN)) { mpfr_set(r25872, r25860, MPFR_RNDN); } else { mpfr_set(r25872, r25871, MPFR_RNDN); }; return mpfr_get_d(r25872, MPFR_RNDN); } static mpfr_t r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889, r25890, r25891, r25892, r25893, r25894, r25895; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25873); mpfr_init_set_str(r25874, "-0.00010281873558823207", 10, MPFR_RNDN); mpfr_init(r25875); mpfr_init(r25876); mpfr_init(r25877); mpfr_init_set_str(r25878, "1", 10, MPFR_RNDN); mpfr_init(r25879); mpfr_init(r25880); mpfr_init(r25881); mpfr_init(r25882); mpfr_init(r25883); mpfr_init_set_str(r25884, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r25885, "2", 10, MPFR_RNDN); mpfr_init(r25886); mpfr_init(r25887); mpfr_init_set_str(r25888, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r25889, "3", 10, MPFR_RNDN); mpfr_init(r25890); mpfr_init(r25891); mpfr_init(r25892); mpfr_init(r25893); mpfr_init(r25894); mpfr_init(r25895); } double f_dm(double x) { mpfr_set_d(r25873, x, MPFR_RNDN); ; mpfr_set_si(r25875, mpfr_cmp(r25873, r25874) <= 0, MPFR_RNDN); mpfr_add(r25876, r25873, r25873, MPFR_RNDN); mpfr_exp(r25877, r25876, MPFR_RNDN); ; mpfr_sub(r25879, r25877, r25878, MPFR_RNDN); mpfr_exp(r25880, r25873, MPFR_RNDN); mpfr_add(r25881, r25880, r25878, MPFR_RNDN); mpfr_div(r25882, r25879, r25881, MPFR_RNDN); mpfr_div(r25883, r25882, r25873, MPFR_RNDN); ; ; mpfr_pow(r25886, r25873, r25885, MPFR_RNDN); mpfr_mul(r25887, r25884, r25886, MPFR_RNDN); ; ; mpfr_pow(r25890, r25873, r25889, MPFR_RNDN); mpfr_mul(r25891, r25888, r25890, MPFR_RNDN); mpfr_add(r25892, r25891, r25873, MPFR_RNDN); mpfr_add(r25893, r25887, r25892, MPFR_RNDN); mpfr_div(r25894, r25893, r25873, MPFR_RNDN); if (mpfr_get_si(r25875, MPFR_RNDN)) { mpfr_set(r25895, r25883, MPFR_RNDN); } else { mpfr_set(r25895, r25894, MPFR_RNDN); }; return mpfr_get_d(r25895, MPFR_RNDN); }