#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r25524 = b; float r25525 = -r25524; float r25526 = r25524 * r25524; float r25527 = 4; float r25528 = a; float r25529 = r25527 * r25528; float r25530 = c; float r25531 = r25529 * r25530; float r25532 = r25526 - r25531; float r25533 = sqrt(r25532); float r25534 = r25525 + r25533; float r25535 = 2; float r25536 = r25535 * r25528; float r25537 = r25534 / r25536; return r25537; } double f_id(double a, double b, double c) { double r25538 = b; double r25539 = -r25538; double r25540 = r25538 * r25538; double r25541 = 4; double r25542 = a; double r25543 = r25541 * r25542; double r25544 = c; double r25545 = r25543 * r25544; double r25546 = r25540 - r25545; double r25547 = sqrt(r25546); double r25548 = r25539 + r25547; double r25549 = 2; double r25550 = r25549 * r25542; double r25551 = r25548 / r25550; return r25551; } double f_of(float a, float b, float c) { float r25552 = b; float r25553 = -2.364349591187642e+101; bool r25554 = r25552 <= r25553; float r25555 = c; float r25556 = r25555 / r25552; float r25557 = r25552 + r25552; float r25558 = 2; float r25559 = a; float r25560 = r25558 * r25559; float r25561 = r25557 / r25560; float r25562 = r25556 - r25561; float r25563 = 1.125868897500658e-82; bool r25564 = r25552 <= r25563; float r25565 = -r25552; float r25566 = r25552 * r25552; float r25567 = 4; float r25568 = r25567 * r25559; float r25569 = r25555 * r25568; float r25570 = r25566 - r25569; float r25571 = sqrt(r25570); float r25572 = r25565 + r25571; float r25573 = r25572 / r25560; float r25574 = -r25555; float r25575 = r25574 / r25552; float r25576 = r25564 ? r25573 : r25575; float r25577 = r25554 ? r25562 : r25576; return r25577; } double f_od(double a, double b, double c) { double r25578 = b; double r25579 = -2.364349591187642e+101; bool r25580 = r25578 <= r25579; double r25581 = c; double r25582 = r25581 / r25578; double r25583 = r25578 + r25578; double r25584 = 2; double r25585 = a; double r25586 = r25584 * r25585; double r25587 = r25583 / r25586; double r25588 = r25582 - r25587; double r25589 = 1.125868897500658e-82; bool r25590 = r25578 <= r25589; double r25591 = -r25578; double r25592 = r25578 * r25578; double r25593 = 4; double r25594 = r25593 * r25585; double r25595 = r25581 * r25594; double r25596 = r25592 - r25595; double r25597 = sqrt(r25596); double r25598 = r25591 + r25597; double r25599 = r25598 / r25586; double r25600 = -r25581; double r25601 = r25600 / r25578; double r25602 = r25590 ? r25599 : r25601; double r25603 = r25580 ? r25588 : r25602; return r25603; } 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 r25604, r25605, r25606, r25607, r25608, r25609, r25610, r25611, r25612, r25613, r25614, r25615, r25616, r25617; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init(r25604); mpfr_init(r25605); mpfr_init(r25606); mpfr_init_set_str(r25607, "4", 10, MPFR_RNDN); mpfr_init(r25608); mpfr_init(r25609); mpfr_init(r25610); mpfr_init(r25611); mpfr_init(r25612); mpfr_init(r25613); mpfr_init(r25614); mpfr_init_set_str(r25615, "2", 10, MPFR_RNDN); mpfr_init(r25616); mpfr_init(r25617); } double f_im(double a, double b, double c) { mpfr_set_d(r25604, b, MPFR_RNDN); mpfr_neg(r25605, r25604, MPFR_RNDN); mpfr_mul(r25606, r25604, r25604, MPFR_RNDN); ; mpfr_set_d(r25608, a, MPFR_RNDN); mpfr_mul(r25609, r25607, r25608, MPFR_RNDN); mpfr_set_d(r25610, c, MPFR_RNDN); mpfr_mul(r25611, r25609, r25610, MPFR_RNDN); mpfr_sub(r25612, r25606, r25611, MPFR_RNDN); mpfr_sqrt(r25613, r25612, MPFR_RNDN); mpfr_add(r25614, r25605, r25613, MPFR_RNDN); ; mpfr_mul(r25616, r25615, r25608, MPFR_RNDN); mpfr_div(r25617, r25614, r25616, MPFR_RNDN); return mpfr_get_d(r25617, MPFR_RNDN); } static mpfr_t r25618, r25619, r25620, r25621, r25622, r25623, r25624, r25625, r25626, r25627, r25628, r25629, r25630, r25631, r25632, r25633, r25634, r25635, r25636, r25637, r25638, r25639, r25640, r25641, r25642, r25643; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init(r25618); mpfr_init_set_str(r25619, "-2.364349591187642e+101", 10, MPFR_RNDN); mpfr_init(r25620); mpfr_init(r25621); mpfr_init(r25622); mpfr_init(r25623); mpfr_init_set_str(r25624, "2", 10, MPFR_RNDN); mpfr_init(r25625); mpfr_init(r25626); mpfr_init(r25627); mpfr_init(r25628); mpfr_init_set_str(r25629, "1.125868897500658e-82", 10, MPFR_RNDN); mpfr_init(r25630); mpfr_init(r25631); mpfr_init(r25632); mpfr_init_set_str(r25633, "4", 10, MPFR_RNDN); mpfr_init(r25634); mpfr_init(r25635); mpfr_init(r25636); mpfr_init(r25637); mpfr_init(r25638); mpfr_init(r25639); mpfr_init(r25640); mpfr_init(r25641); mpfr_init(r25642); mpfr_init(r25643); } double f_fm(double a, double b, double c) { mpfr_set_d(r25618, b, MPFR_RNDN); ; mpfr_set_si(r25620, mpfr_cmp(r25618, r25619) <= 0, MPFR_RNDN); mpfr_set_d(r25621, c, MPFR_RNDN); mpfr_div(r25622, r25621, r25618, MPFR_RNDN); mpfr_add(r25623, r25618, r25618, MPFR_RNDN); ; mpfr_set_d(r25625, a, MPFR_RNDN); mpfr_mul(r25626, r25624, r25625, MPFR_RNDN); mpfr_div(r25627, r25623, r25626, MPFR_RNDN); mpfr_sub(r25628, r25622, r25627, MPFR_RNDN); ; mpfr_set_si(r25630, mpfr_cmp(r25618, r25629) <= 0, MPFR_RNDN); mpfr_neg(r25631, r25618, MPFR_RNDN); mpfr_mul(r25632, r25618, r25618, MPFR_RNDN); ; mpfr_mul(r25634, r25633, r25625, MPFR_RNDN); mpfr_mul(r25635, r25621, r25634, MPFR_RNDN); mpfr_sub(r25636, r25632, r25635, MPFR_RNDN); mpfr_sqrt(r25637, r25636, MPFR_RNDN); mpfr_add(r25638, r25631, r25637, MPFR_RNDN); mpfr_div(r25639, r25638, r25626, MPFR_RNDN); mpfr_neg(r25640, r25621, MPFR_RNDN); mpfr_div(r25641, r25640, r25618, MPFR_RNDN); if (mpfr_get_si(r25630, MPFR_RNDN)) { mpfr_set(r25642, r25639, MPFR_RNDN); } else { mpfr_set(r25642, r25641, MPFR_RNDN); }; if (mpfr_get_si(r25620, MPFR_RNDN)) { mpfr_set(r25643, r25628, MPFR_RNDN); } else { mpfr_set(r25643, r25642, MPFR_RNDN); }; return mpfr_get_d(r25643, MPFR_RNDN); } static mpfr_t r25644, r25645, r25646, r25647, r25648, r25649, r25650, r25651, r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659, r25660, r25661, r25662, r25663, r25664, r25665, r25666, r25667, r25668, r25669; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init(r25644); mpfr_init_set_str(r25645, "-2.364349591187642e+101", 10, MPFR_RNDN); mpfr_init(r25646); mpfr_init(r25647); mpfr_init(r25648); mpfr_init(r25649); mpfr_init_set_str(r25650, "2", 10, MPFR_RNDN); mpfr_init(r25651); mpfr_init(r25652); mpfr_init(r25653); mpfr_init(r25654); mpfr_init_set_str(r25655, "1.125868897500658e-82", 10, MPFR_RNDN); mpfr_init(r25656); mpfr_init(r25657); mpfr_init(r25658); mpfr_init_set_str(r25659, "4", 10, MPFR_RNDN); mpfr_init(r25660); mpfr_init(r25661); mpfr_init(r25662); mpfr_init(r25663); mpfr_init(r25664); mpfr_init(r25665); mpfr_init(r25666); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); } double f_dm(double a, double b, double c) { mpfr_set_d(r25644, b, MPFR_RNDN); ; mpfr_set_si(r25646, mpfr_cmp(r25644, r25645) <= 0, MPFR_RNDN); mpfr_set_d(r25647, c, MPFR_RNDN); mpfr_div(r25648, r25647, r25644, MPFR_RNDN); mpfr_add(r25649, r25644, r25644, MPFR_RNDN); ; mpfr_set_d(r25651, a, MPFR_RNDN); mpfr_mul(r25652, r25650, r25651, MPFR_RNDN); mpfr_div(r25653, r25649, r25652, MPFR_RNDN); mpfr_sub(r25654, r25648, r25653, MPFR_RNDN); ; mpfr_set_si(r25656, mpfr_cmp(r25644, r25655) <= 0, MPFR_RNDN); mpfr_neg(r25657, r25644, MPFR_RNDN); mpfr_mul(r25658, r25644, r25644, MPFR_RNDN); ; mpfr_mul(r25660, r25659, r25651, MPFR_RNDN); mpfr_mul(r25661, r25647, r25660, MPFR_RNDN); mpfr_sub(r25662, r25658, r25661, MPFR_RNDN); mpfr_sqrt(r25663, r25662, MPFR_RNDN); mpfr_add(r25664, r25657, r25663, MPFR_RNDN); mpfr_div(r25665, r25664, r25652, MPFR_RNDN); mpfr_neg(r25666, r25647, MPFR_RNDN); mpfr_div(r25667, r25666, r25644, MPFR_RNDN); if (mpfr_get_si(r25656, MPFR_RNDN)) { mpfr_set(r25668, r25665, MPFR_RNDN); } else { mpfr_set(r25668, r25667, MPFR_RNDN); }; if (mpfr_get_si(r25646, MPFR_RNDN)) { mpfr_set(r25669, r25654, MPFR_RNDN); } else { mpfr_set(r25669, r25668, MPFR_RNDN); }; return mpfr_get_d(r25669, MPFR_RNDN); }