#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 = c; float r17365 = r17363 * r17364; float r17366 = b; float r17367 = -r17366; float r17368 = r17366 * r17366; float r17369 = r17368 - r17365; float r17370 = sqrt(r17369); float r17371 = r17367 - r17370; float r17372 = r17365 / r17371; float r17373 = 2.0f; float r17374 = r17373 * r17362; float r17375 = r17372 / r17374; return r17375; } double f_od(double a, double b, double c) { double r17376 = 4.0; double r17377 = a; double r17378 = r17376 * r17377; double r17379 = c; double r17380 = r17378 * r17379; double r17381 = b; double r17382 = -r17381; double r17383 = r17381 * r17381; double r17384 = r17383 - r17380; double r17385 = sqrt(r17384); double r17386 = r17382 - r17385; double r17387 = r17380 / r17386; double r17388 = 2.0; double r17389 = r17388 * r17377; double r17390 = r17387 / r17389; return r17390; } 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 r17391, r17392, r17393, r17394, r17395, r17396, r17397, r17398, r17399, r17400, r17401, r17402, r17403, r17404; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17391); mpfr_init(r17392); mpfr_init(r17393); mpfr_init_set_str(r17394, "4", 10, MPFR_RNDN); mpfr_init(r17395); mpfr_init(r17396); mpfr_init(r17397); mpfr_init(r17398); mpfr_init(r17399); mpfr_init(r17400); mpfr_init(r17401); mpfr_init_set_str(r17402, "2", 10, MPFR_RNDN); mpfr_init(r17403); mpfr_init(r17404); } double f_im(double a, double b, double c) { mpfr_set_d(r17391, b, MPFR_RNDN); mpfr_neg(r17392, r17391, MPFR_RNDN); mpfr_sqr(r17393, r17391, MPFR_RNDN); ; mpfr_set_d(r17395, a, MPFR_RNDN); mpfr_mul(r17396, r17394, r17395, MPFR_RNDN); mpfr_set_d(r17397, c, MPFR_RNDN); mpfr_mul(r17398, r17396, r17397, MPFR_RNDN); mpfr_sub(r17399, r17393, r17398, MPFR_RNDN); mpfr_sqrt(r17400, r17399, MPFR_RNDN); mpfr_add(r17401, r17392, r17400, MPFR_RNDN); ; mpfr_mul(r17403, r17402, r17395, MPFR_RNDN); mpfr_div(r17404, r17401, r17403, MPFR_RNDN); return mpfr_get_d(r17404, MPFR_RNDN); } static mpfr_t r17405, r17406, r17407, r17408, r17409, r17410, r17411, r17412, r17413, r17414, r17415, r17416, r17417, r17418, r17419; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17405, "4", 10, MPFR_RNDN); mpfr_init(r17406); mpfr_init(r17407); mpfr_init(r17408); mpfr_init(r17409); mpfr_init(r17410); mpfr_init(r17411); mpfr_init(r17412); mpfr_init(r17413); mpfr_init(r17414); mpfr_init(r17415); mpfr_init(r17416); mpfr_init_set_str(r17417, "2", 10, MPFR_RNDN); mpfr_init(r17418); mpfr_init(r17419); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17406, a, MPFR_RNDN); mpfr_mul(r17407, r17405, r17406, MPFR_RNDN); mpfr_set_d(r17408, c, MPFR_RNDN); mpfr_mul(r17409, r17407, r17408, MPFR_RNDN); mpfr_set_d(r17410, b, MPFR_RNDN); mpfr_neg(r17411, r17410, MPFR_RNDN); mpfr_sqr(r17412, r17410, MPFR_RNDN); mpfr_sub(r17413, r17412, r17409, MPFR_RNDN); mpfr_sqrt(r17414, r17413, MPFR_RNDN); mpfr_sub(r17415, r17411, r17414, MPFR_RNDN); mpfr_div(r17416, r17409, r17415, MPFR_RNDN); ; mpfr_mul(r17418, r17417, r17406, MPFR_RNDN); mpfr_div(r17419, r17416, r17418, MPFR_RNDN); return mpfr_get_d(r17419, MPFR_RNDN); } static mpfr_t r17420, r17421, r17422, r17423, r17424, r17425, r17426, r17427, r17428, r17429, r17430, r17431, r17432, r17433, r17434; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17420, "4", 10, MPFR_RNDN); mpfr_init(r17421); mpfr_init(r17422); mpfr_init(r17423); mpfr_init(r17424); mpfr_init(r17425); mpfr_init(r17426); mpfr_init(r17427); mpfr_init(r17428); mpfr_init(r17429); mpfr_init(r17430); mpfr_init(r17431); mpfr_init_set_str(r17432, "2", 10, MPFR_RNDN); mpfr_init(r17433); mpfr_init(r17434); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17421, a, MPFR_RNDN); mpfr_mul(r17422, r17420, r17421, MPFR_RNDN); mpfr_set_d(r17423, c, MPFR_RNDN); mpfr_mul(r17424, r17422, r17423, MPFR_RNDN); mpfr_set_d(r17425, b, MPFR_RNDN); mpfr_neg(r17426, r17425, MPFR_RNDN); mpfr_sqr(r17427, r17425, MPFR_RNDN); mpfr_sub(r17428, r17427, r17424, MPFR_RNDN); mpfr_sqrt(r17429, r17428, MPFR_RNDN); mpfr_sub(r17430, r17426, r17429, MPFR_RNDN); mpfr_div(r17431, r17424, r17430, MPFR_RNDN); ; mpfr_mul(r17433, r17432, r17421, MPFR_RNDN); mpfr_div(r17434, r17431, r17433, MPFR_RNDN); return mpfr_get_d(r17434, MPFR_RNDN); }