#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r21626 = x; float r21627 = r21626 * r21626; float r21628 = 1; float r21629 = r21627 + r21628; float r21630 = r21626 / r21629; return r21630; } double f_id(double x) { double r21631 = x; double r21632 = r21631 * r21631; double r21633 = 1; double r21634 = r21632 + r21633; double r21635 = r21631 / r21634; return r21635; } double f_of(float x) { float r21636 = x; float r21637 = -76951863.5324409; bool r21638 = r21636 <= r21637; float r21639 = 1; float r21640 = 5; float r21641 = pow(r21636, r21640); float r21642 = r21639 / r21641; float r21643 = r21639 / r21636; float r21644 = r21642 + r21643; float r21645 = 3; float r21646 = pow(r21636, r21645); float r21647 = r21639 / r21646; float r21648 = r21644 - r21647; float r21649 = 35688.954472079415; bool r21650 = r21636 <= r21649; float r21651 = fma(r21636, r21636, r21639); float r21652 = sqrt(r21651); float r21653 = r21636 / r21652; float r21654 = r21653 / r21652; float r21655 = r21650 ? r21654 : r21648; float r21656 = r21638 ? r21648 : r21655; return r21656; } double f_od(double x) { double r21657 = x; double r21658 = -76951863.5324409; bool r21659 = r21657 <= r21658; double r21660 = 1; double r21661 = 5; double r21662 = pow(r21657, r21661); double r21663 = r21660 / r21662; double r21664 = r21660 / r21657; double r21665 = r21663 + r21664; double r21666 = 3; double r21667 = pow(r21657, r21666); double r21668 = r21660 / r21667; double r21669 = r21665 - r21668; double r21670 = 35688.954472079415; bool r21671 = r21657 <= r21670; double r21672 = fma(r21657, r21657, r21660); double r21673 = sqrt(r21672); double r21674 = r21657 / r21673; double r21675 = r21674 / r21673; double r21676 = r21671 ? r21675 : r21669; double r21677 = r21659 ? r21669 : r21676; return r21677; } 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 r21678, r21679, r21680, r21681, r21682; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21678); mpfr_init(r21679); mpfr_init_set_str(r21680, "1", 10, MPFR_RNDN); mpfr_init(r21681); mpfr_init(r21682); } double f_im(double x) { mpfr_set_d(r21678, x, MPFR_RNDN); mpfr_mul(r21679, r21678, r21678, MPFR_RNDN); ; mpfr_add(r21681, r21679, r21680, MPFR_RNDN); mpfr_div(r21682, r21678, r21681, MPFR_RNDN); return mpfr_get_d(r21682, MPFR_RNDN); } static mpfr_t r21683, r21684, r21685, r21686, r21687, r21688, r21689, r21690, r21691, r21692, r21693, r21694, r21695, r21696, r21697, r21698, r21699, r21700, r21701, r21702, r21703; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21683); mpfr_init_set_str(r21684, "-76951863.5324409", 10, MPFR_RNDN); mpfr_init(r21685); mpfr_init_set_str(r21686, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21687, "5", 10, MPFR_RNDN); mpfr_init(r21688); mpfr_init(r21689); mpfr_init(r21690); mpfr_init(r21691); mpfr_init_set_str(r21692, "3", 10, MPFR_RNDN); mpfr_init(r21693); mpfr_init(r21694); mpfr_init(r21695); mpfr_init_set_str(r21696, "35688.954472079415", 10, MPFR_RNDN); mpfr_init(r21697); mpfr_init(r21698); mpfr_init(r21699); mpfr_init(r21700); mpfr_init(r21701); mpfr_init(r21702); mpfr_init(r21703); } double f_fm(double x) { mpfr_set_d(r21683, x, MPFR_RNDN); ; mpfr_set_si(r21685, mpfr_cmp(r21683, r21684) <= 0, MPFR_RNDN); ; ; mpfr_pow(r21688, r21683, r21687, MPFR_RNDN); mpfr_div(r21689, r21686, r21688, MPFR_RNDN); mpfr_div(r21690, r21686, r21683, MPFR_RNDN); mpfr_add(r21691, r21689, r21690, MPFR_RNDN); ; mpfr_pow(r21693, r21683, r21692, MPFR_RNDN); mpfr_div(r21694, r21686, r21693, MPFR_RNDN); mpfr_sub(r21695, r21691, r21694, MPFR_RNDN); ; mpfr_set_si(r21697, mpfr_cmp(r21683, r21696) <= 0, MPFR_RNDN); mpfr_fma(r21698, r21683, r21683, r21686, MPFR_RNDN); mpfr_sqrt(r21699, r21698, MPFR_RNDN); mpfr_div(r21700, r21683, r21699, MPFR_RNDN); mpfr_div(r21701, r21700, r21699, MPFR_RNDN); if (mpfr_get_si(r21697, MPFR_RNDN)) { mpfr_set(r21702, r21701, MPFR_RNDN); } else { mpfr_set(r21702, r21695, MPFR_RNDN); }; if (mpfr_get_si(r21685, MPFR_RNDN)) { mpfr_set(r21703, r21695, MPFR_RNDN); } else { mpfr_set(r21703, r21702, MPFR_RNDN); }; return mpfr_get_d(r21703, MPFR_RNDN); } static mpfr_t r21704, r21705, r21706, r21707, r21708, r21709, r21710, r21711, r21712, r21713, r21714, r21715, r21716, r21717, r21718, r21719, r21720, r21721, r21722, r21723, r21724; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21704); mpfr_init_set_str(r21705, "-76951863.5324409", 10, MPFR_RNDN); mpfr_init(r21706); mpfr_init_set_str(r21707, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21708, "5", 10, MPFR_RNDN); mpfr_init(r21709); mpfr_init(r21710); mpfr_init(r21711); mpfr_init(r21712); mpfr_init_set_str(r21713, "3", 10, MPFR_RNDN); mpfr_init(r21714); mpfr_init(r21715); mpfr_init(r21716); mpfr_init_set_str(r21717, "35688.954472079415", 10, MPFR_RNDN); mpfr_init(r21718); mpfr_init(r21719); mpfr_init(r21720); mpfr_init(r21721); mpfr_init(r21722); mpfr_init(r21723); mpfr_init(r21724); } double f_dm(double x) { mpfr_set_d(r21704, x, MPFR_RNDN); ; mpfr_set_si(r21706, mpfr_cmp(r21704, r21705) <= 0, MPFR_RNDN); ; ; mpfr_pow(r21709, r21704, r21708, MPFR_RNDN); mpfr_div(r21710, r21707, r21709, MPFR_RNDN); mpfr_div(r21711, r21707, r21704, MPFR_RNDN); mpfr_add(r21712, r21710, r21711, MPFR_RNDN); ; mpfr_pow(r21714, r21704, r21713, MPFR_RNDN); mpfr_div(r21715, r21707, r21714, MPFR_RNDN); mpfr_sub(r21716, r21712, r21715, MPFR_RNDN); ; mpfr_set_si(r21718, mpfr_cmp(r21704, r21717) <= 0, MPFR_RNDN); mpfr_fma(r21719, r21704, r21704, r21707, MPFR_RNDN); mpfr_sqrt(r21720, r21719, MPFR_RNDN); mpfr_div(r21721, r21704, r21720, MPFR_RNDN); mpfr_div(r21722, r21721, r21720, MPFR_RNDN); if (mpfr_get_si(r21718, MPFR_RNDN)) { mpfr_set(r21723, r21722, MPFR_RNDN); } else { mpfr_set(r21723, r21716, MPFR_RNDN); }; if (mpfr_get_si(r21706, MPFR_RNDN)) { mpfr_set(r21724, r21716, MPFR_RNDN); } else { mpfr_set(r21724, r21723, MPFR_RNDN); }; return mpfr_get_d(r21724, MPFR_RNDN); }