#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r21679 = b; float r21680 = -r21679; float r21681 = r21679 * r21679; float r21682 = 4; float r21683 = a; float r21684 = r21682 * r21683; float r21685 = c; float r21686 = r21684 * r21685; float r21687 = r21681 - r21686; float r21688 = sqrt(r21687); float r21689 = r21680 + r21688; float r21690 = 2; float r21691 = r21690 * r21683; float r21692 = r21689 / r21691; return r21692; } double f_id(double a, double b, double c) { double r21693 = b; double r21694 = -r21693; double r21695 = r21693 * r21693; double r21696 = 4; double r21697 = a; double r21698 = r21696 * r21697; double r21699 = c; double r21700 = r21698 * r21699; double r21701 = r21695 - r21700; double r21702 = sqrt(r21701); double r21703 = r21694 + r21702; double r21704 = 2; double r21705 = r21704 * r21697; double r21706 = r21703 / r21705; return r21706; } double f_of(float a, float b, float c) { float r21707 = 1; float r21708 = 2; float r21709 = r21707 / r21708; float r21710 = 4; float r21711 = c; float r21712 = r21710 * r21711; float r21713 = b; float r21714 = -r21713; float r21715 = -r21711; float r21716 = a; float r21717 = r21710 * r21716; float r21718 = r21713 * r21713; float r21719 = fma(r21715, r21717, r21718); float r21720 = sqrt(r21719); float r21721 = r21714 - r21720; float r21722 = r21712 / r21721; float r21723 = r21709 * r21722; return r21723; } double f_od(double a, double b, double c) { double r21724 = 1; double r21725 = 2; double r21726 = r21724 / r21725; double r21727 = 4; double r21728 = c; double r21729 = r21727 * r21728; double r21730 = b; double r21731 = -r21730; double r21732 = -r21728; double r21733 = a; double r21734 = r21727 * r21733; double r21735 = r21730 * r21730; double r21736 = fma(r21732, r21734, r21735); double r21737 = sqrt(r21736); double r21738 = r21731 - r21737; double r21739 = r21729 / r21738; double r21740 = r21726 * r21739; return r21740; } 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 r21741, r21742, r21743, r21744, r21745, r21746, r21747, r21748, r21749, r21750, r21751, r21752, r21753, r21754; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r21741); mpfr_init(r21742); mpfr_init(r21743); mpfr_init_set_str(r21744, "4", 10, MPFR_RNDN); mpfr_init(r21745); mpfr_init(r21746); mpfr_init(r21747); mpfr_init(r21748); mpfr_init(r21749); mpfr_init(r21750); mpfr_init(r21751); mpfr_init_set_str(r21752, "2", 10, MPFR_RNDN); mpfr_init(r21753); mpfr_init(r21754); } double f_im(double a, double b, double c) { mpfr_set_d(r21741, b, MPFR_RNDN); mpfr_neg(r21742, r21741, MPFR_RNDN); mpfr_mul(r21743, r21741, r21741, MPFR_RNDN); ; mpfr_set_d(r21745, a, MPFR_RNDN); mpfr_mul(r21746, r21744, r21745, MPFR_RNDN); mpfr_set_d(r21747, c, MPFR_RNDN); mpfr_mul(r21748, r21746, r21747, MPFR_RNDN); mpfr_sub(r21749, r21743, r21748, MPFR_RNDN); mpfr_sqrt(r21750, r21749, MPFR_RNDN); mpfr_add(r21751, r21742, r21750, MPFR_RNDN); ; mpfr_mul(r21753, r21752, r21745, MPFR_RNDN); mpfr_div(r21754, r21751, r21753, MPFR_RNDN); return mpfr_get_d(r21754, MPFR_RNDN); } static mpfr_t r21755, r21756, r21757, r21758, r21759, r21760, r21761, r21762, r21763, r21764, r21765, r21766, r21767, r21768, r21769, r21770, r21771; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init_set_str(r21755, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21756, "2", 10, MPFR_RNDN); mpfr_init(r21757); mpfr_init_set_str(r21758, "4", 10, MPFR_RNDN); mpfr_init(r21759); mpfr_init(r21760); mpfr_init(r21761); mpfr_init(r21762); mpfr_init(r21763); mpfr_init(r21764); mpfr_init(r21765); mpfr_init(r21766); mpfr_init(r21767); mpfr_init(r21768); mpfr_init(r21769); mpfr_init(r21770); mpfr_init(r21771); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r21757, r21755, r21756, MPFR_RNDN); ; mpfr_set_d(r21759, c, MPFR_RNDN); mpfr_mul(r21760, r21758, r21759, MPFR_RNDN); mpfr_set_d(r21761, b, MPFR_RNDN); mpfr_neg(r21762, r21761, MPFR_RNDN); mpfr_neg(r21763, r21759, MPFR_RNDN); mpfr_set_d(r21764, a, MPFR_RNDN); mpfr_mul(r21765, r21758, r21764, MPFR_RNDN); mpfr_mul(r21766, r21761, r21761, MPFR_RNDN); mpfr_fma(r21767, r21763, r21765, r21766, MPFR_RNDN); mpfr_sqrt(r21768, r21767, MPFR_RNDN); mpfr_sub(r21769, r21762, r21768, MPFR_RNDN); mpfr_div(r21770, r21760, r21769, MPFR_RNDN); mpfr_mul(r21771, r21757, r21770, MPFR_RNDN); return mpfr_get_d(r21771, MPFR_RNDN); } static mpfr_t r21772, r21773, r21774, r21775, r21776, r21777, r21778, r21779, r21780, r21781, r21782, r21783, r21784, r21785, r21786, r21787, r21788; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init_set_str(r21772, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21773, "2", 10, MPFR_RNDN); mpfr_init(r21774); mpfr_init_set_str(r21775, "4", 10, MPFR_RNDN); mpfr_init(r21776); mpfr_init(r21777); mpfr_init(r21778); mpfr_init(r21779); mpfr_init(r21780); mpfr_init(r21781); mpfr_init(r21782); mpfr_init(r21783); mpfr_init(r21784); mpfr_init(r21785); mpfr_init(r21786); mpfr_init(r21787); mpfr_init(r21788); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r21774, r21772, r21773, MPFR_RNDN); ; mpfr_set_d(r21776, c, MPFR_RNDN); mpfr_mul(r21777, r21775, r21776, MPFR_RNDN); mpfr_set_d(r21778, b, MPFR_RNDN); mpfr_neg(r21779, r21778, MPFR_RNDN); mpfr_neg(r21780, r21776, MPFR_RNDN); mpfr_set_d(r21781, a, MPFR_RNDN); mpfr_mul(r21782, r21775, r21781, MPFR_RNDN); mpfr_mul(r21783, r21778, r21778, MPFR_RNDN); mpfr_fma(r21784, r21780, r21782, r21783, MPFR_RNDN); mpfr_sqrt(r21785, r21784, MPFR_RNDN); mpfr_sub(r21786, r21779, r21785, MPFR_RNDN); mpfr_div(r21787, r21777, r21786, MPFR_RNDN); mpfr_mul(r21788, r21774, r21787, MPFR_RNDN); return mpfr_get_d(r21788, MPFR_RNDN); }