#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22650 = x; float r22651 = r22650 * r22650; float r22652 = 1; float r22653 = r22651 + r22652; float r22654 = r22650 / r22653; return r22654; } double f_id(double x) { double r22655 = x; double r22656 = r22655 * r22655; double r22657 = 1; double r22658 = r22656 + r22657; double r22659 = r22655 / r22658; return r22659; } double f_of(float x) { float r22660 = x; float r22661 = -7.340577108609785e+19; bool r22662 = r22660 <= r22661; float r22663 = 1; float r22664 = 5; float r22665 = pow(r22660, r22664); float r22666 = r22663 / r22665; float r22667 = r22663 / r22660; float r22668 = r22666 + r22667; float r22669 = 3; float r22670 = pow(r22660, r22669); float r22671 = r22663 / r22670; float r22672 = r22668 - r22671; float r22673 = 8160780.502918356; bool r22674 = r22660 <= r22673; float r22675 = r22660 * r22660; float r22676 = r22675 + r22663; float r22677 = r22660 / r22676; float r22678 = r22674 ? r22677 : r22672; float r22679 = r22662 ? r22672 : r22678; return r22679; } double f_od(double x) { double r22680 = x; double r22681 = -7.340577108609785e+19; bool r22682 = r22680 <= r22681; double r22683 = 1; double r22684 = 5; double r22685 = pow(r22680, r22684); double r22686 = r22683 / r22685; double r22687 = r22683 / r22680; double r22688 = r22686 + r22687; double r22689 = 3; double r22690 = pow(r22680, r22689); double r22691 = r22683 / r22690; double r22692 = r22688 - r22691; double r22693 = 8160780.502918356; bool r22694 = r22680 <= r22693; double r22695 = r22680 * r22680; double r22696 = r22695 + r22683; double r22697 = r22680 / r22696; double r22698 = r22694 ? r22697 : r22692; double r22699 = r22682 ? r22692 : r22698; return r22699; } 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 r22700, r22701, r22702, r22703, r22704; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22700); mpfr_init(r22701); mpfr_init_set_str(r22702, "1", 10, MPFR_RNDN); mpfr_init(r22703); mpfr_init(r22704); } double f_im(double x) { mpfr_set_d(r22700, x, MPFR_RNDN); mpfr_mul(r22701, r22700, r22700, MPFR_RNDN); ; mpfr_add(r22703, r22701, r22702, MPFR_RNDN); mpfr_div(r22704, r22700, r22703, MPFR_RNDN); return mpfr_get_d(r22704, MPFR_RNDN); } static mpfr_t r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22705); mpfr_init_set_str(r22706, "-7.340577108609785e+19", 10, MPFR_RNDN); mpfr_init(r22707); mpfr_init_set_str(r22708, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22709, "5", 10, MPFR_RNDN); mpfr_init(r22710); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init_set_str(r22714, "3", 10, MPFR_RNDN); mpfr_init(r22715); mpfr_init(r22716); mpfr_init(r22717); mpfr_init_set_str(r22718, "8160780.502918356", 10, MPFR_RNDN); mpfr_init(r22719); mpfr_init(r22720); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); } double f_fm(double x) { mpfr_set_d(r22705, x, MPFR_RNDN); ; mpfr_set_si(r22707, mpfr_cmp(r22705, r22706) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22710, r22705, r22709, MPFR_RNDN); mpfr_div(r22711, r22708, r22710, MPFR_RNDN); mpfr_div(r22712, r22708, r22705, MPFR_RNDN); mpfr_add(r22713, r22711, r22712, MPFR_RNDN); ; mpfr_pow(r22715, r22705, r22714, MPFR_RNDN); mpfr_div(r22716, r22708, r22715, MPFR_RNDN); mpfr_sub(r22717, r22713, r22716, MPFR_RNDN); ; mpfr_set_si(r22719, mpfr_cmp(r22705, r22718) <= 0, MPFR_RNDN); mpfr_mul(r22720, r22705, r22705, MPFR_RNDN); mpfr_add(r22721, r22720, r22708, MPFR_RNDN); mpfr_div(r22722, r22705, r22721, MPFR_RNDN); if (mpfr_get_si(r22719, MPFR_RNDN)) { mpfr_set(r22723, r22722, MPFR_RNDN); } else { mpfr_set(r22723, r22717, MPFR_RNDN); }; if (mpfr_get_si(r22707, MPFR_RNDN)) { mpfr_set(r22724, r22717, MPFR_RNDN); } else { mpfr_set(r22724, r22723, MPFR_RNDN); }; return mpfr_get_d(r22724, MPFR_RNDN); } static mpfr_t r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22725); mpfr_init_set_str(r22726, "-7.340577108609785e+19", 10, MPFR_RNDN); mpfr_init(r22727); mpfr_init_set_str(r22728, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22729, "5", 10, MPFR_RNDN); mpfr_init(r22730); mpfr_init(r22731); mpfr_init(r22732); mpfr_init(r22733); mpfr_init_set_str(r22734, "3", 10, MPFR_RNDN); mpfr_init(r22735); mpfr_init(r22736); mpfr_init(r22737); mpfr_init_set_str(r22738, "8160780.502918356", 10, MPFR_RNDN); mpfr_init(r22739); mpfr_init(r22740); mpfr_init(r22741); mpfr_init(r22742); mpfr_init(r22743); mpfr_init(r22744); } double f_dm(double x) { mpfr_set_d(r22725, x, MPFR_RNDN); ; mpfr_set_si(r22727, mpfr_cmp(r22725, r22726) <= 0, MPFR_RNDN); ; ; mpfr_pow(r22730, r22725, r22729, MPFR_RNDN); mpfr_div(r22731, r22728, r22730, MPFR_RNDN); mpfr_div(r22732, r22728, r22725, MPFR_RNDN); mpfr_add(r22733, r22731, r22732, MPFR_RNDN); ; mpfr_pow(r22735, r22725, r22734, MPFR_RNDN); mpfr_div(r22736, r22728, r22735, MPFR_RNDN); mpfr_sub(r22737, r22733, r22736, MPFR_RNDN); ; mpfr_set_si(r22739, mpfr_cmp(r22725, r22738) <= 0, MPFR_RNDN); mpfr_mul(r22740, r22725, r22725, MPFR_RNDN); mpfr_add(r22741, r22740, r22728, MPFR_RNDN); mpfr_div(r22742, r22725, r22741, MPFR_RNDN); if (mpfr_get_si(r22739, MPFR_RNDN)) { mpfr_set(r22743, r22742, MPFR_RNDN); } else { mpfr_set(r22743, r22737, MPFR_RNDN); }; if (mpfr_get_si(r22727, MPFR_RNDN)) { mpfr_set(r22744, r22737, MPFR_RNDN); } else { mpfr_set(r22744, r22743, MPFR_RNDN); }; return mpfr_get_d(r22744, MPFR_RNDN); }