#include #include #include #include #include char *name = "Kahan's exp quotient"; double f_if(float x) { float r25986 = x; float r25987 = exp(r25986); float r25988 = 1; float r25989 = r25987 - r25988; float r25990 = r25989 / r25986; return r25990; } double f_id(double x) { double r25991 = x; double r25992 = exp(r25991); double r25993 = 1; double r25994 = r25992 - r25993; double r25995 = r25994 / r25991; return r25995; } double f_of(float x) { float r25996 = x; float r25997 = -0.00015598560706120413; bool r25998 = r25996 <= r25997; float r25999 = exp(r25996); float r26000 = 3; float r26001 = pow(r25999, r26000); float r26002 = 1; float r26003 = r26001 - r26002; float r26004 = r25996 + r25996; float r26005 = exp(r26004); float r26006 = r25999 + r26002; float r26007 = r26005 + r26006; float r26008 = r26003 / r26007; float r26009 = r26008 / r25996; float r26010 = 1/6; float r26011 = 2; float r26012 = pow(r25996, r26011); float r26013 = r26010 * r26012; float r26014 = 1/2; float r26015 = r26014 * r25996; float r26016 = r26002 + r26015; float r26017 = r26013 + r26016; float r26018 = r25998 ? r26009 : r26017; return r26018; } double f_od(double x) { double r26019 = x; double r26020 = -0.00015598560706120413; bool r26021 = r26019 <= r26020; double r26022 = exp(r26019); double r26023 = 3; double r26024 = pow(r26022, r26023); double r26025 = 1; double r26026 = r26024 - r26025; double r26027 = r26019 + r26019; double r26028 = exp(r26027); double r26029 = r26022 + r26025; double r26030 = r26028 + r26029; double r26031 = r26026 / r26030; double r26032 = r26031 / r26019; double r26033 = 1/6; double r26034 = 2; double r26035 = pow(r26019, r26034); double r26036 = r26033 * r26035; double r26037 = 1/2; double r26038 = r26037 * r26019; double r26039 = r26025 + r26038; double r26040 = r26036 + r26039; double r26041 = r26021 ? r26032 : r26040; return r26041; } 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 r26042, r26043, r26044, r26045, r26046; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26042); mpfr_init(r26043); mpfr_init_set_str(r26044, "1", 10, MPFR_RNDN); mpfr_init(r26045); mpfr_init(r26046); } double f_im(double x) { mpfr_set_d(r26042, x, MPFR_RNDN); mpfr_exp(r26043, r26042, MPFR_RNDN); ; mpfr_sub(r26045, r26043, r26044, MPFR_RNDN); mpfr_div(r26046, r26045, r26042, MPFR_RNDN); return mpfr_get_d(r26046, MPFR_RNDN); } static mpfr_t r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26047); mpfr_init_set_str(r26048, "-0.00015598560706120413", 10, MPFR_RNDN); mpfr_init(r26049); mpfr_init(r26050); mpfr_init_set_str(r26051, "3", 10, MPFR_RNDN); mpfr_init(r26052); mpfr_init_set_str(r26053, "1", 10, MPFR_RNDN); mpfr_init(r26054); mpfr_init(r26055); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init_set_str(r26061, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r26062, "2", 10, MPFR_RNDN); mpfr_init(r26063); mpfr_init(r26064); mpfr_init_set_str(r26065, "1/2", 10, MPFR_RNDN); mpfr_init(r26066); mpfr_init(r26067); mpfr_init(r26068); mpfr_init(r26069); } double f_fm(double x) { mpfr_set_d(r26047, x, MPFR_RNDN); ; mpfr_set_si(r26049, mpfr_cmp(r26047, r26048) <= 0, MPFR_RNDN); mpfr_exp(r26050, r26047, MPFR_RNDN); ; mpfr_pow(r26052, r26050, r26051, MPFR_RNDN); ; mpfr_sub(r26054, r26052, r26053, MPFR_RNDN); mpfr_add(r26055, r26047, r26047, MPFR_RNDN); mpfr_exp(r26056, r26055, MPFR_RNDN); mpfr_add(r26057, r26050, r26053, MPFR_RNDN); mpfr_add(r26058, r26056, r26057, MPFR_RNDN); mpfr_div(r26059, r26054, r26058, MPFR_RNDN); mpfr_div(r26060, r26059, r26047, MPFR_RNDN); ; ; mpfr_pow(r26063, r26047, r26062, MPFR_RNDN); mpfr_mul(r26064, r26061, r26063, MPFR_RNDN); ; mpfr_mul(r26066, r26065, r26047, MPFR_RNDN); mpfr_add(r26067, r26053, r26066, MPFR_RNDN); mpfr_add(r26068, r26064, r26067, MPFR_RNDN); if (mpfr_get_si(r26049, MPFR_RNDN)) { mpfr_set(r26069, r26060, MPFR_RNDN); } else { mpfr_set(r26069, r26068, MPFR_RNDN); }; return mpfr_get_d(r26069, MPFR_RNDN); } static mpfr_t r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26070); mpfr_init_set_str(r26071, "-0.00015598560706120413", 10, MPFR_RNDN); mpfr_init(r26072); mpfr_init(r26073); mpfr_init_set_str(r26074, "3", 10, MPFR_RNDN); mpfr_init(r26075); mpfr_init_set_str(r26076, "1", 10, MPFR_RNDN); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init_set_str(r26084, "1/6", 10, MPFR_RNDN); mpfr_init_set_str(r26085, "2", 10, MPFR_RNDN); mpfr_init(r26086); mpfr_init(r26087); mpfr_init_set_str(r26088, "1/2", 10, MPFR_RNDN); mpfr_init(r26089); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); } double f_dm(double x) { mpfr_set_d(r26070, x, MPFR_RNDN); ; mpfr_set_si(r26072, mpfr_cmp(r26070, r26071) <= 0, MPFR_RNDN); mpfr_exp(r26073, r26070, MPFR_RNDN); ; mpfr_pow(r26075, r26073, r26074, MPFR_RNDN); ; mpfr_sub(r26077, r26075, r26076, MPFR_RNDN); mpfr_add(r26078, r26070, r26070, MPFR_RNDN); mpfr_exp(r26079, r26078, MPFR_RNDN); mpfr_add(r26080, r26073, r26076, MPFR_RNDN); mpfr_add(r26081, r26079, r26080, MPFR_RNDN); mpfr_div(r26082, r26077, r26081, MPFR_RNDN); mpfr_div(r26083, r26082, r26070, MPFR_RNDN); ; ; mpfr_pow(r26086, r26070, r26085, MPFR_RNDN); mpfr_mul(r26087, r26084, r26086, MPFR_RNDN); ; mpfr_mul(r26089, r26088, r26070, MPFR_RNDN); mpfr_add(r26090, r26076, r26089, MPFR_RNDN); mpfr_add(r26091, r26087, r26090, MPFR_RNDN); if (mpfr_get_si(r26072, MPFR_RNDN)) { mpfr_set(r26092, r26083, MPFR_RNDN); } else { mpfr_set(r26092, r26091, MPFR_RNDN); }; return mpfr_get_d(r26092, MPFR_RNDN); }