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