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