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