#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22547 = x; float r22548 = r22547 * r22547; float r22549 = 1; float r22550 = r22548 + r22549; float r22551 = r22547 / r22550; return r22551; } double f_id(double x) { double r22552 = x; double r22553 = r22552 * r22552; double r22554 = 1; double r22555 = r22553 + r22554; double r22556 = r22552 / r22555; return r22556; } double f_of(float x) { float r22557 = x; float r22558 = -1335004178484.0273; bool r22559 = r22557 <= r22558; float r22560 = 891429.0791613997; bool r22561 = r22557 <= r22560; float r22562 = !r22561; bool r22563 = r22559 || r22562; float r22564 = 1; float r22565 = 5; float r22566 = pow(r22557, r22565); float r22567 = r22564 / r22566; float r22568 = r22564 / r22557; float r22569 = r22567 + r22568; float r22570 = 3; float r22571 = pow(r22557, r22570); float r22572 = r22564 / r22571; float r22573 = r22569 - r22572; float r22574 = r22557 * r22557; float r22575 = r22574 + r22564; float r22576 = r22557 / r22575; float r22577 = r22563 ? r22573 : r22576; return r22577; } double f_od(double x) { double r22578 = x; double r22579 = -1335004178484.0273; bool r22580 = r22578 <= r22579; double r22581 = 891429.0791613997; bool r22582 = r22578 <= r22581; double r22583 = !r22582; bool r22584 = r22580 || r22583; double r22585 = 1; double r22586 = 5; double r22587 = pow(r22578, r22586); double r22588 = r22585 / r22587; double r22589 = r22585 / r22578; double r22590 = r22588 + r22589; double r22591 = 3; double r22592 = pow(r22578, r22591); double r22593 = r22585 / r22592; double r22594 = r22590 - r22593; double r22595 = r22578 * r22578; double r22596 = r22595 + r22585; double r22597 = r22578 / r22596; double r22598 = r22584 ? r22594 : r22597; return r22598; } 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 r22599, r22600, r22601, r22602, r22603; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22599); mpfr_init(r22600); mpfr_init_set_str(r22601, "1", 10, MPFR_RNDN); mpfr_init(r22602); mpfr_init(r22603); } double f_im(double x) { mpfr_set_d(r22599, x, MPFR_RNDN); mpfr_mul(r22600, r22599, r22599, MPFR_RNDN); ; mpfr_add(r22602, r22600, r22601, MPFR_RNDN); mpfr_div(r22603, r22599, r22602, MPFR_RNDN); return mpfr_get_d(r22603, MPFR_RNDN); } static mpfr_t r22604, r22605, r22606, r22607, r22608, r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619, r22620, r22621, r22622, r22623, r22624; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r22604); mpfr_init_set_str(r22605, "-1335004178484.0273", 10, MPFR_RNDN); mpfr_init(r22606); mpfr_init_set_str(r22607, "891429.0791613997", 10, MPFR_RNDN); mpfr_init(r22608); mpfr_init(r22609); mpfr_init(r22610); mpfr_init_set_str(r22611, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22612, "5", 10, MPFR_RNDN); mpfr_init(r22613); mpfr_init(r22614); mpfr_init(r22615); mpfr_init(r22616); mpfr_init_set_str(r22617, "3", 10, MPFR_RNDN); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); mpfr_init(r22621); mpfr_init(r22622); mpfr_init(r22623); mpfr_init(r22624); } double f_fm(double x) { mpfr_set_d(r22604, x, MPFR_RNDN); ; mpfr_set_si(r22606, mpfr_cmp(r22604, r22605) <= 0, MPFR_RNDN); ; mpfr_set_si(r22608, mpfr_cmp(r22604, r22607) <= 0, MPFR_RNDN); mpfr_set_si(r22609, !mpfr_get_si(r22608, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r22610, mpfr_get_si(r22606, MPFR_RNDN) || mpfr_get_si(r22609, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r22613, r22604, r22612, MPFR_RNDN); mpfr_div(r22614, r22611, r22613, MPFR_RNDN); mpfr_div(r22615, r22611, r22604, MPFR_RNDN); mpfr_add(r22616, r22614, r22615, MPFR_RNDN); ; mpfr_pow(r22618, r22604, r22617, MPFR_RNDN); mpfr_div(r22619, r22611, r22618, MPFR_RNDN); mpfr_sub(r22620, r22616, r22619, MPFR_RNDN); mpfr_mul(r22621, r22604, r22604, MPFR_RNDN); mpfr_add(r22622, r22621, r22611, MPFR_RNDN); mpfr_div(r22623, r22604, r22622, MPFR_RNDN); if (mpfr_get_si(r22610, MPFR_RNDN)) { mpfr_set(r22624, r22620, MPFR_RNDN); } else { mpfr_set(r22624, r22623, MPFR_RNDN); }; return mpfr_get_d(r22624, MPFR_RNDN); } static mpfr_t r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635, r22636, r22637, r22638, r22639, r22640, r22641, r22642, r22643, r22644, r22645; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22625); mpfr_init_set_str(r22626, "-1335004178484.0273", 10, MPFR_RNDN); mpfr_init(r22627); mpfr_init_set_str(r22628, "891429.0791613997", 10, MPFR_RNDN); mpfr_init(r22629); mpfr_init(r22630); mpfr_init(r22631); mpfr_init_set_str(r22632, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22633, "5", 10, MPFR_RNDN); mpfr_init(r22634); mpfr_init(r22635); mpfr_init(r22636); mpfr_init(r22637); mpfr_init_set_str(r22638, "3", 10, MPFR_RNDN); mpfr_init(r22639); mpfr_init(r22640); mpfr_init(r22641); mpfr_init(r22642); mpfr_init(r22643); mpfr_init(r22644); mpfr_init(r22645); } double f_dm(double x) { mpfr_set_d(r22625, x, MPFR_RNDN); ; mpfr_set_si(r22627, mpfr_cmp(r22625, r22626) <= 0, MPFR_RNDN); ; mpfr_set_si(r22629, mpfr_cmp(r22625, r22628) <= 0, MPFR_RNDN); mpfr_set_si(r22630, !mpfr_get_si(r22629, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r22631, mpfr_get_si(r22627, MPFR_RNDN) || mpfr_get_si(r22630, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r22634, r22625, r22633, MPFR_RNDN); mpfr_div(r22635, r22632, r22634, MPFR_RNDN); mpfr_div(r22636, r22632, r22625, MPFR_RNDN); mpfr_add(r22637, r22635, r22636, MPFR_RNDN); ; mpfr_pow(r22639, r22625, r22638, MPFR_RNDN); mpfr_div(r22640, r22632, r22639, MPFR_RNDN); mpfr_sub(r22641, r22637, r22640, MPFR_RNDN); mpfr_mul(r22642, r22625, r22625, MPFR_RNDN); mpfr_add(r22643, r22642, r22632, MPFR_RNDN); mpfr_div(r22644, r22625, r22643, MPFR_RNDN); if (mpfr_get_si(r22631, MPFR_RNDN)) { mpfr_set(r22645, r22641, MPFR_RNDN); } else { mpfr_set(r22645, r22644, MPFR_RNDN); }; return mpfr_get_d(r22645, MPFR_RNDN); }