#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r21098 = b; float r21099 = -r21098; float r21100 = r21098 * r21098; float r21101 = 4; float r21102 = a; float r21103 = r21101 * r21102; float r21104 = c; float r21105 = r21103 * r21104; float r21106 = r21100 - r21105; float r21107 = sqrt(r21106); float r21108 = r21099 + r21107; float r21109 = 2; float r21110 = r21109 * r21102; float r21111 = r21108 / r21110; return r21111; } double f_id(double a, double b, double c) { double r21112 = b; double r21113 = -r21112; double r21114 = r21112 * r21112; double r21115 = 4; double r21116 = a; double r21117 = r21115 * r21116; double r21118 = c; double r21119 = r21117 * r21118; double r21120 = r21114 - r21119; double r21121 = sqrt(r21120); double r21122 = r21113 + r21121; double r21123 = 2; double r21124 = r21123 * r21116; double r21125 = r21122 / r21124; return r21125; } double f_of(float a, float b, float c) { float r21126 = 1; float r21127 = 2; float r21128 = r21126 / r21127; float r21129 = c; float r21130 = 4; float r21131 = r21129 * r21130; float r21132 = b; float r21133 = -r21132; float r21134 = r21132 * r21132; float r21135 = a; float r21136 = r21130 * r21135; float r21137 = r21129 * r21136; float r21138 = r21134 - r21137; float r21139 = sqrt(r21138); float r21140 = r21133 - r21139; float r21141 = r21131 / r21140; float r21142 = r21128 * r21141; return r21142; } double f_od(double a, double b, double c) { double r21143 = 1; double r21144 = 2; double r21145 = r21143 / r21144; double r21146 = c; double r21147 = 4; double r21148 = r21146 * r21147; double r21149 = b; double r21150 = -r21149; double r21151 = r21149 * r21149; double r21152 = a; double r21153 = r21147 * r21152; double r21154 = r21146 * r21153; double r21155 = r21151 - r21154; double r21156 = sqrt(r21155); double r21157 = r21150 - r21156; double r21158 = r21148 / r21157; double r21159 = r21145 * r21158; return r21159; } 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 r21160, r21161, r21162, r21163, r21164, r21165, r21166, r21167, r21168, r21169, r21170, r21171, r21172, r21173; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21160); mpfr_init(r21161); mpfr_init(r21162); mpfr_init_set_str(r21163, "4", 10, MPFR_RNDN); mpfr_init(r21164); mpfr_init(r21165); mpfr_init(r21166); mpfr_init(r21167); mpfr_init(r21168); mpfr_init(r21169); mpfr_init(r21170); mpfr_init_set_str(r21171, "2", 10, MPFR_RNDN); mpfr_init(r21172); mpfr_init(r21173); } double f_im(double a, double b, double c) { mpfr_set_d(r21160, b, MPFR_RNDN); mpfr_neg(r21161, r21160, MPFR_RNDN); mpfr_mul(r21162, r21160, r21160, MPFR_RNDN); ; mpfr_set_d(r21164, a, MPFR_RNDN); mpfr_mul(r21165, r21163, r21164, MPFR_RNDN); mpfr_set_d(r21166, c, MPFR_RNDN); mpfr_mul(r21167, r21165, r21166, MPFR_RNDN); mpfr_sub(r21168, r21162, r21167, MPFR_RNDN); mpfr_sqrt(r21169, r21168, MPFR_RNDN); mpfr_add(r21170, r21161, r21169, MPFR_RNDN); ; mpfr_mul(r21172, r21171, r21164, MPFR_RNDN); mpfr_div(r21173, r21170, r21172, MPFR_RNDN); return mpfr_get_d(r21173, MPFR_RNDN); } static mpfr_t r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21174, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21175, "2", 10, MPFR_RNDN); mpfr_init(r21176); mpfr_init(r21177); mpfr_init_set_str(r21178, "4", 10, MPFR_RNDN); mpfr_init(r21179); mpfr_init(r21180); mpfr_init(r21181); mpfr_init(r21182); mpfr_init(r21183); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r21176, r21174, r21175, MPFR_RNDN); mpfr_set_d(r21177, c, MPFR_RNDN); ; mpfr_mul(r21179, r21177, r21178, MPFR_RNDN); mpfr_set_d(r21180, b, MPFR_RNDN); mpfr_neg(r21181, r21180, MPFR_RNDN); mpfr_mul(r21182, r21180, r21180, MPFR_RNDN); mpfr_set_d(r21183, a, MPFR_RNDN); mpfr_mul(r21184, r21178, r21183, MPFR_RNDN); mpfr_mul(r21185, r21177, r21184, MPFR_RNDN); mpfr_sub(r21186, r21182, r21185, MPFR_RNDN); mpfr_sqrt(r21187, r21186, MPFR_RNDN); mpfr_sub(r21188, r21181, r21187, MPFR_RNDN); mpfr_div(r21189, r21179, r21188, MPFR_RNDN); mpfr_mul(r21190, r21176, r21189, MPFR_RNDN); return mpfr_get_d(r21190, MPFR_RNDN); } static mpfr_t r21191, r21192, r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21191, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21192, "2", 10, MPFR_RNDN); mpfr_init(r21193); mpfr_init(r21194); mpfr_init_set_str(r21195, "4", 10, MPFR_RNDN); mpfr_init(r21196); mpfr_init(r21197); mpfr_init(r21198); mpfr_init(r21199); mpfr_init(r21200); mpfr_init(r21201); mpfr_init(r21202); mpfr_init(r21203); mpfr_init(r21204); mpfr_init(r21205); mpfr_init(r21206); mpfr_init(r21207); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r21193, r21191, r21192, MPFR_RNDN); mpfr_set_d(r21194, c, MPFR_RNDN); ; mpfr_mul(r21196, r21194, r21195, MPFR_RNDN); mpfr_set_d(r21197, b, MPFR_RNDN); mpfr_neg(r21198, r21197, MPFR_RNDN); mpfr_mul(r21199, r21197, r21197, MPFR_RNDN); mpfr_set_d(r21200, a, MPFR_RNDN); mpfr_mul(r21201, r21195, r21200, MPFR_RNDN); mpfr_mul(r21202, r21194, r21201, MPFR_RNDN); mpfr_sub(r21203, r21199, r21202, MPFR_RNDN); mpfr_sqrt(r21204, r21203, MPFR_RNDN); mpfr_sub(r21205, r21198, r21204, MPFR_RNDN); mpfr_div(r21206, r21196, r21205, MPFR_RNDN); mpfr_mul(r21207, r21193, r21206, MPFR_RNDN); return mpfr_get_d(r21207, MPFR_RNDN); }