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