#include #include #include #include #include char *name = "Quadratic roots, wide range"; double f_if(float a, float b, float c) { float r24360 = b; float r24361 = -r24360; float r24362 = r24360 * r24360; float r24363 = 4; float r24364 = a; float r24365 = r24363 * r24364; float r24366 = c; float r24367 = r24365 * r24366; float r24368 = r24362 - r24367; float r24369 = sqrt(r24368); float r24370 = r24361 + r24369; float r24371 = 2; float r24372 = r24371 * r24364; float r24373 = r24370 / r24372; return r24373; } double f_id(double a, double b, double c) { double r24374 = b; double r24375 = -r24374; double r24376 = r24374 * r24374; double r24377 = 4; double r24378 = a; double r24379 = r24377 * r24378; double r24380 = c; double r24381 = r24379 * r24380; double r24382 = r24376 - r24381; double r24383 = sqrt(r24382); double r24384 = r24375 + r24383; double r24385 = 2; double r24386 = r24385 * r24378; double r24387 = r24384 / r24386; return r24387; } double f_of(float a, float b, float c) { float r24388 = c; float r24389 = 2; float r24390 = r24388 / r24389; float r24391 = 4; float r24392 = b; float r24393 = -r24392; float r24394 = -r24391; float r24395 = a; float r24396 = r24388 * r24395; float r24397 = r24392 * r24392; float r24398 = fma(r24394, r24396, r24397); float r24399 = sqrt(r24398); float r24400 = r24393 - r24399; float r24401 = r24391 / r24400; float r24402 = r24390 * r24401; return r24402; } double f_od(double a, double b, double c) { double r24403 = c; double r24404 = 2; double r24405 = r24403 / r24404; double r24406 = 4; double r24407 = b; double r24408 = -r24407; double r24409 = -r24406; double r24410 = a; double r24411 = r24403 * r24410; double r24412 = r24407 * r24407; double r24413 = fma(r24409, r24411, r24412); double r24414 = sqrt(r24413); double r24415 = r24408 - r24414; double r24416 = r24406 / r24415; double r24417 = r24405 * r24416; return r24417; } 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 r24418, r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init_set_str(r24421, "4", 10, MPFR_RNDN); mpfr_init(r24422); mpfr_init(r24423); mpfr_init(r24424); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init_set_str(r24429, "2", 10, MPFR_RNDN); mpfr_init(r24430); mpfr_init(r24431); } double f_im(double a, double b, double c) { mpfr_set_d(r24418, b, MPFR_RNDN); mpfr_neg(r24419, r24418, MPFR_RNDN); mpfr_mul(r24420, r24418, r24418, MPFR_RNDN); ; mpfr_set_d(r24422, a, MPFR_RNDN); mpfr_mul(r24423, r24421, r24422, MPFR_RNDN); mpfr_set_d(r24424, c, MPFR_RNDN); mpfr_mul(r24425, r24423, r24424, MPFR_RNDN); mpfr_sub(r24426, r24420, r24425, MPFR_RNDN); mpfr_sqrt(r24427, r24426, MPFR_RNDN); mpfr_add(r24428, r24419, r24427, MPFR_RNDN); ; mpfr_mul(r24430, r24429, r24422, MPFR_RNDN); mpfr_div(r24431, r24428, r24430, MPFR_RNDN); return mpfr_get_d(r24431, MPFR_RNDN); } static mpfr_t r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r24432); mpfr_init_set_str(r24433, "2", 10, MPFR_RNDN); mpfr_init(r24434); mpfr_init_set_str(r24435, "4", 10, MPFR_RNDN); mpfr_init(r24436); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); } double f_fm(double a, double b, double c) { mpfr_set_d(r24432, c, MPFR_RNDN); ; mpfr_div(r24434, r24432, r24433, MPFR_RNDN); ; mpfr_set_d(r24436, b, MPFR_RNDN); mpfr_neg(r24437, r24436, MPFR_RNDN); mpfr_neg(r24438, r24435, MPFR_RNDN); mpfr_set_d(r24439, a, MPFR_RNDN); mpfr_mul(r24440, r24432, r24439, MPFR_RNDN); mpfr_mul(r24441, r24436, r24436, MPFR_RNDN); mpfr_fma(r24442, r24438, r24440, r24441, MPFR_RNDN); mpfr_sqrt(r24443, r24442, MPFR_RNDN); mpfr_sub(r24444, r24437, r24443, MPFR_RNDN); mpfr_div(r24445, r24435, r24444, MPFR_RNDN); mpfr_mul(r24446, r24434, r24445, MPFR_RNDN); return mpfr_get_d(r24446, MPFR_RNDN); } static mpfr_t r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r24447); mpfr_init_set_str(r24448, "2", 10, MPFR_RNDN); mpfr_init(r24449); mpfr_init_set_str(r24450, "4", 10, MPFR_RNDN); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); } double f_dm(double a, double b, double c) { mpfr_set_d(r24447, c, MPFR_RNDN); ; mpfr_div(r24449, r24447, r24448, MPFR_RNDN); ; mpfr_set_d(r24451, b, MPFR_RNDN); mpfr_neg(r24452, r24451, MPFR_RNDN); mpfr_neg(r24453, r24450, MPFR_RNDN); mpfr_set_d(r24454, a, MPFR_RNDN); mpfr_mul(r24455, r24447, r24454, MPFR_RNDN); mpfr_mul(r24456, r24451, r24451, MPFR_RNDN); mpfr_fma(r24457, r24453, r24455, r24456, MPFR_RNDN); mpfr_sqrt(r24458, r24457, MPFR_RNDN); mpfr_sub(r24459, r24452, r24458, MPFR_RNDN); mpfr_div(r24460, r24450, r24459, MPFR_RNDN); mpfr_mul(r24461, r24449, r24460, MPFR_RNDN); return mpfr_get_d(r24461, MPFR_RNDN); }