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