#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r24379 = b; float r24380 = -r24379; float r24381 = r24379 * r24379; float r24382 = 4; float r24383 = a; float r24384 = r24382 * r24383; float r24385 = c; float r24386 = r24384 * r24385; float r24387 = r24381 - r24386; float r24388 = sqrt(r24387); float r24389 = r24380 + r24388; float r24390 = 2; float r24391 = r24390 * r24383; float r24392 = r24389 / r24391; return r24392; } double f_id(double a, double b, double c) { double r24393 = b; double r24394 = -r24393; double r24395 = r24393 * r24393; double r24396 = 4; double r24397 = a; double r24398 = r24396 * r24397; double r24399 = c; double r24400 = r24398 * r24399; double r24401 = r24395 - r24400; double r24402 = sqrt(r24401); double r24403 = r24394 + r24402; double r24404 = 2; double r24405 = r24404 * r24397; double r24406 = r24403 / r24405; return r24406; } double f_of(float a, float b, float c) { float r24407 = c; float r24408 = 2; float r24409 = r24407 / r24408; float r24410 = 4; float r24411 = b; float r24412 = -r24411; float r24413 = r24411 * r24411; float r24414 = a; float r24415 = r24410 * r24414; float r24416 = r24407 * r24415; float r24417 = r24413 - r24416; float r24418 = sqrt(r24417); float r24419 = r24412 - r24418; float r24420 = r24410 / r24419; float r24421 = r24409 * r24420; return r24421; } double f_od(double a, double b, double c) { double r24422 = c; double r24423 = 2; double r24424 = r24422 / r24423; double r24425 = 4; double r24426 = b; double r24427 = -r24426; double r24428 = r24426 * r24426; double r24429 = a; double r24430 = r24425 * r24429; double r24431 = r24422 * r24430; double r24432 = r24428 - r24431; double r24433 = sqrt(r24432); double r24434 = r24427 - r24433; double r24435 = r24425 / r24434; double r24436 = r24424 * r24435; return r24436; } 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 r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init_set_str(r24440, "4", 10, MPFR_RNDN); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init_set_str(r24448, "2", 10, MPFR_RNDN); mpfr_init(r24449); mpfr_init(r24450); } double f_im(double a, double b, double c) { mpfr_set_d(r24437, b, MPFR_RNDN); mpfr_neg(r24438, r24437, MPFR_RNDN); mpfr_mul(r24439, r24437, r24437, MPFR_RNDN); ; mpfr_set_d(r24441, a, MPFR_RNDN); mpfr_mul(r24442, r24440, r24441, MPFR_RNDN); mpfr_set_d(r24443, c, MPFR_RNDN); mpfr_mul(r24444, r24442, r24443, MPFR_RNDN); mpfr_sub(r24445, r24439, r24444, MPFR_RNDN); mpfr_sqrt(r24446, r24445, MPFR_RNDN); mpfr_add(r24447, r24438, r24446, MPFR_RNDN); ; mpfr_mul(r24449, r24448, r24441, MPFR_RNDN); mpfr_div(r24450, r24447, r24449, MPFR_RNDN); return mpfr_get_d(r24450, MPFR_RNDN); } static mpfr_t r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24451); mpfr_init_set_str(r24452, "2", 10, MPFR_RNDN); mpfr_init(r24453); mpfr_init_set_str(r24454, "4", 10, MPFR_RNDN); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); mpfr_init(r24464); mpfr_init(r24465); } double f_fm(double a, double b, double c) { mpfr_set_d(r24451, c, MPFR_RNDN); ; mpfr_div(r24453, r24451, r24452, MPFR_RNDN); ; mpfr_set_d(r24455, b, MPFR_RNDN); mpfr_neg(r24456, r24455, MPFR_RNDN); mpfr_mul(r24457, r24455, r24455, MPFR_RNDN); mpfr_set_d(r24458, a, MPFR_RNDN); mpfr_mul(r24459, r24454, r24458, MPFR_RNDN); mpfr_mul(r24460, r24451, r24459, MPFR_RNDN); mpfr_sub(r24461, r24457, r24460, MPFR_RNDN); mpfr_sqrt(r24462, r24461, MPFR_RNDN); mpfr_sub(r24463, r24456, r24462, MPFR_RNDN); mpfr_div(r24464, r24454, r24463, MPFR_RNDN); mpfr_mul(r24465, r24453, r24464, MPFR_RNDN); return mpfr_get_d(r24465, MPFR_RNDN); } static mpfr_t r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24466); mpfr_init_set_str(r24467, "2", 10, MPFR_RNDN); mpfr_init(r24468); mpfr_init_set_str(r24469, "4", 10, MPFR_RNDN); mpfr_init(r24470); mpfr_init(r24471); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); } double f_dm(double a, double b, double c) { mpfr_set_d(r24466, c, MPFR_RNDN); ; mpfr_div(r24468, r24466, r24467, MPFR_RNDN); ; mpfr_set_d(r24470, b, MPFR_RNDN); mpfr_neg(r24471, r24470, MPFR_RNDN); mpfr_mul(r24472, r24470, r24470, MPFR_RNDN); mpfr_set_d(r24473, a, MPFR_RNDN); mpfr_mul(r24474, r24469, r24473, MPFR_RNDN); mpfr_mul(r24475, r24466, r24474, MPFR_RNDN); mpfr_sub(r24476, r24472, r24475, MPFR_RNDN); mpfr_sqrt(r24477, r24476, MPFR_RNDN); mpfr_sub(r24478, r24471, r24477, MPFR_RNDN); mpfr_div(r24479, r24469, r24478, MPFR_RNDN); mpfr_mul(r24480, r24468, r24479, MPFR_RNDN); return mpfr_get_d(r24480, MPFR_RNDN); }