#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17206 = b; float r17207 = -r17206; float r17208 = r17206 * r17206; float r17209 = 4.0f; float r17210 = a; float r17211 = r17209 * r17210; float r17212 = c; float r17213 = r17211 * r17212; float r17214 = r17208 - r17213; float r17215 = sqrt(r17214); float r17216 = r17207 + r17215; float r17217 = 2.0f; float r17218 = r17217 * r17210; float r17219 = r17216 / r17218; return r17219; } double f_id(double a, double b, double c) { double r17220 = b; double r17221 = -r17220; double r17222 = r17220 * r17220; double r17223 = 4.0; double r17224 = a; double r17225 = r17223 * r17224; double r17226 = c; double r17227 = r17225 * r17226; double r17228 = r17222 - r17227; double r17229 = sqrt(r17228); double r17230 = r17221 + r17229; double r17231 = 2.0; double r17232 = r17231 * r17224; double r17233 = r17230 / r17232; return r17233; } double f_of(float a, float b, float c) { float r17234 = 4.0f; float r17235 = a; float r17236 = r17234 * r17235; float r17237 = c; float r17238 = r17236 * r17237; float r17239 = b; float r17240 = -r17239; float r17241 = r17239 * r17239; float r17242 = r17241 - r17238; float r17243 = sqrt(r17242); float r17244 = r17240 - r17243; float r17245 = r17238 / r17244; float r17246 = 2.0f; float r17247 = r17246 * r17235; float r17248 = r17245 / r17247; return r17248; } double f_od(double a, double b, double c) { double r17249 = 4.0; double r17250 = a; double r17251 = r17249 * r17250; double r17252 = c; double r17253 = r17251 * r17252; double r17254 = b; double r17255 = -r17254; double r17256 = r17254 * r17254; double r17257 = r17256 - r17253; double r17258 = sqrt(r17257); double r17259 = r17255 - r17258; double r17260 = r17253 / r17259; double r17261 = 2.0; double r17262 = r17261 * r17250; double r17263 = r17260 / r17262; return r17263; } 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 r17264, r17265, r17266, r17267, r17268, r17269, r17270, r17271, r17272, r17273, r17274, r17275, r17276, r17277; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17264); mpfr_init(r17265); mpfr_init(r17266); mpfr_init_set_str(r17267, "4", 10, MPFR_RNDN); mpfr_init(r17268); mpfr_init(r17269); mpfr_init(r17270); mpfr_init(r17271); mpfr_init(r17272); mpfr_init(r17273); mpfr_init(r17274); mpfr_init_set_str(r17275, "2", 10, MPFR_RNDN); mpfr_init(r17276); mpfr_init(r17277); } double f_im(double a, double b, double c) { mpfr_set_d(r17264, b, MPFR_RNDN); mpfr_neg(r17265, r17264, MPFR_RNDN); mpfr_sqr(r17266, r17264, MPFR_RNDN); ; mpfr_set_d(r17268, a, MPFR_RNDN); mpfr_mul(r17269, r17267, r17268, MPFR_RNDN); mpfr_set_d(r17270, c, MPFR_RNDN); mpfr_mul(r17271, r17269, r17270, MPFR_RNDN); mpfr_sub(r17272, r17266, r17271, MPFR_RNDN); mpfr_sqrt(r17273, r17272, MPFR_RNDN); mpfr_add(r17274, r17265, r17273, MPFR_RNDN); ; mpfr_mul(r17276, r17275, r17268, MPFR_RNDN); mpfr_div(r17277, r17274, r17276, MPFR_RNDN); return mpfr_get_d(r17277, MPFR_RNDN); } static mpfr_t r17278, r17279, r17280, r17281, r17282, r17283, r17284, r17285, r17286, r17287, r17288, r17289, r17290, r17291, r17292; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17278, "4", 10, MPFR_RNDN); mpfr_init(r17279); mpfr_init(r17280); mpfr_init(r17281); mpfr_init(r17282); mpfr_init(r17283); mpfr_init(r17284); mpfr_init(r17285); mpfr_init(r17286); mpfr_init(r17287); mpfr_init(r17288); mpfr_init(r17289); mpfr_init_set_str(r17290, "2", 10, MPFR_RNDN); mpfr_init(r17291); mpfr_init(r17292); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17279, a, MPFR_RNDN); mpfr_mul(r17280, r17278, r17279, MPFR_RNDN); mpfr_set_d(r17281, c, MPFR_RNDN); mpfr_mul(r17282, r17280, r17281, MPFR_RNDN); mpfr_set_d(r17283, b, MPFR_RNDN); mpfr_neg(r17284, r17283, MPFR_RNDN); mpfr_sqr(r17285, r17283, MPFR_RNDN); mpfr_sub(r17286, r17285, r17282, MPFR_RNDN); mpfr_sqrt(r17287, r17286, MPFR_RNDN); mpfr_sub(r17288, r17284, r17287, MPFR_RNDN); mpfr_div(r17289, r17282, r17288, MPFR_RNDN); ; mpfr_mul(r17291, r17290, r17279, MPFR_RNDN); mpfr_div(r17292, r17289, r17291, MPFR_RNDN); return mpfr_get_d(r17292, MPFR_RNDN); } static mpfr_t r17293, r17294, r17295, r17296, r17297, r17298, r17299, r17300, r17301, r17302, r17303, r17304, r17305, r17306, r17307; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17293, "4", 10, MPFR_RNDN); mpfr_init(r17294); mpfr_init(r17295); mpfr_init(r17296); mpfr_init(r17297); mpfr_init(r17298); mpfr_init(r17299); mpfr_init(r17300); mpfr_init(r17301); mpfr_init(r17302); mpfr_init(r17303); mpfr_init(r17304); mpfr_init_set_str(r17305, "2", 10, MPFR_RNDN); mpfr_init(r17306); mpfr_init(r17307); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17294, a, MPFR_RNDN); mpfr_mul(r17295, r17293, r17294, MPFR_RNDN); mpfr_set_d(r17296, c, MPFR_RNDN); mpfr_mul(r17297, r17295, r17296, MPFR_RNDN); mpfr_set_d(r17298, b, MPFR_RNDN); mpfr_neg(r17299, r17298, MPFR_RNDN); mpfr_sqr(r17300, r17298, MPFR_RNDN); mpfr_sub(r17301, r17300, r17297, MPFR_RNDN); mpfr_sqrt(r17302, r17301, MPFR_RNDN); mpfr_sub(r17303, r17299, r17302, MPFR_RNDN); mpfr_div(r17304, r17297, r17303, MPFR_RNDN); ; mpfr_mul(r17306, r17305, r17294, MPFR_RNDN); mpfr_div(r17307, r17304, r17306, MPFR_RNDN); return mpfr_get_d(r17307, MPFR_RNDN); }