#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 = b; float r17252 = -r17251; float r17253 = r17251 * r17251; float r17254 = 4.0f; float r17255 = a; float r17256 = r17254 * r17255; float r17257 = c; float r17258 = r17256 * r17257; float r17259 = r17253 - r17258; float r17260 = sqrt(r17259); float r17261 = r17252 - r17260; float r17262 = 2.0f; float r17263 = r17262 / r17254; float r17264 = r17263 / r17257; float r17265 = r17261 * r17264; float r17266 = r17250 / r17265; return r17266; } double f_od(double a, double b, double c) { double r17267 = 1.0; double r17268 = b; double r17269 = -r17268; double r17270 = r17268 * r17268; double r17271 = 4.0; double r17272 = a; double r17273 = r17271 * r17272; double r17274 = c; double r17275 = r17273 * r17274; double r17276 = r17270 - r17275; double r17277 = sqrt(r17276); double r17278 = r17269 - r17277; double r17279 = 2.0; double r17280 = r17279 / r17271; double r17281 = r17280 / r17274; double r17282 = r17278 * r17281; double r17283 = r17267 / 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(r17299); mpfr_init(r17300); mpfr_init(r17301); mpfr_init_set_str(r17302, "4", 10, MPFR_RNDN); mpfr_init(r17303); mpfr_init(r17304); mpfr_init(r17305); mpfr_init(r17306); mpfr_init(r17307); mpfr_init(r17308); mpfr_init(r17309); mpfr_init_set_str(r17310, "2", 10, MPFR_RNDN); mpfr_init(r17311); mpfr_init(r17312); mpfr_init(r17313); mpfr_init(r17314); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17299, b, MPFR_RNDN); mpfr_neg(r17300, r17299, MPFR_RNDN); mpfr_sqr(r17301, r17299, MPFR_RNDN); ; mpfr_set_d(r17303, a, MPFR_RNDN); mpfr_mul(r17304, r17302, r17303, MPFR_RNDN); mpfr_set_d(r17305, c, MPFR_RNDN); mpfr_mul(r17306, r17304, r17305, MPFR_RNDN); mpfr_sub(r17307, r17301, r17306, MPFR_RNDN); mpfr_sqrt(r17308, r17307, MPFR_RNDN); mpfr_sub(r17309, r17300, r17308, MPFR_RNDN); ; mpfr_div(r17311, r17310, r17302, MPFR_RNDN); mpfr_div(r17312, r17311, r17305, MPFR_RNDN); mpfr_mul(r17313, r17309, r17312, MPFR_RNDN); mpfr_div(r17314, r17298, 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(r17316); mpfr_init(r17317); mpfr_init(r17318); mpfr_init_set_str(r17319, "4", 10, MPFR_RNDN); mpfr_init(r17320); mpfr_init(r17321); mpfr_init(r17322); mpfr_init(r17323); mpfr_init(r17324); mpfr_init(r17325); mpfr_init(r17326); mpfr_init_set_str(r17327, "2", 10, MPFR_RNDN); mpfr_init(r17328); mpfr_init(r17329); mpfr_init(r17330); mpfr_init(r17331); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17316, b, MPFR_RNDN); mpfr_neg(r17317, r17316, MPFR_RNDN); mpfr_sqr(r17318, r17316, MPFR_RNDN); ; mpfr_set_d(r17320, a, MPFR_RNDN); mpfr_mul(r17321, r17319, r17320, MPFR_RNDN); mpfr_set_d(r17322, c, MPFR_RNDN); mpfr_mul(r17323, r17321, r17322, MPFR_RNDN); mpfr_sub(r17324, r17318, r17323, MPFR_RNDN); mpfr_sqrt(r17325, r17324, MPFR_RNDN); mpfr_sub(r17326, r17317, r17325, MPFR_RNDN); ; mpfr_div(r17328, r17327, r17319, MPFR_RNDN); mpfr_div(r17329, r17328, r17322, MPFR_RNDN); mpfr_mul(r17330, r17326, r17329, MPFR_RNDN); mpfr_div(r17331, r17315, r17330, MPFR_RNDN); return mpfr_get_d(r17331, MPFR_RNDN); }