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