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