#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r8617 = b; float r8618 = -r8617; float r8619 = r8617 * r8617; float r8620 = 4; float r8621 = a; float r8622 = r8620 * r8621; float r8623 = c; float r8624 = r8622 * r8623; float r8625 = r8619 - r8624; float r8626 = sqrt(r8625); float r8627 = r8618 + r8626; float r8628 = 2; float r8629 = r8628 * r8621; float r8630 = r8627 / r8629; return r8630; } double f_id(double a, double b, double c) { double r8631 = b; double r8632 = -r8631; double r8633 = r8631 * r8631; double r8634 = 4; double r8635 = a; double r8636 = r8634 * r8635; double r8637 = c; double r8638 = r8636 * r8637; double r8639 = r8633 - r8638; double r8640 = sqrt(r8639); double r8641 = r8632 + r8640; double r8642 = 2; double r8643 = r8642 * r8635; double r8644 = r8641 / r8643; return r8644; } double f_of(float a, float b, float c) { float r8645 = c; float r8646 = 2; float r8647 = 4; float r8648 = r8646 / r8647; float r8649 = r8645 / r8648; float r8650 = 1; float r8651 = b; float r8652 = -r8651; float r8653 = r8651 * r8651; float r8654 = a; float r8655 = r8647 * r8654; float r8656 = r8655 * r8645; float r8657 = r8653 - r8656; float r8658 = sqrt(r8657); float r8659 = r8652 - r8658; float r8660 = r8650 / r8659; float r8661 = r8649 * r8660; return r8661; } double f_od(double a, double b, double c) { double r8662 = c; double r8663 = 2; double r8664 = 4; double r8665 = r8663 / r8664; double r8666 = r8662 / r8665; double r8667 = 1; double r8668 = b; double r8669 = -r8668; double r8670 = r8668 * r8668; double r8671 = a; double r8672 = r8664 * r8671; double r8673 = r8672 * r8662; double r8674 = r8670 - r8673; double r8675 = sqrt(r8674); double r8676 = r8669 - r8675; double r8677 = r8667 / r8676; double r8678 = r8666 * r8677; return r8678; } 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 r8679, r8680, r8681, r8682, r8683, r8684, r8685, r8686, r8687, r8688, r8689, r8690, r8691, r8692; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r8679); mpfr_init(r8680); mpfr_init(r8681); mpfr_init_set_str(r8682, "4", 10, MPFR_RNDN); mpfr_init(r8683); mpfr_init(r8684); mpfr_init(r8685); mpfr_init(r8686); mpfr_init(r8687); mpfr_init(r8688); mpfr_init(r8689); mpfr_init_set_str(r8690, "2", 10, MPFR_RNDN); mpfr_init(r8691); mpfr_init(r8692); } double f_im(double a, double b, double c) { mpfr_set_d(r8679, b, MPFR_RNDN); mpfr_neg(r8680, r8679, MPFR_RNDN); mpfr_mul(r8681, r8679, r8679, MPFR_RNDN); ; mpfr_set_d(r8683, a, MPFR_RNDN); mpfr_mul(r8684, r8682, r8683, MPFR_RNDN); mpfr_set_d(r8685, c, MPFR_RNDN); mpfr_mul(r8686, r8684, r8685, MPFR_RNDN); mpfr_sub(r8687, r8681, r8686, MPFR_RNDN); mpfr_sqrt(r8688, r8687, MPFR_RNDN); mpfr_add(r8689, r8680, r8688, MPFR_RNDN); ; mpfr_mul(r8691, r8690, r8683, MPFR_RNDN); mpfr_div(r8692, r8689, r8691, MPFR_RNDN); return mpfr_get_d(r8692, MPFR_RNDN); } static mpfr_t r8693, r8694, r8695, r8696, r8697, r8698, r8699, r8700, r8701, r8702, r8703, r8704, r8705, r8706, r8707, r8708, r8709; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r8693); mpfr_init_set_str(r8694, "2", 10, MPFR_RNDN); mpfr_init_set_str(r8695, "4", 10, MPFR_RNDN); mpfr_init(r8696); mpfr_init(r8697); mpfr_init_set_str(r8698, "1", 10, MPFR_RNDN); mpfr_init(r8699); mpfr_init(r8700); mpfr_init(r8701); mpfr_init(r8702); mpfr_init(r8703); mpfr_init(r8704); mpfr_init(r8705); mpfr_init(r8706); mpfr_init(r8707); mpfr_init(r8708); mpfr_init(r8709); } double f_fm(double a, double b, double c) { mpfr_set_d(r8693, c, MPFR_RNDN); ; ; mpfr_div(r8696, r8694, r8695, MPFR_RNDN); mpfr_div(r8697, r8693, r8696, MPFR_RNDN); ; mpfr_set_d(r8699, b, MPFR_RNDN); mpfr_neg(r8700, r8699, MPFR_RNDN); mpfr_mul(r8701, r8699, r8699, MPFR_RNDN); mpfr_set_d(r8702, a, MPFR_RNDN); mpfr_mul(r8703, r8695, r8702, MPFR_RNDN); mpfr_mul(r8704, r8703, r8693, MPFR_RNDN); mpfr_sub(r8705, r8701, r8704, MPFR_RNDN); mpfr_sqrt(r8706, r8705, MPFR_RNDN); mpfr_sub(r8707, r8700, r8706, MPFR_RNDN); mpfr_div(r8708, r8698, r8707, MPFR_RNDN); mpfr_mul(r8709, r8697, r8708, MPFR_RNDN); return mpfr_get_d(r8709, MPFR_RNDN); } static mpfr_t r8710, r8711, r8712, r8713, r8714, r8715, r8716, r8717, r8718, r8719, r8720, r8721, r8722, r8723, r8724, r8725, r8726; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r8710); mpfr_init_set_str(r8711, "2", 10, MPFR_RNDN); mpfr_init_set_str(r8712, "4", 10, MPFR_RNDN); mpfr_init(r8713); mpfr_init(r8714); mpfr_init_set_str(r8715, "1", 10, MPFR_RNDN); mpfr_init(r8716); mpfr_init(r8717); mpfr_init(r8718); mpfr_init(r8719); mpfr_init(r8720); mpfr_init(r8721); mpfr_init(r8722); mpfr_init(r8723); mpfr_init(r8724); mpfr_init(r8725); mpfr_init(r8726); } double f_dm(double a, double b, double c) { mpfr_set_d(r8710, c, MPFR_RNDN); ; ; mpfr_div(r8713, r8711, r8712, MPFR_RNDN); mpfr_div(r8714, r8710, r8713, MPFR_RNDN); ; mpfr_set_d(r8716, b, MPFR_RNDN); mpfr_neg(r8717, r8716, MPFR_RNDN); mpfr_mul(r8718, r8716, r8716, MPFR_RNDN); mpfr_set_d(r8719, a, MPFR_RNDN); mpfr_mul(r8720, r8712, r8719, MPFR_RNDN); mpfr_mul(r8721, r8720, r8710, MPFR_RNDN); mpfr_sub(r8722, r8718, r8721, MPFR_RNDN); mpfr_sqrt(r8723, r8722, MPFR_RNDN); mpfr_sub(r8724, r8717, r8723, MPFR_RNDN); mpfr_div(r8725, r8715, r8724, MPFR_RNDN); mpfr_mul(r8726, r8714, r8725, MPFR_RNDN); return mpfr_get_d(r8726, MPFR_RNDN); }