#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r25476 = x; float r25477 = r25476 * r25476; float r25478 = 1; float r25479 = r25477 + r25478; float r25480 = r25476 / r25479; return r25480; } double f_id(double x) { double r25481 = x; double r25482 = r25481 * r25481; double r25483 = 1; double r25484 = r25482 + r25483; double r25485 = r25481 / r25484; return r25485; } double f_of(float x) { float r25486 = x; float r25487 = -7.130874545194846e+21; bool r25488 = r25486 <= r25487; float r25489 = 403.8751728915121; bool r25490 = r25486 <= r25489; float r25491 = !r25490; bool r25492 = r25488 || r25491; float r25493 = 1; float r25494 = 5; float r25495 = pow(r25486, r25494); float r25496 = r25493 / r25495; float r25497 = r25493 / r25486; float r25498 = r25496 + r25497; float r25499 = 3; float r25500 = pow(r25486, r25499); float r25501 = r25493 / r25500; float r25502 = r25498 - r25501; float r25503 = r25486 * r25486; float r25504 = r25503 + r25493; float r25505 = r25486 / r25504; float r25506 = r25492 ? r25502 : r25505; return r25506; } double f_od(double x) { double r25507 = x; double r25508 = -7.130874545194846e+21; bool r25509 = r25507 <= r25508; double r25510 = 403.8751728915121; bool r25511 = r25507 <= r25510; double r25512 = !r25511; bool r25513 = r25509 || r25512; double r25514 = 1; double r25515 = 5; double r25516 = pow(r25507, r25515); double r25517 = r25514 / r25516; double r25518 = r25514 / r25507; double r25519 = r25517 + r25518; double r25520 = 3; double r25521 = pow(r25507, r25520); double r25522 = r25514 / r25521; double r25523 = r25519 - r25522; double r25524 = r25507 * r25507; double r25525 = r25524 + r25514; double r25526 = r25507 / r25525; double r25527 = r25513 ? r25523 : r25526; return r25527; } 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 r25528, r25529, r25530, r25531, r25532; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25528); mpfr_init(r25529); mpfr_init_set_str(r25530, "1", 10, MPFR_RNDN); mpfr_init(r25531); mpfr_init(r25532); } double f_im(double x) { mpfr_set_d(r25528, x, MPFR_RNDN); mpfr_mul(r25529, r25528, r25528, MPFR_RNDN); ; mpfr_add(r25531, r25529, r25530, MPFR_RNDN); mpfr_div(r25532, r25528, r25531, MPFR_RNDN); return mpfr_get_d(r25532, MPFR_RNDN); } static mpfr_t r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25533); mpfr_init_set_str(r25534, "-7.130874545194846e+21", 10, MPFR_RNDN); mpfr_init(r25535); mpfr_init_set_str(r25536, "403.8751728915121", 10, MPFR_RNDN); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init_set_str(r25540, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25541, "5", 10, MPFR_RNDN); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init_set_str(r25546, "3", 10, MPFR_RNDN); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); } double f_fm(double x) { mpfr_set_d(r25533, x, MPFR_RNDN); ; mpfr_set_si(r25535, mpfr_cmp(r25533, r25534) <= 0, MPFR_RNDN); ; mpfr_set_si(r25537, mpfr_cmp(r25533, r25536) <= 0, MPFR_RNDN); mpfr_set_si(r25538, !mpfr_get_si(r25537, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25539, mpfr_get_si(r25535, MPFR_RNDN) || mpfr_get_si(r25538, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25542, r25533, r25541, MPFR_RNDN); mpfr_div(r25543, r25540, r25542, MPFR_RNDN); mpfr_div(r25544, r25540, r25533, MPFR_RNDN); mpfr_add(r25545, r25543, r25544, MPFR_RNDN); ; mpfr_pow(r25547, r25533, r25546, MPFR_RNDN); mpfr_div(r25548, r25540, r25547, MPFR_RNDN); mpfr_sub(r25549, r25545, r25548, MPFR_RNDN); mpfr_mul(r25550, r25533, r25533, MPFR_RNDN); mpfr_add(r25551, r25550, r25540, MPFR_RNDN); mpfr_div(r25552, r25533, r25551, MPFR_RNDN); if (mpfr_get_si(r25539, MPFR_RNDN)) { mpfr_set(r25553, r25549, MPFR_RNDN); } else { mpfr_set(r25553, r25552, MPFR_RNDN); }; return mpfr_get_d(r25553, MPFR_RNDN); } static mpfr_t r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25554); mpfr_init_set_str(r25555, "-7.130874545194846e+21", 10, MPFR_RNDN); mpfr_init(r25556); mpfr_init_set_str(r25557, "403.8751728915121", 10, MPFR_RNDN); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init_set_str(r25561, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25562, "5", 10, MPFR_RNDN); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init_set_str(r25567, "3", 10, MPFR_RNDN); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); mpfr_init(r25574); } double f_dm(double x) { mpfr_set_d(r25554, x, MPFR_RNDN); ; mpfr_set_si(r25556, mpfr_cmp(r25554, r25555) <= 0, MPFR_RNDN); ; mpfr_set_si(r25558, mpfr_cmp(r25554, r25557) <= 0, MPFR_RNDN); mpfr_set_si(r25559, !mpfr_get_si(r25558, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25560, mpfr_get_si(r25556, MPFR_RNDN) || mpfr_get_si(r25559, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25563, r25554, r25562, MPFR_RNDN); mpfr_div(r25564, r25561, r25563, MPFR_RNDN); mpfr_div(r25565, r25561, r25554, MPFR_RNDN); mpfr_add(r25566, r25564, r25565, MPFR_RNDN); ; mpfr_pow(r25568, r25554, r25567, MPFR_RNDN); mpfr_div(r25569, r25561, r25568, MPFR_RNDN); mpfr_sub(r25570, r25566, r25569, MPFR_RNDN); mpfr_mul(r25571, r25554, r25554, MPFR_RNDN); mpfr_add(r25572, r25571, r25561, MPFR_RNDN); mpfr_div(r25573, r25554, r25572, MPFR_RNDN); if (mpfr_get_si(r25560, MPFR_RNDN)) { mpfr_set(r25574, r25570, MPFR_RNDN); } else { mpfr_set(r25574, r25573, MPFR_RNDN); }; return mpfr_get_d(r25574, MPFR_RNDN); }