#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17551 = b; float r17552 = -r17551; float r17553 = r17551 * r17551; float r17554 = 4.0f; float r17555 = a; float r17556 = r17554 * r17555; float r17557 = c; float r17558 = r17556 * r17557; float r17559 = r17553 - r17558; float r17560 = sqrt(r17559); float r17561 = r17552 + r17560; float r17562 = 2.0f; float r17563 = r17562 * r17555; float r17564 = r17561 / r17563; return r17564; } double f_id(double a, double b, double c) { double r17565 = b; double r17566 = -r17565; double r17567 = r17565 * r17565; double r17568 = 4.0; double r17569 = a; double r17570 = r17568 * r17569; double r17571 = c; double r17572 = r17570 * r17571; double r17573 = r17567 - r17572; double r17574 = sqrt(r17573); double r17575 = r17566 + r17574; double r17576 = 2.0; double r17577 = r17576 * r17569; double r17578 = r17575 / r17577; return r17578; } double f_of(float a, float b, float c) { float r17579 = 4.0f; float r17580 = a; float r17581 = r17579 * r17580; float r17582 = c; float r17583 = r17581 * r17582; float r17584 = b; float r17585 = -r17584; float r17586 = r17584 * r17584; float r17587 = r17586 - r17583; float r17588 = sqrt(r17587); float r17589 = r17585 - r17588; float r17590 = r17583 / r17589; float r17591 = 2.0f; float r17592 = r17591 * r17580; float r17593 = r17590 / r17592; return r17593; } double f_od(double a, double b, double c) { double r17594 = 4.0; double r17595 = a; double r17596 = r17594 * r17595; double r17597 = c; double r17598 = r17596 * r17597; double r17599 = b; double r17600 = -r17599; double r17601 = r17599 * r17599; double r17602 = r17601 - r17598; double r17603 = sqrt(r17602); double r17604 = r17600 - r17603; double r17605 = r17598 / r17604; double r17606 = 2.0; double r17607 = r17606 * r17595; double r17608 = r17605 / r17607; return r17608; } 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 r17609, r17610, r17611, r17612, r17613, r17614, r17615, r17616, r17617, r17618, r17619, r17620, r17621, r17622; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17609); mpfr_init(r17610); mpfr_init(r17611); mpfr_init_set_str(r17612, "4", 10, MPFR_RNDN); mpfr_init(r17613); mpfr_init(r17614); mpfr_init(r17615); mpfr_init(r17616); mpfr_init(r17617); mpfr_init(r17618); mpfr_init(r17619); mpfr_init_set_str(r17620, "2", 10, MPFR_RNDN); mpfr_init(r17621); mpfr_init(r17622); } double f_im(double a, double b, double c) { mpfr_set_d(r17609, b, MPFR_RNDN); mpfr_neg(r17610, r17609, MPFR_RNDN); mpfr_sqr(r17611, r17609, MPFR_RNDN); ; mpfr_set_d(r17613, a, MPFR_RNDN); mpfr_mul(r17614, r17612, r17613, MPFR_RNDN); mpfr_set_d(r17615, c, MPFR_RNDN); mpfr_mul(r17616, r17614, r17615, MPFR_RNDN); mpfr_sub(r17617, r17611, r17616, MPFR_RNDN); mpfr_sqrt(r17618, r17617, MPFR_RNDN); mpfr_add(r17619, r17610, r17618, MPFR_RNDN); ; mpfr_mul(r17621, r17620, r17613, MPFR_RNDN); mpfr_div(r17622, r17619, r17621, MPFR_RNDN); return mpfr_get_d(r17622, MPFR_RNDN); } static mpfr_t r17623, r17624, r17625, r17626, r17627, r17628, r17629, r17630, r17631, r17632, r17633, r17634, r17635, r17636, r17637; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17623, "4", 10, MPFR_RNDN); mpfr_init(r17624); mpfr_init(r17625); mpfr_init(r17626); mpfr_init(r17627); mpfr_init(r17628); mpfr_init(r17629); mpfr_init(r17630); mpfr_init(r17631); mpfr_init(r17632); mpfr_init(r17633); mpfr_init(r17634); mpfr_init_set_str(r17635, "2", 10, MPFR_RNDN); mpfr_init(r17636); mpfr_init(r17637); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17624, a, MPFR_RNDN); mpfr_mul(r17625, r17623, r17624, MPFR_RNDN); mpfr_set_d(r17626, c, MPFR_RNDN); mpfr_mul(r17627, r17625, r17626, MPFR_RNDN); mpfr_set_d(r17628, b, MPFR_RNDN); mpfr_neg(r17629, r17628, MPFR_RNDN); mpfr_sqr(r17630, r17628, MPFR_RNDN); mpfr_sub(r17631, r17630, r17627, MPFR_RNDN); mpfr_sqrt(r17632, r17631, MPFR_RNDN); mpfr_sub(r17633, r17629, r17632, MPFR_RNDN); mpfr_div(r17634, r17627, r17633, MPFR_RNDN); ; mpfr_mul(r17636, r17635, r17624, MPFR_RNDN); mpfr_div(r17637, r17634, r17636, MPFR_RNDN); return mpfr_get_d(r17637, MPFR_RNDN); } static mpfr_t r17638, r17639, r17640, r17641, r17642, r17643, r17644, r17645, r17646, r17647, r17648, r17649, r17650, r17651, r17652; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17638, "4", 10, MPFR_RNDN); mpfr_init(r17639); mpfr_init(r17640); mpfr_init(r17641); mpfr_init(r17642); mpfr_init(r17643); mpfr_init(r17644); mpfr_init(r17645); mpfr_init(r17646); mpfr_init(r17647); mpfr_init(r17648); mpfr_init(r17649); mpfr_init_set_str(r17650, "2", 10, MPFR_RNDN); mpfr_init(r17651); mpfr_init(r17652); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17639, a, MPFR_RNDN); mpfr_mul(r17640, r17638, r17639, MPFR_RNDN); mpfr_set_d(r17641, c, MPFR_RNDN); mpfr_mul(r17642, r17640, r17641, MPFR_RNDN); mpfr_set_d(r17643, b, MPFR_RNDN); mpfr_neg(r17644, r17643, MPFR_RNDN); mpfr_sqr(r17645, r17643, MPFR_RNDN); mpfr_sub(r17646, r17645, r17642, MPFR_RNDN); mpfr_sqrt(r17647, r17646, MPFR_RNDN); mpfr_sub(r17648, r17644, r17647, MPFR_RNDN); mpfr_div(r17649, r17642, r17648, MPFR_RNDN); ; mpfr_mul(r17651, r17650, r17639, MPFR_RNDN); mpfr_div(r17652, r17649, r17651, MPFR_RNDN); return mpfr_get_d(r17652, MPFR_RNDN); }