#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17583 = b; float r17584 = -r17583; float r17585 = r17583 * r17583; float r17586 = 4.0f; float r17587 = a; float r17588 = r17586 * r17587; float r17589 = c; float r17590 = r17588 * r17589; float r17591 = r17585 - r17590; float r17592 = sqrt(r17591); float r17593 = r17584 + r17592; float r17594 = 2.0f; float r17595 = r17594 * r17587; float r17596 = r17593 / r17595; return r17596; } double f_id(double a, double b, double c) { double r17597 = b; double r17598 = -r17597; double r17599 = r17597 * r17597; double r17600 = 4.0; double r17601 = a; double r17602 = r17600 * r17601; double r17603 = c; double r17604 = r17602 * r17603; double r17605 = r17599 - r17604; double r17606 = sqrt(r17605); double r17607 = r17598 + r17606; double r17608 = 2.0; double r17609 = r17608 * r17601; double r17610 = r17607 / r17609; return r17610; } double f_of(float a, float b, float c) { float r17611 = 4.0f; float r17612 = a; float r17613 = r17611 * r17612; float r17614 = c; float r17615 = r17613 * r17614; float r17616 = b; float r17617 = -r17616; float r17618 = r17616 * r17616; float r17619 = r17618 - r17615; float r17620 = sqrt(r17619); float r17621 = r17617 - r17620; float r17622 = r17615 / r17621; float r17623 = 2.0f; float r17624 = r17623 * r17612; float r17625 = r17622 / r17624; return r17625; } double f_od(double a, double b, double c) { double r17626 = 4.0; double r17627 = a; double r17628 = r17626 * r17627; double r17629 = c; double r17630 = r17628 * r17629; double r17631 = b; double r17632 = -r17631; double r17633 = r17631 * r17631; double r17634 = r17633 - r17630; double r17635 = sqrt(r17634); double r17636 = r17632 - r17635; double r17637 = r17630 / r17636; double r17638 = 2.0; double r17639 = r17638 * r17627; double r17640 = r17637 / r17639; return r17640; } 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 r17641, r17642, r17643, r17644, r17645, r17646, r17647, r17648, r17649, r17650, r17651, r17652, r17653, r17654; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17641); mpfr_init(r17642); mpfr_init(r17643); mpfr_init_set_str(r17644, "4", 10, MPFR_RNDN); mpfr_init(r17645); mpfr_init(r17646); mpfr_init(r17647); mpfr_init(r17648); mpfr_init(r17649); mpfr_init(r17650); mpfr_init(r17651); mpfr_init_set_str(r17652, "2", 10, MPFR_RNDN); mpfr_init(r17653); mpfr_init(r17654); } double f_im(double a, double b, double c) { mpfr_set_d(r17641, b, MPFR_RNDN); mpfr_neg(r17642, r17641, MPFR_RNDN); mpfr_sqr(r17643, r17641, MPFR_RNDN); ; mpfr_set_d(r17645, a, MPFR_RNDN); mpfr_mul(r17646, r17644, r17645, MPFR_RNDN); mpfr_set_d(r17647, c, MPFR_RNDN); mpfr_mul(r17648, r17646, r17647, MPFR_RNDN); mpfr_sub(r17649, r17643, r17648, MPFR_RNDN); mpfr_sqrt(r17650, r17649, MPFR_RNDN); mpfr_add(r17651, r17642, r17650, MPFR_RNDN); ; mpfr_mul(r17653, r17652, r17645, MPFR_RNDN); mpfr_div(r17654, r17651, r17653, MPFR_RNDN); return mpfr_get_d(r17654, MPFR_RNDN); } static mpfr_t r17655, r17656, r17657, r17658, r17659, r17660, r17661, r17662, r17663, r17664, r17665, r17666, r17667, r17668, r17669; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17655, "4", 10, MPFR_RNDN); mpfr_init(r17656); mpfr_init(r17657); mpfr_init(r17658); mpfr_init(r17659); mpfr_init(r17660); mpfr_init(r17661); mpfr_init(r17662); mpfr_init(r17663); mpfr_init(r17664); mpfr_init(r17665); mpfr_init(r17666); mpfr_init_set_str(r17667, "2", 10, MPFR_RNDN); mpfr_init(r17668); mpfr_init(r17669); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17656, a, MPFR_RNDN); mpfr_mul(r17657, r17655, r17656, MPFR_RNDN); mpfr_set_d(r17658, c, MPFR_RNDN); mpfr_mul(r17659, r17657, r17658, MPFR_RNDN); mpfr_set_d(r17660, b, MPFR_RNDN); mpfr_neg(r17661, r17660, MPFR_RNDN); mpfr_sqr(r17662, r17660, MPFR_RNDN); mpfr_sub(r17663, r17662, r17659, MPFR_RNDN); mpfr_sqrt(r17664, r17663, MPFR_RNDN); mpfr_sub(r17665, r17661, r17664, MPFR_RNDN); mpfr_div(r17666, r17659, r17665, MPFR_RNDN); ; mpfr_mul(r17668, r17667, r17656, MPFR_RNDN); mpfr_div(r17669, r17666, r17668, MPFR_RNDN); return mpfr_get_d(r17669, MPFR_RNDN); } static mpfr_t r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17670, "4", 10, MPFR_RNDN); mpfr_init(r17671); mpfr_init(r17672); mpfr_init(r17673); mpfr_init(r17674); mpfr_init(r17675); mpfr_init(r17676); mpfr_init(r17677); mpfr_init(r17678); mpfr_init(r17679); mpfr_init(r17680); mpfr_init(r17681); mpfr_init_set_str(r17682, "2", 10, MPFR_RNDN); mpfr_init(r17683); mpfr_init(r17684); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17671, a, MPFR_RNDN); mpfr_mul(r17672, r17670, r17671, MPFR_RNDN); mpfr_set_d(r17673, c, MPFR_RNDN); mpfr_mul(r17674, r17672, r17673, MPFR_RNDN); mpfr_set_d(r17675, b, MPFR_RNDN); mpfr_neg(r17676, r17675, MPFR_RNDN); mpfr_sqr(r17677, r17675, MPFR_RNDN); mpfr_sub(r17678, r17677, r17674, MPFR_RNDN); mpfr_sqrt(r17679, r17678, MPFR_RNDN); mpfr_sub(r17680, r17676, r17679, MPFR_RNDN); mpfr_div(r17681, r17674, r17680, MPFR_RNDN); ; mpfr_mul(r17683, r17682, r17671, MPFR_RNDN); mpfr_div(r17684, r17681, r17683, MPFR_RNDN); return mpfr_get_d(r17684, MPFR_RNDN); }