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