#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17282 = b; float r17283 = -r17282; float r17284 = r17282 * r17282; float r17285 = 4.0f; float r17286 = a; float r17287 = r17285 * r17286; float r17288 = c; float r17289 = r17287 * r17288; float r17290 = r17284 - r17289; float r17291 = sqrt(r17290); float r17292 = r17283 + r17291; float r17293 = 2.0f; float r17294 = r17293 * r17286; float r17295 = r17292 / r17294; return r17295; } double f_id(double a, double b, double c) { double r17296 = b; double r17297 = -r17296; double r17298 = r17296 * r17296; double r17299 = 4.0; double r17300 = a; double r17301 = r17299 * r17300; double r17302 = c; double r17303 = r17301 * r17302; double r17304 = r17298 - r17303; double r17305 = sqrt(r17304); double r17306 = r17297 + r17305; double r17307 = 2.0; double r17308 = r17307 * r17300; double r17309 = r17306 / r17308; return r17309; } double f_of(float a, float b, float c) { float r17310 = 1.0f; float r17311 = 2.0f; float r17312 = r17310 / r17311; float r17313 = b; float r17314 = -r17313; float r17315 = r17313 * r17313; float r17316 = a; float r17317 = c; float r17318 = 4.0f; float r17319 = r17317 * r17318; float r17320 = r17316 * r17319; float r17321 = r17315 - r17320; float r17322 = sqrt(r17321); float r17323 = r17314 + r17322; float r17324 = r17323 / r17316; float r17325 = r17312 * r17324; return r17325; } double f_od(double a, double b, double c) { double r17326 = 1.0; double r17327 = 2.0; double r17328 = r17326 / r17327; double r17329 = b; double r17330 = -r17329; double r17331 = r17329 * r17329; double r17332 = a; double r17333 = c; double r17334 = 4.0; double r17335 = r17333 * r17334; double r17336 = r17332 * r17335; double r17337 = r17331 - r17336; double r17338 = sqrt(r17337); double r17339 = r17330 + r17338; double r17340 = r17339 / r17332; double r17341 = r17328 * r17340; return r17341; } 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 r17342, r17343, r17344, r17345, r17346, r17347, r17348, r17349, r17350, r17351, r17352, r17353, r17354, r17355; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17342); mpfr_init(r17343); mpfr_init(r17344); mpfr_init_set_str(r17345, "4", 10, MPFR_RNDN); mpfr_init(r17346); mpfr_init(r17347); mpfr_init(r17348); mpfr_init(r17349); mpfr_init(r17350); mpfr_init(r17351); mpfr_init(r17352); mpfr_init_set_str(r17353, "2", 10, MPFR_RNDN); mpfr_init(r17354); mpfr_init(r17355); } double f_im(double a, double b, double c) { mpfr_set_d(r17342, b, MPFR_RNDN); mpfr_neg(r17343, r17342, MPFR_RNDN); mpfr_sqr(r17344, r17342, MPFR_RNDN); ; mpfr_set_d(r17346, a, MPFR_RNDN); mpfr_mul(r17347, r17345, r17346, MPFR_RNDN); mpfr_set_d(r17348, c, MPFR_RNDN); mpfr_mul(r17349, r17347, r17348, MPFR_RNDN); mpfr_sub(r17350, r17344, r17349, MPFR_RNDN); mpfr_sqrt(r17351, r17350, MPFR_RNDN); mpfr_add(r17352, r17343, r17351, MPFR_RNDN); ; mpfr_mul(r17354, r17353, r17346, MPFR_RNDN); mpfr_div(r17355, r17352, r17354, MPFR_RNDN); return mpfr_get_d(r17355, MPFR_RNDN); } static mpfr_t r17356, r17357, r17358, r17359, r17360, r17361, r17362, r17363, r17364, r17365, r17366, r17367, r17368, r17369, r17370, r17371; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17356, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17357, "2", 10, MPFR_RNDN); mpfr_init(r17358); mpfr_init(r17359); mpfr_init(r17360); mpfr_init(r17361); mpfr_init(r17362); mpfr_init(r17363); mpfr_init_set_str(r17364, "4", 10, MPFR_RNDN); mpfr_init(r17365); mpfr_init(r17366); mpfr_init(r17367); mpfr_init(r17368); mpfr_init(r17369); mpfr_init(r17370); mpfr_init(r17371); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r17358, r17356, r17357, MPFR_RNDN); mpfr_set_d(r17359, b, MPFR_RNDN); mpfr_neg(r17360, r17359, MPFR_RNDN); mpfr_sqr(r17361, r17359, MPFR_RNDN); mpfr_set_d(r17362, a, MPFR_RNDN); mpfr_set_d(r17363, c, MPFR_RNDN); ; mpfr_mul(r17365, r17363, r17364, MPFR_RNDN); mpfr_mul(r17366, r17362, r17365, MPFR_RNDN); mpfr_sub(r17367, r17361, r17366, MPFR_RNDN); mpfr_sqrt(r17368, r17367, MPFR_RNDN); mpfr_add(r17369, r17360, r17368, MPFR_RNDN); mpfr_div(r17370, r17369, r17362, MPFR_RNDN); mpfr_mul(r17371, r17358, r17370, MPFR_RNDN); return mpfr_get_d(r17371, MPFR_RNDN); } static mpfr_t r17372, r17373, r17374, r17375, r17376, r17377, r17378, r17379, r17380, r17381, r17382, r17383, r17384, r17385, r17386, r17387; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17372, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17373, "2", 10, MPFR_RNDN); mpfr_init(r17374); mpfr_init(r17375); mpfr_init(r17376); mpfr_init(r17377); mpfr_init(r17378); mpfr_init(r17379); mpfr_init_set_str(r17380, "4", 10, MPFR_RNDN); mpfr_init(r17381); mpfr_init(r17382); mpfr_init(r17383); mpfr_init(r17384); mpfr_init(r17385); mpfr_init(r17386); mpfr_init(r17387); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r17374, r17372, r17373, MPFR_RNDN); mpfr_set_d(r17375, b, MPFR_RNDN); mpfr_neg(r17376, r17375, MPFR_RNDN); mpfr_sqr(r17377, r17375, MPFR_RNDN); mpfr_set_d(r17378, a, MPFR_RNDN); mpfr_set_d(r17379, c, MPFR_RNDN); ; mpfr_mul(r17381, r17379, r17380, MPFR_RNDN); mpfr_mul(r17382, r17378, r17381, MPFR_RNDN); mpfr_sub(r17383, r17377, r17382, MPFR_RNDN); mpfr_sqrt(r17384, r17383, MPFR_RNDN); mpfr_add(r17385, r17376, r17384, MPFR_RNDN); mpfr_div(r17386, r17385, r17378, MPFR_RNDN); mpfr_mul(r17387, r17374, r17386, MPFR_RNDN); return mpfr_get_d(r17387, MPFR_RNDN); }