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