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