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