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