#include #include #include #include #include char *name = "expq2 (section 3.11)"; double f_if(float x) { float r22420 = x; float r22421 = exp(r22420); float r22422 = 1; float r22423 = r22421 - r22422; float r22424 = r22421 / r22423; return r22424; } double f_id(double x) { double r22425 = x; double r22426 = exp(r22425); double r22427 = 1; double r22428 = r22426 - r22427; double r22429 = r22426 / r22428; return r22429; } double f_of(float x) { float r22430 = x; float r22431 = -r22430; float r22432 = exp(r22431); float r22433 = 1.0015501802622442; bool r22434 = r22432 <= r22433; float r22435 = 1/2; float r22436 = 1; float r22437 = r22436 / r22430; float r22438 = 1/12; float r22439 = r22438 * r22430; float r22440 = r22437 + r22439; float r22441 = r22435 + r22440; float r22442 = exp(r22430); float r22443 = r22430 + r22430; float r22444 = exp(r22443); float r22445 = r22444 - r22436; float r22446 = exp(r22445); float r22447 = log(r22446); float r22448 = r22442 + r22436; float r22449 = r22447 / r22448; float r22450 = r22442 / r22449; float r22451 = r22434 ? r22441 : r22450; return r22451; } double f_od(double x) { double r22452 = x; double r22453 = -r22452; double r22454 = exp(r22453); double r22455 = 1.0015501802622442; bool r22456 = r22454 <= r22455; double r22457 = 1/2; double r22458 = 1; double r22459 = r22458 / r22452; double r22460 = 1/12; double r22461 = r22460 * r22452; double r22462 = r22459 + r22461; double r22463 = r22457 + r22462; double r22464 = exp(r22452); double r22465 = r22452 + r22452; double r22466 = exp(r22465); double r22467 = r22466 - r22458; double r22468 = exp(r22467); double r22469 = log(r22468); double r22470 = r22464 + r22458; double r22471 = r22469 / r22470; double r22472 = r22464 / r22471; double r22473 = r22456 ? r22463 : r22472; return r22473; } 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 r22474, r22475, r22476, r22477, r22478; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22474); mpfr_init(r22475); mpfr_init_set_str(r22476, "1", 10, MPFR_RNDN); mpfr_init(r22477); mpfr_init(r22478); } double f_im(double x) { mpfr_set_d(r22474, x, MPFR_RNDN); mpfr_exp(r22475, r22474, MPFR_RNDN); ; mpfr_sub(r22477, r22475, r22476, MPFR_RNDN); mpfr_div(r22478, r22475, r22477, MPFR_RNDN); return mpfr_get_d(r22478, MPFR_RNDN); } static mpfr_t r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499, r22500; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init_set_str(r22482, "1.0015501802622442", 10, MPFR_RNDN); mpfr_init(r22483); mpfr_init_set_str(r22484, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r22485, "1", 10, MPFR_RNDN); mpfr_init(r22486); mpfr_init_set_str(r22487, "1/12", 10, MPFR_RNDN); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); mpfr_init(r22500); } double f_fm(double x) { mpfr_set_d(r22479, x, MPFR_RNDN); mpfr_neg(r22480, r22479, MPFR_RNDN); mpfr_exp(r22481, r22480, MPFR_RNDN); ; mpfr_set_si(r22483, mpfr_cmp(r22481, r22482) <= 0, MPFR_RNDN); ; ; mpfr_div(r22486, r22485, r22479, MPFR_RNDN); ; mpfr_mul(r22488, r22487, r22479, MPFR_RNDN); mpfr_add(r22489, r22486, r22488, MPFR_RNDN); mpfr_add(r22490, r22484, r22489, MPFR_RNDN); mpfr_exp(r22491, r22479, MPFR_RNDN); mpfr_add(r22492, r22479, r22479, MPFR_RNDN); mpfr_exp(r22493, r22492, MPFR_RNDN); mpfr_sub(r22494, r22493, r22485, MPFR_RNDN); mpfr_exp(r22495, r22494, MPFR_RNDN); mpfr_log(r22496, r22495, MPFR_RNDN); mpfr_add(r22497, r22491, r22485, MPFR_RNDN); mpfr_div(r22498, r22496, r22497, MPFR_RNDN); mpfr_div(r22499, r22491, r22498, MPFR_RNDN); if (mpfr_get_si(r22483, MPFR_RNDN)) { mpfr_set(r22500, r22490, MPFR_RNDN); } else { mpfr_set(r22500, r22499, MPFR_RNDN); }; return mpfr_get_d(r22500, MPFR_RNDN); } static mpfr_t r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519, r22520, r22521, r22522; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init_set_str(r22504, "1.0015501802622442", 10, MPFR_RNDN); mpfr_init(r22505); mpfr_init_set_str(r22506, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r22507, "1", 10, MPFR_RNDN); mpfr_init(r22508); mpfr_init_set_str(r22509, "1/12", 10, MPFR_RNDN); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init(r22513); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); mpfr_init(r22520); mpfr_init(r22521); mpfr_init(r22522); } double f_dm(double x) { mpfr_set_d(r22501, x, MPFR_RNDN); mpfr_neg(r22502, r22501, MPFR_RNDN); mpfr_exp(r22503, r22502, MPFR_RNDN); ; mpfr_set_si(r22505, mpfr_cmp(r22503, r22504) <= 0, MPFR_RNDN); ; ; mpfr_div(r22508, r22507, r22501, MPFR_RNDN); ; mpfr_mul(r22510, r22509, r22501, MPFR_RNDN); mpfr_add(r22511, r22508, r22510, MPFR_RNDN); mpfr_add(r22512, r22506, r22511, MPFR_RNDN); mpfr_exp(r22513, r22501, MPFR_RNDN); mpfr_add(r22514, r22501, r22501, MPFR_RNDN); mpfr_exp(r22515, r22514, MPFR_RNDN); mpfr_sub(r22516, r22515, r22507, MPFR_RNDN); mpfr_exp(r22517, r22516, MPFR_RNDN); mpfr_log(r22518, r22517, MPFR_RNDN); mpfr_add(r22519, r22513, r22507, MPFR_RNDN); mpfr_div(r22520, r22518, r22519, MPFR_RNDN); mpfr_div(r22521, r22513, r22520, MPFR_RNDN); if (mpfr_get_si(r22505, MPFR_RNDN)) { mpfr_set(r22522, r22512, MPFR_RNDN); } else { mpfr_set(r22522, r22521, MPFR_RNDN); }; return mpfr_get_d(r22522, MPFR_RNDN); }