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