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