#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r21466 = b; float r21467 = -r21466; float r21468 = r21466 * r21466; float r21469 = 4; float r21470 = a; float r21471 = r21469 * r21470; float r21472 = c; float r21473 = r21471 * r21472; float r21474 = r21468 - r21473; float r21475 = sqrt(r21474); float r21476 = r21467 + r21475; float r21477 = 2; float r21478 = r21477 * r21470; float r21479 = r21476 / r21478; return r21479; } double f_id(double a, double b, double c) { double r21480 = b; double r21481 = -r21480; double r21482 = r21480 * r21480; double r21483 = 4; double r21484 = a; double r21485 = r21483 * r21484; double r21486 = c; double r21487 = r21485 * r21486; double r21488 = r21482 - r21487; double r21489 = sqrt(r21488); double r21490 = r21481 + r21489; double r21491 = 2; double r21492 = r21491 * r21484; double r21493 = r21490 / r21492; return r21493; } double f_of(float a, float b, float c) { float r21494 = 1; float r21495 = b; float r21496 = -r21495; float r21497 = a; float r21498 = 4; float r21499 = r21497 * r21498; float r21500 = c; float r21501 = -r21500; float r21502 = r21495 * r21495; float r21503 = fma(r21499, r21501, r21502); float r21504 = sqrt(r21503); float r21505 = r21496 - r21504; float r21506 = 2; float r21507 = r21506 / r21500; float r21508 = r21494 / r21498; float r21509 = r21507 * r21508; float r21510 = r21505 * r21509; float r21511 = r21494 / r21510; return r21511; } double f_od(double a, double b, double c) { double r21512 = 1; double r21513 = b; double r21514 = -r21513; double r21515 = a; double r21516 = 4; double r21517 = r21515 * r21516; double r21518 = c; double r21519 = -r21518; double r21520 = r21513 * r21513; double r21521 = fma(r21517, r21519, r21520); double r21522 = sqrt(r21521); double r21523 = r21514 - r21522; double r21524 = 2; double r21525 = r21524 / r21518; double r21526 = r21512 / r21516; double r21527 = r21525 * r21526; double r21528 = r21523 * r21527; double r21529 = r21512 / r21528; return r21529; } 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 r21530, r21531, r21532, r21533, r21534, r21535, r21536, r21537, r21538, r21539, r21540, r21541, r21542, r21543; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21530); mpfr_init(r21531); mpfr_init(r21532); mpfr_init_set_str(r21533, "4", 10, MPFR_RNDN); mpfr_init(r21534); mpfr_init(r21535); mpfr_init(r21536); mpfr_init(r21537); mpfr_init(r21538); mpfr_init(r21539); mpfr_init(r21540); mpfr_init_set_str(r21541, "2", 10, MPFR_RNDN); mpfr_init(r21542); mpfr_init(r21543); } double f_im(double a, double b, double c) { mpfr_set_d(r21530, b, MPFR_RNDN); mpfr_neg(r21531, r21530, MPFR_RNDN); mpfr_mul(r21532, r21530, r21530, MPFR_RNDN); ; mpfr_set_d(r21534, a, MPFR_RNDN); mpfr_mul(r21535, r21533, r21534, MPFR_RNDN); mpfr_set_d(r21536, c, MPFR_RNDN); mpfr_mul(r21537, r21535, r21536, MPFR_RNDN); mpfr_sub(r21538, r21532, r21537, MPFR_RNDN); mpfr_sqrt(r21539, r21538, MPFR_RNDN); mpfr_add(r21540, r21531, r21539, MPFR_RNDN); ; mpfr_mul(r21542, r21541, r21534, MPFR_RNDN); mpfr_div(r21543, r21540, r21542, MPFR_RNDN); return mpfr_get_d(r21543, MPFR_RNDN); } static mpfr_t r21544, r21545, r21546, r21547, r21548, r21549, r21550, r21551, r21552, r21553, r21554, r21555, r21556, r21557, r21558, r21559, r21560, r21561; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21544, "1", 10, MPFR_RNDN); mpfr_init(r21545); mpfr_init(r21546); mpfr_init(r21547); mpfr_init_set_str(r21548, "4", 10, MPFR_RNDN); mpfr_init(r21549); mpfr_init(r21550); mpfr_init(r21551); mpfr_init(r21552); mpfr_init(r21553); mpfr_init(r21554); mpfr_init(r21555); mpfr_init_set_str(r21556, "2", 10, MPFR_RNDN); mpfr_init(r21557); mpfr_init(r21558); mpfr_init(r21559); mpfr_init(r21560); mpfr_init(r21561); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r21545, b, MPFR_RNDN); mpfr_neg(r21546, r21545, MPFR_RNDN); mpfr_set_d(r21547, a, MPFR_RNDN); ; mpfr_mul(r21549, r21547, r21548, MPFR_RNDN); mpfr_set_d(r21550, c, MPFR_RNDN); mpfr_neg(r21551, r21550, MPFR_RNDN); mpfr_mul(r21552, r21545, r21545, MPFR_RNDN); mpfr_fma(r21553, r21549, r21551, r21552, MPFR_RNDN); mpfr_sqrt(r21554, r21553, MPFR_RNDN); mpfr_sub(r21555, r21546, r21554, MPFR_RNDN); ; mpfr_div(r21557, r21556, r21550, MPFR_RNDN); mpfr_div(r21558, r21544, r21548, MPFR_RNDN); mpfr_mul(r21559, r21557, r21558, MPFR_RNDN); mpfr_mul(r21560, r21555, r21559, MPFR_RNDN); mpfr_div(r21561, r21544, r21560, MPFR_RNDN); return mpfr_get_d(r21561, MPFR_RNDN); } static mpfr_t r21562, r21563, r21564, r21565, r21566, r21567, r21568, r21569, r21570, r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21562, "1", 10, MPFR_RNDN); mpfr_init(r21563); mpfr_init(r21564); mpfr_init(r21565); mpfr_init_set_str(r21566, "4", 10, MPFR_RNDN); mpfr_init(r21567); mpfr_init(r21568); mpfr_init(r21569); mpfr_init(r21570); mpfr_init(r21571); mpfr_init(r21572); mpfr_init(r21573); mpfr_init_set_str(r21574, "2", 10, MPFR_RNDN); mpfr_init(r21575); mpfr_init(r21576); mpfr_init(r21577); mpfr_init(r21578); mpfr_init(r21579); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r21563, b, MPFR_RNDN); mpfr_neg(r21564, r21563, MPFR_RNDN); mpfr_set_d(r21565, a, MPFR_RNDN); ; mpfr_mul(r21567, r21565, r21566, MPFR_RNDN); mpfr_set_d(r21568, c, MPFR_RNDN); mpfr_neg(r21569, r21568, MPFR_RNDN); mpfr_mul(r21570, r21563, r21563, MPFR_RNDN); mpfr_fma(r21571, r21567, r21569, r21570, MPFR_RNDN); mpfr_sqrt(r21572, r21571, MPFR_RNDN); mpfr_sub(r21573, r21564, r21572, MPFR_RNDN); ; mpfr_div(r21575, r21574, r21568, MPFR_RNDN); mpfr_div(r21576, r21562, r21566, MPFR_RNDN); mpfr_mul(r21577, r21575, r21576, MPFR_RNDN); mpfr_mul(r21578, r21573, r21577, MPFR_RNDN); mpfr_div(r21579, r21562, r21578, MPFR_RNDN); return mpfr_get_d(r21579, MPFR_RNDN); }