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