#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r24522 = b; float r24523 = -r24522; float r24524 = r24522 * r24522; float r24525 = 4; float r24526 = a; float r24527 = r24525 * r24526; float r24528 = c; float r24529 = r24527 * r24528; float r24530 = r24524 - r24529; float r24531 = sqrt(r24530); float r24532 = r24523 + r24531; float r24533 = 2; float r24534 = r24533 * r24526; float r24535 = r24532 / r24534; return r24535; } double f_id(double a, double b, double c) { double r24536 = b; double r24537 = -r24536; double r24538 = r24536 * r24536; double r24539 = 4; double r24540 = a; double r24541 = r24539 * r24540; double r24542 = c; double r24543 = r24541 * r24542; double r24544 = r24538 - r24543; double r24545 = sqrt(r24544); double r24546 = r24537 + r24545; double r24547 = 2; double r24548 = r24547 * r24540; double r24549 = r24546 / r24548; return r24549; } double f_of(float a, float b, float c) { float r24550 = c; float r24551 = a; float r24552 = r24550 * r24551; float r24553 = 4; float r24554 = -r24553; float r24555 = r24552 * r24554; float r24556 = r24553 * r24551; float r24557 = -r24550; float r24558 = b; float r24559 = r24558 * r24558; float r24560 = fma(r24556, r24557, r24559); float r24561 = sqrt(r24560); float r24562 = 3; float r24563 = pow(r24561, r24562); float r24564 = pow(r24558, r24562); float r24565 = r24563 + r24564; float r24566 = r24555 / r24565; float r24567 = r24561 * r24561; float r24568 = r24561 * r24558; float r24569 = r24559 - r24568; float r24570 = r24567 + r24569; float r24571 = r24566 * r24570; float r24572 = 2; float r24573 = r24572 * r24551; float r24574 = r24571 / r24573; return r24574; } double f_od(double a, double b, double c) { double r24575 = c; double r24576 = a; double r24577 = r24575 * r24576; double r24578 = 4; double r24579 = -r24578; double r24580 = r24577 * r24579; double r24581 = r24578 * r24576; double r24582 = -r24575; double r24583 = b; double r24584 = r24583 * r24583; double r24585 = fma(r24581, r24582, r24584); double r24586 = sqrt(r24585); double r24587 = 3; double r24588 = pow(r24586, r24587); double r24589 = pow(r24583, r24587); double r24590 = r24588 + r24589; double r24591 = r24580 / r24590; double r24592 = r24586 * r24586; double r24593 = r24586 * r24583; double r24594 = r24584 - r24593; double r24595 = r24592 + r24594; double r24596 = r24591 * r24595; double r24597 = 2; double r24598 = r24597 * r24576; double r24599 = r24596 / r24598; return r24599; } 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 r24600, r24601, r24602, r24603, r24604, r24605, r24606, r24607, r24608, r24609, r24610, r24611, r24612, r24613; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24600); mpfr_init(r24601); mpfr_init(r24602); mpfr_init_set_str(r24603, "4", 10, MPFR_RNDN); mpfr_init(r24604); mpfr_init(r24605); mpfr_init(r24606); mpfr_init(r24607); mpfr_init(r24608); mpfr_init(r24609); mpfr_init(r24610); mpfr_init_set_str(r24611, "2", 10, MPFR_RNDN); mpfr_init(r24612); mpfr_init(r24613); } double f_im(double a, double b, double c) { mpfr_set_d(r24600, b, MPFR_RNDN); mpfr_neg(r24601, r24600, MPFR_RNDN); mpfr_mul(r24602, r24600, r24600, MPFR_RNDN); ; mpfr_set_d(r24604, a, MPFR_RNDN); mpfr_mul(r24605, r24603, r24604, MPFR_RNDN); mpfr_set_d(r24606, c, MPFR_RNDN); mpfr_mul(r24607, r24605, r24606, MPFR_RNDN); mpfr_sub(r24608, r24602, r24607, MPFR_RNDN); mpfr_sqrt(r24609, r24608, MPFR_RNDN); mpfr_add(r24610, r24601, r24609, MPFR_RNDN); ; mpfr_mul(r24612, r24611, r24604, MPFR_RNDN); mpfr_div(r24613, r24610, r24612, MPFR_RNDN); return mpfr_get_d(r24613, MPFR_RNDN); } static mpfr_t r24614, r24615, r24616, r24617, r24618, r24619, r24620, r24621, r24622, r24623, r24624, r24625, r24626, r24627, r24628, r24629, r24630, r24631, r24632, r24633, r24634, r24635, r24636, r24637, r24638; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24614); mpfr_init(r24615); mpfr_init(r24616); mpfr_init_set_str(r24617, "4", 10, MPFR_RNDN); mpfr_init(r24618); mpfr_init(r24619); mpfr_init(r24620); mpfr_init(r24621); mpfr_init(r24622); mpfr_init(r24623); mpfr_init(r24624); mpfr_init(r24625); mpfr_init_set_str(r24626, "3", 10, MPFR_RNDN); mpfr_init(r24627); mpfr_init(r24628); mpfr_init(r24629); mpfr_init(r24630); mpfr_init(r24631); mpfr_init(r24632); mpfr_init(r24633); mpfr_init(r24634); mpfr_init(r24635); mpfr_init_set_str(r24636, "2", 10, MPFR_RNDN); mpfr_init(r24637); mpfr_init(r24638); } double f_fm(double a, double b, double c) { mpfr_set_d(r24614, c, MPFR_RNDN); mpfr_set_d(r24615, a, MPFR_RNDN); mpfr_mul(r24616, r24614, r24615, MPFR_RNDN); ; mpfr_neg(r24618, r24617, MPFR_RNDN); mpfr_mul(r24619, r24616, r24618, MPFR_RNDN); mpfr_mul(r24620, r24617, r24615, MPFR_RNDN); mpfr_neg(r24621, r24614, MPFR_RNDN); mpfr_set_d(r24622, b, MPFR_RNDN); mpfr_mul(r24623, r24622, r24622, MPFR_RNDN); mpfr_fma(r24624, r24620, r24621, r24623, MPFR_RNDN); mpfr_sqrt(r24625, r24624, MPFR_RNDN); ; mpfr_pow(r24627, r24625, r24626, MPFR_RNDN); mpfr_pow(r24628, r24622, r24626, MPFR_RNDN); mpfr_add(r24629, r24627, r24628, MPFR_RNDN); mpfr_div(r24630, r24619, r24629, MPFR_RNDN); mpfr_mul(r24631, r24625, r24625, MPFR_RNDN); mpfr_mul(r24632, r24625, r24622, MPFR_RNDN); mpfr_sub(r24633, r24623, r24632, MPFR_RNDN); mpfr_add(r24634, r24631, r24633, MPFR_RNDN); mpfr_mul(r24635, r24630, r24634, MPFR_RNDN); ; mpfr_mul(r24637, r24636, r24615, MPFR_RNDN); mpfr_div(r24638, r24635, r24637, MPFR_RNDN); return mpfr_get_d(r24638, MPFR_RNDN); } static mpfr_t r24639, r24640, r24641, r24642, r24643, r24644, r24645, r24646, r24647, r24648, r24649, r24650, r24651, r24652, r24653, r24654, r24655, r24656, r24657, r24658, r24659, r24660, r24661, r24662, r24663; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24639); mpfr_init(r24640); mpfr_init(r24641); mpfr_init_set_str(r24642, "4", 10, MPFR_RNDN); mpfr_init(r24643); mpfr_init(r24644); mpfr_init(r24645); mpfr_init(r24646); mpfr_init(r24647); mpfr_init(r24648); mpfr_init(r24649); mpfr_init(r24650); mpfr_init_set_str(r24651, "3", 10, MPFR_RNDN); mpfr_init(r24652); mpfr_init(r24653); mpfr_init(r24654); mpfr_init(r24655); mpfr_init(r24656); mpfr_init(r24657); mpfr_init(r24658); mpfr_init(r24659); mpfr_init(r24660); mpfr_init_set_str(r24661, "2", 10, MPFR_RNDN); mpfr_init(r24662); mpfr_init(r24663); } double f_dm(double a, double b, double c) { mpfr_set_d(r24639, c, MPFR_RNDN); mpfr_set_d(r24640, a, MPFR_RNDN); mpfr_mul(r24641, r24639, r24640, MPFR_RNDN); ; mpfr_neg(r24643, r24642, MPFR_RNDN); mpfr_mul(r24644, r24641, r24643, MPFR_RNDN); mpfr_mul(r24645, r24642, r24640, MPFR_RNDN); mpfr_neg(r24646, r24639, MPFR_RNDN); mpfr_set_d(r24647, b, MPFR_RNDN); mpfr_mul(r24648, r24647, r24647, MPFR_RNDN); mpfr_fma(r24649, r24645, r24646, r24648, MPFR_RNDN); mpfr_sqrt(r24650, r24649, MPFR_RNDN); ; mpfr_pow(r24652, r24650, r24651, MPFR_RNDN); mpfr_pow(r24653, r24647, r24651, MPFR_RNDN); mpfr_add(r24654, r24652, r24653, MPFR_RNDN); mpfr_div(r24655, r24644, r24654, MPFR_RNDN); mpfr_mul(r24656, r24650, r24650, MPFR_RNDN); mpfr_mul(r24657, r24650, r24647, MPFR_RNDN); mpfr_sub(r24658, r24648, r24657, MPFR_RNDN); mpfr_add(r24659, r24656, r24658, MPFR_RNDN); mpfr_mul(r24660, r24655, r24659, MPFR_RNDN); ; mpfr_mul(r24662, r24661, r24640, MPFR_RNDN); mpfr_div(r24663, r24660, r24662, MPFR_RNDN); return mpfr_get_d(r24663, MPFR_RNDN); }