#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r8483 = x; float r8484 = r8483 * r8483; float r8485 = 1.0f; float r8486 = r8484 + r8485; float r8487 = r8483 / r8486; return r8487; } double f_id(double x) { double r8488 = x; double r8489 = r8488 * r8488; double r8490 = 1.0; double r8491 = r8489 + r8490; double r8492 = r8488 / r8491; return r8492; } double f_of(float x) { float r8493 = x; float r8494 = -685695.7792154074f; bool r8495 = r8493 <= r8494; float r8496 = 1104.59988535652f; bool r8497 = r8493 <= r8496; float r8498 = !r8497; bool r8499 = r8495 || r8498; float r8500 = 1.0f; float r8501 = r8500 / r8493; float r8502 = 5.0f; float r8503 = pow(r8493, r8502); float r8504 = r8500 / r8503; float r8505 = r8501 + r8504; float r8506 = 3.0f; float r8507 = pow(r8493, r8506); float r8508 = r8500 / r8507; float r8509 = r8505 - r8508; float r8510 = r8493 * r8493; float r8511 = r8510 + r8500; float r8512 = r8493 / r8511; float r8513 = r8499 ? r8509 : r8512; return r8513; } double f_od(double x) { double r8514 = x; double r8515 = -685695.7792154074; bool r8516 = r8514 <= r8515; double r8517 = 1104.59988535652; bool r8518 = r8514 <= r8517; double r8519 = !r8518; bool r8520 = r8516 || r8519; double r8521 = 1.0; double r8522 = r8521 / r8514; double r8523 = 5.0; double r8524 = pow(r8514, r8523); double r8525 = r8521 / r8524; double r8526 = r8522 + r8525; double r8527 = 3.0; double r8528 = pow(r8514, r8527); double r8529 = r8521 / r8528; double r8530 = r8526 - r8529; double r8531 = r8514 * r8514; double r8532 = r8531 + r8521; double r8533 = r8514 / r8532; double r8534 = r8520 ? r8530 : r8533; return r8534; } 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 r8535, r8536, r8537, r8538, r8539; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r8535); mpfr_init(r8536); mpfr_init_set_str(r8537, "1", 10, MPFR_RNDN); mpfr_init(r8538); mpfr_init(r8539); } double f_im(double x) { mpfr_set_d(r8535, x, MPFR_RNDN); mpfr_mul(r8536, r8535, r8535, MPFR_RNDN); ; mpfr_add(r8538, r8536, r8537, MPFR_RNDN); mpfr_div(r8539, r8535, r8538, MPFR_RNDN); return mpfr_get_d(r8539, MPFR_RNDN); } static mpfr_t r8540, r8541, r8542, r8543, r8544, r8545, r8546, r8547, r8548, r8549, r8550, r8551, r8552, r8553, r8554, r8555, r8556, r8557, r8558, r8559, r8560; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r8540); mpfr_init_set_str(r8541, "-685695.7792154074", 10, MPFR_RNDN); mpfr_init(r8542); mpfr_init_set_str(r8543, "1104.59988535652", 10, MPFR_RNDN); mpfr_init(r8544); mpfr_init(r8545); mpfr_init(r8546); mpfr_init_set_str(r8547, "1", 10, MPFR_RNDN); mpfr_init(r8548); mpfr_init_set_str(r8549, "5", 10, MPFR_RNDN); mpfr_init(r8550); mpfr_init(r8551); mpfr_init(r8552); mpfr_init_set_str(r8553, "3", 10, MPFR_RNDN); mpfr_init(r8554); mpfr_init(r8555); mpfr_init(r8556); mpfr_init(r8557); mpfr_init(r8558); mpfr_init(r8559); mpfr_init(r8560); } double f_fm(double x) { mpfr_set_d(r8540, x, MPFR_RNDN); ; mpfr_set_si(r8542, mpfr_cmp(r8540, r8541) <= 0, MPFR_RNDN); ; mpfr_set_si(r8544, mpfr_cmp(r8540, r8543) <= 0, MPFR_RNDN); mpfr_set_si(r8545, !mpfr_get_si(r8544, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r8546, mpfr_get_si(r8542, MPFR_RNDN) || mpfr_get_si(r8545, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r8548, r8547, r8540, MPFR_RNDN); ; mpfr_pow(r8550, r8540, r8549, MPFR_RNDN); mpfr_div(r8551, r8547, r8550, MPFR_RNDN); mpfr_add(r8552, r8548, r8551, MPFR_RNDN); ; mpfr_pow(r8554, r8540, r8553, MPFR_RNDN); mpfr_div(r8555, r8547, r8554, MPFR_RNDN); mpfr_sub(r8556, r8552, r8555, MPFR_RNDN); mpfr_mul(r8557, r8540, r8540, MPFR_RNDN); mpfr_add(r8558, r8557, r8547, MPFR_RNDN); mpfr_div(r8559, r8540, r8558, MPFR_RNDN); if (mpfr_get_si(r8546, MPFR_RNDN)) { mpfr_set(r8560, r8556, MPFR_RNDN); } else { mpfr_set(r8560, r8559, MPFR_RNDN); }; return mpfr_get_d(r8560, MPFR_RNDN); } static mpfr_t r8561, r8562, r8563, r8564, r8565, r8566, r8567, r8568, r8569, r8570, r8571, r8572, r8573, r8574, r8575, r8576, r8577, r8578, r8579, r8580, r8581; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r8561); mpfr_init_set_str(r8562, "-685695.7792154074", 10, MPFR_RNDN); mpfr_init(r8563); mpfr_init_set_str(r8564, "1104.59988535652", 10, MPFR_RNDN); mpfr_init(r8565); mpfr_init(r8566); mpfr_init(r8567); mpfr_init_set_str(r8568, "1", 10, MPFR_RNDN); mpfr_init(r8569); mpfr_init_set_str(r8570, "5", 10, MPFR_RNDN); mpfr_init(r8571); mpfr_init(r8572); mpfr_init(r8573); mpfr_init_set_str(r8574, "3", 10, MPFR_RNDN); mpfr_init(r8575); mpfr_init(r8576); mpfr_init(r8577); mpfr_init(r8578); mpfr_init(r8579); mpfr_init(r8580); mpfr_init(r8581); } double f_dm(double x) { mpfr_set_d(r8561, x, MPFR_RNDN); ; mpfr_set_si(r8563, mpfr_cmp(r8561, r8562) <= 0, MPFR_RNDN); ; mpfr_set_si(r8565, mpfr_cmp(r8561, r8564) <= 0, MPFR_RNDN); mpfr_set_si(r8566, !mpfr_get_si(r8565, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r8567, mpfr_get_si(r8563, MPFR_RNDN) || mpfr_get_si(r8566, MPFR_RNDN), MPFR_RNDN); ; mpfr_div(r8569, r8568, r8561, MPFR_RNDN); ; mpfr_pow(r8571, r8561, r8570, MPFR_RNDN); mpfr_div(r8572, r8568, r8571, MPFR_RNDN); mpfr_add(r8573, r8569, r8572, MPFR_RNDN); ; mpfr_pow(r8575, r8561, r8574, MPFR_RNDN); mpfr_div(r8576, r8568, r8575, MPFR_RNDN); mpfr_sub(r8577, r8573, r8576, MPFR_RNDN); mpfr_mul(r8578, r8561, r8561, MPFR_RNDN); mpfr_add(r8579, r8578, r8568, MPFR_RNDN); mpfr_div(r8580, r8561, r8579, MPFR_RNDN); if (mpfr_get_si(r8567, MPFR_RNDN)) { mpfr_set(r8581, r8577, MPFR_RNDN); } else { mpfr_set(r8581, r8580, MPFR_RNDN); }; return mpfr_get_d(r8581, MPFR_RNDN); }