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