#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22507 = x; float r22508 = r22507 * r22507; float r22509 = 1; float r22510 = r22508 + r22509; float r22511 = r22507 / r22510; return r22511; } double f_id(double x) { double r22512 = x; double r22513 = r22512 * r22512; double r22514 = 1; double r22515 = r22513 + r22514; double r22516 = r22512 / r22515; return r22516; } double f_of(float x) { float r22517 = x; float r22518 = -19409101.03762743; bool r22519 = r22517 <= r22518; float r22520 = 1; float r22521 = 5; float r22522 = pow(r22517, r22521); float r22523 = r22520 / r22522; float r22524 = r22520 / r22517; float r22525 = r22523 + r22524; float r22526 = 3; float r22527 = pow(r22517, r22526); float r22528 = r22520 / r22527; float r22529 = r22525 - r22528; float r22530 = 144528077.66840905; bool r22531 = r22517 <= r22530; float r22532 = r22517 * r22517; float r22533 = r22532 + r22520; float r22534 = r22517 / r22533; float r22535 = r22531 ? r22534 : r22529; float r22536 = r22519 ? r22529 : r22535; return r22536; } double f_od(double x) { double r22537 = x; double r22538 = -19409101.03762743; bool r22539 = r22537 <= r22538; double r22540 = 1; double r22541 = 5; double r22542 = pow(r22537, r22541); double r22543 = r22540 / r22542; double r22544 = r22540 / r22537; double r22545 = r22543 + r22544; double r22546 = 3; double r22547 = pow(r22537, r22546); double r22548 = r22540 / r22547; double r22549 = r22545 - r22548; double r22550 = 144528077.66840905; bool r22551 = r22537 <= r22550; double r22552 = r22537 * r22537; double r22553 = r22552 + r22540; double r22554 = r22537 / r22553; double r22555 = r22551 ? r22554 : r22549; double r22556 = r22539 ? r22549 : r22555; return r22556; } 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 r22557, r22558, r22559, r22560, r22561; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22557); mpfr_init(r22558); mpfr_init_set_str(r22559, "1", 10, MPFR_RNDN); mpfr_init(r22560); mpfr_init(r22561); } double f_im(double x) { mpfr_set_d(r22557, x, MPFR_RNDN); mpfr_mul(r22558, r22557, r22557, MPFR_RNDN); ; mpfr_add(r22560, r22558, r22559, MPFR_RNDN); mpfr_div(r22561, r22557, r22560, MPFR_RNDN); return mpfr_get_d(r22561, MPFR_RNDN); } static mpfr_t r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22562); mpfr_init_set_str(r22563, "-19409101.03762743", 10, MPFR_RNDN); mpfr_init(r22564); mpfr_init_set_str(r22565, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22566, "5", 10, MPFR_RNDN); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init(r22570); mpfr_init_set_str(r22571, "3", 10, MPFR_RNDN); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init_set_str(r22575, "144528077.66840905", 10, MPFR_RNDN); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); } double f_fm(double x) { mpfr_set_d(r22562, x, MPFR_RNDN); ; mpfr_set_si(r22564, mpfr_cmp(r22562, r22563) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22567, r22562, r22566, MPFR_RNDN); mpfr_div(r22568, r22565, r22567, MPFR_RNDN); mpfr_div(r22569, r22565, r22562, MPFR_RNDN); mpfr_add(r22570, r22568, r22569, MPFR_RNDN); ; mpfr_pow(r22572, r22562, r22571, MPFR_RNDN); mpfr_div(r22573, r22565, r22572, MPFR_RNDN); mpfr_sub(r22574, r22570, r22573, MPFR_RNDN); ; mpfr_set_si(r22576, mpfr_cmp(r22562, r22575) <= 0, MPFR_RNDN); mpfr_mul(r22577, r22562, r22562, MPFR_RNDN); mpfr_add(r22578, r22577, r22565, MPFR_RNDN); mpfr_div(r22579, r22562, r22578, MPFR_RNDN); if (mpfr_get_si(r22576, MPFR_RNDN)) { mpfr_set(r22580, r22579, MPFR_RNDN); } else { mpfr_set(r22580, r22574, MPFR_RNDN); }; if (mpfr_get_si(r22564, MPFR_RNDN)) { mpfr_set(r22581, r22574, MPFR_RNDN); } else { mpfr_set(r22581, r22580, MPFR_RNDN); }; return mpfr_get_d(r22581, MPFR_RNDN); } static mpfr_t r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22582); mpfr_init_set_str(r22583, "-19409101.03762743", 10, MPFR_RNDN); mpfr_init(r22584); mpfr_init_set_str(r22585, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22586, "5", 10, MPFR_RNDN); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init_set_str(r22591, "3", 10, MPFR_RNDN); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init_set_str(r22595, "144528077.66840905", 10, MPFR_RNDN); mpfr_init(r22596); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); } double f_dm(double x) { mpfr_set_d(r22582, x, MPFR_RNDN); ; mpfr_set_si(r22584, mpfr_cmp(r22582, r22583) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22587, r22582, r22586, MPFR_RNDN); mpfr_div(r22588, r22585, r22587, MPFR_RNDN); mpfr_div(r22589, r22585, r22582, MPFR_RNDN); mpfr_add(r22590, r22588, r22589, MPFR_RNDN); ; mpfr_pow(r22592, r22582, r22591, MPFR_RNDN); mpfr_div(r22593, r22585, r22592, MPFR_RNDN); mpfr_sub(r22594, r22590, r22593, MPFR_RNDN); ; mpfr_set_si(r22596, mpfr_cmp(r22582, r22595) <= 0, MPFR_RNDN); mpfr_mul(r22597, r22582, r22582, MPFR_RNDN); mpfr_add(r22598, r22597, r22585, MPFR_RNDN); mpfr_div(r22599, r22582, r22598, MPFR_RNDN); if (mpfr_get_si(r22596, MPFR_RNDN)) { mpfr_set(r22600, r22599, MPFR_RNDN); } else { mpfr_set(r22600, r22594, MPFR_RNDN); }; if (mpfr_get_si(r22584, MPFR_RNDN)) { mpfr_set(r22601, r22594, MPFR_RNDN); } else { mpfr_set(r22601, r22600, MPFR_RNDN); }; return mpfr_get_d(r22601, MPFR_RNDN); }