#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17250 = b; float r17251 = -r17250; float r17252 = r17250 * r17250; float r17253 = 4.0f; float r17254 = a; float r17255 = r17253 * r17254; float r17256 = c; float r17257 = r17255 * r17256; float r17258 = r17252 - r17257; float r17259 = sqrt(r17258); float r17260 = r17251 + r17259; float r17261 = 2.0f; float r17262 = r17261 * r17254; float r17263 = r17260 / r17262; return r17263; } double f_id(double a, double b, double c) { double r17264 = b; double r17265 = -r17264; double r17266 = r17264 * r17264; double r17267 = 4.0; double r17268 = a; double r17269 = r17267 * r17268; double r17270 = c; double r17271 = r17269 * r17270; double r17272 = r17266 - r17271; double r17273 = sqrt(r17272); double r17274 = r17265 + r17273; double r17275 = 2.0; double r17276 = r17275 * r17268; double r17277 = r17274 / r17276; return r17277; } double f_of(float a, float b, float c) { float r17278 = 1.0f; float r17279 = 2.0f; float r17280 = r17278 / r17279; float r17281 = 4.0f; float r17282 = c; float r17283 = r17281 * r17282; float r17284 = b; float r17285 = -r17284; float r17286 = r17284 * r17284; float r17287 = a; float r17288 = r17282 * r17287; float r17289 = r17288 * r17281; float r17290 = r17286 - r17289; float r17291 = sqrt(r17290); float r17292 = r17285 - r17291; float r17293 = r17283 / r17292; float r17294 = r17280 * r17293; return r17294; } double f_od(double a, double b, double c) { double r17295 = 1.0; double r17296 = 2.0; double r17297 = r17295 / r17296; double r17298 = 4.0; double r17299 = c; double r17300 = r17298 * r17299; double r17301 = b; double r17302 = -r17301; double r17303 = r17301 * r17301; double r17304 = a; double r17305 = r17299 * r17304; double r17306 = r17305 * r17298; double r17307 = r17303 - r17306; double r17308 = sqrt(r17307); double r17309 = r17302 - r17308; double r17310 = r17300 / r17309; double r17311 = r17297 * r17310; return r17311; } 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 r17312, r17313, r17314, r17315, r17316, r17317, r17318, r17319, r17320, r17321, r17322, r17323, r17324, r17325; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17312); mpfr_init(r17313); mpfr_init(r17314); mpfr_init_set_str(r17315, "4", 10, MPFR_RNDN); mpfr_init(r17316); mpfr_init(r17317); mpfr_init(r17318); mpfr_init(r17319); mpfr_init(r17320); mpfr_init(r17321); mpfr_init(r17322); mpfr_init_set_str(r17323, "2", 10, MPFR_RNDN); mpfr_init(r17324); mpfr_init(r17325); } double f_im(double a, double b, double c) { mpfr_set_d(r17312, b, MPFR_RNDN); mpfr_neg(r17313, r17312, MPFR_RNDN); mpfr_sqr(r17314, r17312, MPFR_RNDN); ; mpfr_set_d(r17316, a, MPFR_RNDN); mpfr_mul(r17317, r17315, r17316, MPFR_RNDN); mpfr_set_d(r17318, c, MPFR_RNDN); mpfr_mul(r17319, r17317, r17318, MPFR_RNDN); mpfr_sub(r17320, r17314, r17319, MPFR_RNDN); mpfr_sqrt(r17321, r17320, MPFR_RNDN); mpfr_add(r17322, r17313, r17321, MPFR_RNDN); ; mpfr_mul(r17324, r17323, r17316, MPFR_RNDN); mpfr_div(r17325, r17322, r17324, MPFR_RNDN); return mpfr_get_d(r17325, MPFR_RNDN); } static mpfr_t r17326, r17327, r17328, r17329, r17330, r17331, r17332, r17333, r17334, r17335, r17336, r17337, r17338, r17339, r17340, r17341, r17342; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17326, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17327, "2", 10, MPFR_RNDN); mpfr_init(r17328); mpfr_init_set_str(r17329, "4", 10, MPFR_RNDN); mpfr_init(r17330); mpfr_init(r17331); mpfr_init(r17332); mpfr_init(r17333); mpfr_init(r17334); mpfr_init(r17335); mpfr_init(r17336); mpfr_init(r17337); mpfr_init(r17338); mpfr_init(r17339); mpfr_init(r17340); mpfr_init(r17341); mpfr_init(r17342); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r17328, r17326, r17327, MPFR_RNDN); ; mpfr_set_d(r17330, c, MPFR_RNDN); mpfr_mul(r17331, r17329, r17330, MPFR_RNDN); mpfr_set_d(r17332, b, MPFR_RNDN); mpfr_neg(r17333, r17332, MPFR_RNDN); mpfr_sqr(r17334, r17332, MPFR_RNDN); mpfr_set_d(r17335, a, MPFR_RNDN); mpfr_mul(r17336, r17330, r17335, MPFR_RNDN); mpfr_mul(r17337, r17336, r17329, MPFR_RNDN); mpfr_sub(r17338, r17334, r17337, MPFR_RNDN); mpfr_sqrt(r17339, r17338, MPFR_RNDN); mpfr_sub(r17340, r17333, r17339, MPFR_RNDN); mpfr_div(r17341, r17331, r17340, MPFR_RNDN); mpfr_mul(r17342, r17328, r17341, MPFR_RNDN); return mpfr_get_d(r17342, MPFR_RNDN); } static mpfr_t r17343, r17344, r17345, r17346, r17347, r17348, r17349, r17350, r17351, r17352, r17353, r17354, r17355, r17356, r17357, r17358, r17359; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17343, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17344, "2", 10, MPFR_RNDN); mpfr_init(r17345); mpfr_init_set_str(r17346, "4", 10, MPFR_RNDN); mpfr_init(r17347); mpfr_init(r17348); mpfr_init(r17349); mpfr_init(r17350); mpfr_init(r17351); mpfr_init(r17352); mpfr_init(r17353); mpfr_init(r17354); mpfr_init(r17355); mpfr_init(r17356); mpfr_init(r17357); mpfr_init(r17358); mpfr_init(r17359); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r17345, r17343, r17344, MPFR_RNDN); ; mpfr_set_d(r17347, c, MPFR_RNDN); mpfr_mul(r17348, r17346, r17347, MPFR_RNDN); mpfr_set_d(r17349, b, MPFR_RNDN); mpfr_neg(r17350, r17349, MPFR_RNDN); mpfr_sqr(r17351, r17349, MPFR_RNDN); mpfr_set_d(r17352, a, MPFR_RNDN); mpfr_mul(r17353, r17347, r17352, MPFR_RNDN); mpfr_mul(r17354, r17353, r17346, MPFR_RNDN); mpfr_sub(r17355, r17351, r17354, MPFR_RNDN); mpfr_sqrt(r17356, r17355, MPFR_RNDN); mpfr_sub(r17357, r17350, r17356, MPFR_RNDN); mpfr_div(r17358, r17348, r17357, MPFR_RNDN); mpfr_mul(r17359, r17345, r17358, MPFR_RNDN); return mpfr_get_d(r17359, MPFR_RNDN); }