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