#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 = 1.0f; float r17612 = 2.0f; float r17613 = r17611 / r17612; float r17614 = 4.0f; float r17615 = c; float r17616 = r17614 * r17615; float r17617 = b; float r17618 = -r17617; float r17619 = r17617 * r17617; float r17620 = a; float r17621 = r17615 * r17620; float r17622 = r17621 * r17614; float r17623 = r17619 - r17622; float r17624 = sqrt(r17623); float r17625 = r17618 - r17624; float r17626 = r17616 / r17625; float r17627 = r17613 * r17626; return r17627; } double f_od(double a, double b, double c) { double r17628 = 1.0; double r17629 = 2.0; double r17630 = r17628 / r17629; double r17631 = 4.0; double r17632 = c; double r17633 = r17631 * r17632; double r17634 = b; double r17635 = -r17634; double r17636 = r17634 * r17634; double r17637 = a; double r17638 = r17632 * r17637; double r17639 = r17638 * r17631; double r17640 = r17636 - r17639; double r17641 = sqrt(r17640); double r17642 = r17635 - r17641; double r17643 = r17633 / r17642; double r17644 = r17630 * r17643; return r17644; } 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 r17645, r17646, r17647, r17648, r17649, r17650, r17651, r17652, r17653, r17654, r17655, r17656, r17657, r17658; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17645); mpfr_init(r17646); mpfr_init(r17647); mpfr_init_set_str(r17648, "4", 10, MPFR_RNDN); mpfr_init(r17649); mpfr_init(r17650); mpfr_init(r17651); mpfr_init(r17652); mpfr_init(r17653); mpfr_init(r17654); mpfr_init(r17655); mpfr_init_set_str(r17656, "2", 10, MPFR_RNDN); mpfr_init(r17657); mpfr_init(r17658); } double f_im(double a, double b, double c) { mpfr_set_d(r17645, b, MPFR_RNDN); mpfr_neg(r17646, r17645, MPFR_RNDN); mpfr_sqr(r17647, r17645, MPFR_RNDN); ; mpfr_set_d(r17649, a, MPFR_RNDN); mpfr_mul(r17650, r17648, r17649, MPFR_RNDN); mpfr_set_d(r17651, c, MPFR_RNDN); mpfr_mul(r17652, r17650, r17651, MPFR_RNDN); mpfr_sub(r17653, r17647, r17652, MPFR_RNDN); mpfr_sqrt(r17654, r17653, MPFR_RNDN); mpfr_add(r17655, r17646, r17654, MPFR_RNDN); ; mpfr_mul(r17657, r17656, r17649, MPFR_RNDN); mpfr_div(r17658, r17655, r17657, MPFR_RNDN); return mpfr_get_d(r17658, MPFR_RNDN); } static mpfr_t r17659, r17660, r17661, r17662, r17663, r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17659, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17660, "2", 10, MPFR_RNDN); mpfr_init(r17661); mpfr_init_set_str(r17662, "4", 10, MPFR_RNDN); mpfr_init(r17663); mpfr_init(r17664); mpfr_init(r17665); mpfr_init(r17666); mpfr_init(r17667); mpfr_init(r17668); mpfr_init(r17669); mpfr_init(r17670); mpfr_init(r17671); mpfr_init(r17672); mpfr_init(r17673); mpfr_init(r17674); mpfr_init(r17675); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r17661, r17659, r17660, MPFR_RNDN); ; mpfr_set_d(r17663, c, MPFR_RNDN); mpfr_mul(r17664, r17662, r17663, MPFR_RNDN); mpfr_set_d(r17665, b, MPFR_RNDN); mpfr_neg(r17666, r17665, MPFR_RNDN); mpfr_sqr(r17667, r17665, MPFR_RNDN); mpfr_set_d(r17668, a, MPFR_RNDN); mpfr_mul(r17669, r17663, r17668, MPFR_RNDN); mpfr_mul(r17670, r17669, r17662, MPFR_RNDN); mpfr_sub(r17671, r17667, r17670, MPFR_RNDN); mpfr_sqrt(r17672, r17671, MPFR_RNDN); mpfr_sub(r17673, r17666, r17672, MPFR_RNDN); mpfr_div(r17674, r17664, r17673, MPFR_RNDN); mpfr_mul(r17675, r17661, r17674, MPFR_RNDN); return mpfr_get_d(r17675, MPFR_RNDN); } static mpfr_t r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691, r17692; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17676, "1", 10, MPFR_RNDN); mpfr_init_set_str(r17677, "2", 10, MPFR_RNDN); mpfr_init(r17678); mpfr_init_set_str(r17679, "4", 10, MPFR_RNDN); mpfr_init(r17680); mpfr_init(r17681); mpfr_init(r17682); mpfr_init(r17683); mpfr_init(r17684); mpfr_init(r17685); mpfr_init(r17686); mpfr_init(r17687); mpfr_init(r17688); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); mpfr_init(r17692); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r17678, r17676, r17677, MPFR_RNDN); ; mpfr_set_d(r17680, c, MPFR_RNDN); mpfr_mul(r17681, r17679, r17680, MPFR_RNDN); mpfr_set_d(r17682, b, MPFR_RNDN); mpfr_neg(r17683, r17682, MPFR_RNDN); mpfr_sqr(r17684, r17682, MPFR_RNDN); mpfr_set_d(r17685, a, MPFR_RNDN); mpfr_mul(r17686, r17680, r17685, MPFR_RNDN); mpfr_mul(r17687, r17686, r17679, MPFR_RNDN); mpfr_sub(r17688, r17684, r17687, MPFR_RNDN); mpfr_sqrt(r17689, r17688, MPFR_RNDN); mpfr_sub(r17690, r17683, r17689, MPFR_RNDN); mpfr_div(r17691, r17681, r17690, MPFR_RNDN); mpfr_mul(r17692, r17678, r17691, MPFR_RNDN); return mpfr_get_d(r17692, MPFR_RNDN); }