#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17333 = b; float r17334 = -r17333; float r17335 = r17333 * r17333; float r17336 = 4.0f; float r17337 = a; float r17338 = r17336 * r17337; float r17339 = c; float r17340 = r17338 * r17339; float r17341 = r17335 - r17340; float r17342 = sqrt(r17341); float r17343 = r17334 + r17342; float r17344 = 2.0f; float r17345 = r17344 * r17337; float r17346 = r17343 / r17345; return r17346; } double f_id(double a, double b, double c) { double r17347 = b; double r17348 = -r17347; double r17349 = r17347 * r17347; double r17350 = 4.0; double r17351 = a; double r17352 = r17350 * r17351; double r17353 = c; double r17354 = r17352 * r17353; double r17355 = r17349 - r17354; double r17356 = sqrt(r17355); double r17357 = r17348 + r17356; double r17358 = 2.0; double r17359 = r17358 * r17351; double r17360 = r17357 / r17359; return r17360; } double f_of(float a, float b, float c) { float r17361 = 4.0f; float r17362 = a; float r17363 = r17361 * r17362; float r17364 = 1.0f; float r17365 = r17363 / r17364; float r17366 = c; float r17367 = b; float r17368 = -r17367; float r17369 = r17367 * r17367; float r17370 = r17363 * r17366; float r17371 = r17369 - r17370; float r17372 = sqrt(r17371); float r17373 = r17368 - r17372; float r17374 = r17366 / r17373; float r17375 = r17365 * r17374; float r17376 = 2.0f; float r17377 = r17376 * r17362; float r17378 = r17375 / r17377; return r17378; } double f_od(double a, double b, double c) { double r17379 = 4.0; double r17380 = a; double r17381 = r17379 * r17380; double r17382 = 1.0; double r17383 = r17381 / r17382; double r17384 = c; double r17385 = b; double r17386 = -r17385; double r17387 = r17385 * r17385; double r17388 = r17381 * r17384; double r17389 = r17387 - r17388; double r17390 = sqrt(r17389); double r17391 = r17386 - r17390; double r17392 = r17384 / r17391; double r17393 = r17383 * r17392; double r17394 = 2.0; double r17395 = r17394 * r17380; double r17396 = r17393 / r17395; return r17396; } 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 r17397, r17398, r17399, r17400, r17401, r17402, r17403, r17404, r17405, r17406, r17407, r17408, r17409, r17410; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17397); mpfr_init(r17398); mpfr_init(r17399); mpfr_init_set_str(r17400, "4", 10, MPFR_RNDN); mpfr_init(r17401); mpfr_init(r17402); mpfr_init(r17403); mpfr_init(r17404); mpfr_init(r17405); mpfr_init(r17406); mpfr_init(r17407); mpfr_init_set_str(r17408, "2", 10, MPFR_RNDN); mpfr_init(r17409); mpfr_init(r17410); } double f_im(double a, double b, double c) { mpfr_set_d(r17397, b, MPFR_RNDN); mpfr_neg(r17398, r17397, MPFR_RNDN); mpfr_sqr(r17399, r17397, MPFR_RNDN); ; mpfr_set_d(r17401, a, MPFR_RNDN); mpfr_mul(r17402, r17400, r17401, MPFR_RNDN); mpfr_set_d(r17403, c, MPFR_RNDN); mpfr_mul(r17404, r17402, r17403, MPFR_RNDN); mpfr_sub(r17405, r17399, r17404, MPFR_RNDN); mpfr_sqrt(r17406, r17405, MPFR_RNDN); mpfr_add(r17407, r17398, r17406, MPFR_RNDN); ; mpfr_mul(r17409, r17408, r17401, MPFR_RNDN); mpfr_div(r17410, r17407, r17409, MPFR_RNDN); return mpfr_get_d(r17410, MPFR_RNDN); } static mpfr_t r17411, r17412, r17413, r17414, r17415, r17416, r17417, r17418, r17419, r17420, r17421, r17422, r17423, r17424, r17425, r17426, r17427, r17428; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17411, "4", 10, MPFR_RNDN); mpfr_init(r17412); mpfr_init(r17413); mpfr_init_set_str(r17414, "1", 10, MPFR_RNDN); mpfr_init(r17415); mpfr_init(r17416); mpfr_init(r17417); mpfr_init(r17418); mpfr_init(r17419); mpfr_init(r17420); mpfr_init(r17421); mpfr_init(r17422); mpfr_init(r17423); mpfr_init(r17424); mpfr_init(r17425); mpfr_init_set_str(r17426, "2", 10, MPFR_RNDN); mpfr_init(r17427); mpfr_init(r17428); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17412, a, MPFR_RNDN); mpfr_mul(r17413, r17411, r17412, MPFR_RNDN); ; mpfr_div(r17415, r17413, r17414, MPFR_RNDN); mpfr_set_d(r17416, c, MPFR_RNDN); mpfr_set_d(r17417, b, MPFR_RNDN); mpfr_neg(r17418, r17417, MPFR_RNDN); mpfr_sqr(r17419, r17417, MPFR_RNDN); mpfr_mul(r17420, r17413, r17416, MPFR_RNDN); mpfr_sub(r17421, r17419, r17420, MPFR_RNDN); mpfr_sqrt(r17422, r17421, MPFR_RNDN); mpfr_sub(r17423, r17418, r17422, MPFR_RNDN); mpfr_div(r17424, r17416, r17423, MPFR_RNDN); mpfr_mul(r17425, r17415, r17424, MPFR_RNDN); ; mpfr_mul(r17427, r17426, r17412, MPFR_RNDN); mpfr_div(r17428, r17425, r17427, MPFR_RNDN); return mpfr_get_d(r17428, MPFR_RNDN); } static mpfr_t r17429, r17430, r17431, r17432, r17433, r17434, r17435, r17436, r17437, r17438, r17439, r17440, r17441, r17442, r17443, r17444, r17445, r17446; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17429, "4", 10, MPFR_RNDN); mpfr_init(r17430); mpfr_init(r17431); mpfr_init_set_str(r17432, "1", 10, MPFR_RNDN); mpfr_init(r17433); mpfr_init(r17434); mpfr_init(r17435); mpfr_init(r17436); mpfr_init(r17437); mpfr_init(r17438); mpfr_init(r17439); mpfr_init(r17440); mpfr_init(r17441); mpfr_init(r17442); mpfr_init(r17443); mpfr_init_set_str(r17444, "2", 10, MPFR_RNDN); mpfr_init(r17445); mpfr_init(r17446); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17430, a, MPFR_RNDN); mpfr_mul(r17431, r17429, r17430, MPFR_RNDN); ; mpfr_div(r17433, r17431, r17432, MPFR_RNDN); mpfr_set_d(r17434, c, MPFR_RNDN); mpfr_set_d(r17435, b, MPFR_RNDN); mpfr_neg(r17436, r17435, MPFR_RNDN); mpfr_sqr(r17437, r17435, MPFR_RNDN); mpfr_mul(r17438, r17431, r17434, MPFR_RNDN); mpfr_sub(r17439, r17437, r17438, MPFR_RNDN); mpfr_sqrt(r17440, r17439, MPFR_RNDN); mpfr_sub(r17441, r17436, r17440, MPFR_RNDN); mpfr_div(r17442, r17434, r17441, MPFR_RNDN); mpfr_mul(r17443, r17433, r17442, MPFR_RNDN); ; mpfr_mul(r17445, r17444, r17430, MPFR_RNDN); mpfr_div(r17446, r17443, r17445, MPFR_RNDN); return mpfr_get_d(r17446, MPFR_RNDN); }