#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r24339 = b; float r24340 = -r24339; float r24341 = r24339 * r24339; float r24342 = 4; float r24343 = a; float r24344 = r24342 * r24343; float r24345 = c; float r24346 = r24344 * r24345; float r24347 = r24341 - r24346; float r24348 = sqrt(r24347); float r24349 = r24340 + r24348; float r24350 = 2; float r24351 = r24350 * r24343; float r24352 = r24349 / r24351; return r24352; } double f_id(double a, double b, double c) { double r24353 = b; double r24354 = -r24353; double r24355 = r24353 * r24353; double r24356 = 4; double r24357 = a; double r24358 = r24356 * r24357; double r24359 = c; double r24360 = r24358 * r24359; double r24361 = r24355 - r24360; double r24362 = sqrt(r24361); double r24363 = r24354 + r24362; double r24364 = 2; double r24365 = r24364 * r24357; double r24366 = r24363 / r24365; return r24366; } double f_of(float a, float b, float c) { float r24367 = c; float r24368 = 2; float r24369 = r24367 / r24368; float r24370 = 4; float r24371 = b; float r24372 = -r24371; float r24373 = r24371 * r24371; float r24374 = a; float r24375 = r24370 * r24374; float r24376 = r24367 * r24375; float r24377 = r24373 - r24376; float r24378 = sqrt(r24377); float r24379 = r24372 - r24378; float r24380 = r24370 / r24379; float r24381 = r24369 * r24380; return r24381; } double f_od(double a, double b, double c) { double r24382 = c; double r24383 = 2; double r24384 = r24382 / r24383; double r24385 = 4; double r24386 = b; double r24387 = -r24386; double r24388 = r24386 * r24386; double r24389 = a; double r24390 = r24385 * r24389; double r24391 = r24382 * r24390; double r24392 = r24388 - r24391; double r24393 = sqrt(r24392); double r24394 = r24387 - r24393; double r24395 = r24385 / r24394; double r24396 = r24384 * r24395; return r24396; } 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 r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24397); mpfr_init(r24398); mpfr_init(r24399); mpfr_init_set_str(r24400, "4", 10, MPFR_RNDN); mpfr_init(r24401); mpfr_init(r24402); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init(r24406); mpfr_init(r24407); mpfr_init_set_str(r24408, "2", 10, MPFR_RNDN); mpfr_init(r24409); mpfr_init(r24410); } double f_im(double a, double b, double c) { mpfr_set_d(r24397, b, MPFR_RNDN); mpfr_neg(r24398, r24397, MPFR_RNDN); mpfr_mul(r24399, r24397, r24397, MPFR_RNDN); ; mpfr_set_d(r24401, a, MPFR_RNDN); mpfr_mul(r24402, r24400, r24401, MPFR_RNDN); mpfr_set_d(r24403, c, MPFR_RNDN); mpfr_mul(r24404, r24402, r24403, MPFR_RNDN); mpfr_sub(r24405, r24399, r24404, MPFR_RNDN); mpfr_sqrt(r24406, r24405, MPFR_RNDN); mpfr_add(r24407, r24398, r24406, MPFR_RNDN); ; mpfr_mul(r24409, r24408, r24401, MPFR_RNDN); mpfr_div(r24410, r24407, r24409, MPFR_RNDN); return mpfr_get_d(r24410, MPFR_RNDN); } static mpfr_t r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24411); mpfr_init_set_str(r24412, "2", 10, MPFR_RNDN); mpfr_init(r24413); mpfr_init_set_str(r24414, "4", 10, MPFR_RNDN); mpfr_init(r24415); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init(r24421); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); } double f_fm(double a, double b, double c) { mpfr_set_d(r24411, c, MPFR_RNDN); ; mpfr_div(r24413, r24411, r24412, MPFR_RNDN); ; mpfr_set_d(r24415, b, MPFR_RNDN); mpfr_neg(r24416, r24415, MPFR_RNDN); mpfr_mul(r24417, r24415, r24415, MPFR_RNDN); mpfr_set_d(r24418, a, MPFR_RNDN); mpfr_mul(r24419, r24414, r24418, MPFR_RNDN); mpfr_mul(r24420, r24411, r24419, MPFR_RNDN); mpfr_sub(r24421, r24417, r24420, MPFR_RNDN); mpfr_sqrt(r24422, r24421, MPFR_RNDN); mpfr_sub(r24423, r24416, r24422, MPFR_RNDN); mpfr_div(r24424, r24414, r24423, MPFR_RNDN); mpfr_mul(r24425, r24413, r24424, MPFR_RNDN); return mpfr_get_d(r24425, MPFR_RNDN); } static mpfr_t r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24426); mpfr_init_set_str(r24427, "2", 10, MPFR_RNDN); mpfr_init(r24428); mpfr_init_set_str(r24429, "4", 10, MPFR_RNDN); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); } double f_dm(double a, double b, double c) { mpfr_set_d(r24426, c, MPFR_RNDN); ; mpfr_div(r24428, r24426, r24427, MPFR_RNDN); ; mpfr_set_d(r24430, b, MPFR_RNDN); mpfr_neg(r24431, r24430, MPFR_RNDN); mpfr_mul(r24432, r24430, r24430, MPFR_RNDN); mpfr_set_d(r24433, a, MPFR_RNDN); mpfr_mul(r24434, r24429, r24433, MPFR_RNDN); mpfr_mul(r24435, r24426, r24434, MPFR_RNDN); mpfr_sub(r24436, r24432, r24435, MPFR_RNDN); mpfr_sqrt(r24437, r24436, MPFR_RNDN); mpfr_sub(r24438, r24431, r24437, MPFR_RNDN); mpfr_div(r24439, r24429, r24438, MPFR_RNDN); mpfr_mul(r24440, r24428, r24439, MPFR_RNDN); return mpfr_get_d(r24440, MPFR_RNDN); }