#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17631 = b; float r17632 = -r17631; float r17633 = r17631 * r17631; float r17634 = 4.0f; float r17635 = a; float r17636 = r17634 * r17635; float r17637 = c; float r17638 = r17636 * r17637; float r17639 = r17633 - r17638; float r17640 = sqrt(r17639); float r17641 = r17632 + r17640; float r17642 = 2.0f; float r17643 = r17642 * r17635; float r17644 = r17641 / r17643; return r17644; } double f_id(double a, double b, double c) { double r17645 = b; double r17646 = -r17645; double r17647 = r17645 * r17645; double r17648 = 4.0; double r17649 = a; double r17650 = r17648 * r17649; double r17651 = c; double r17652 = r17650 * r17651; double r17653 = r17647 - r17652; double r17654 = sqrt(r17653); double r17655 = r17646 + r17654; double r17656 = 2.0; double r17657 = r17656 * r17649; double r17658 = r17655 / r17657; return r17658; } double f_of(float a, float b, float c) { float r17659 = 4.0f; float r17660 = a; float r17661 = r17659 * r17660; float r17662 = c; float r17663 = r17661 * r17662; float r17664 = b; float r17665 = -r17664; float r17666 = r17664 * r17664; float r17667 = r17666 - r17663; float r17668 = sqrt(r17667); float r17669 = r17665 - r17668; float r17670 = r17663 / r17669; float r17671 = 2.0f; float r17672 = r17671 * r17660; float r17673 = r17670 / r17672; return r17673; } double f_od(double a, double b, double c) { double r17674 = 4.0; double r17675 = a; double r17676 = r17674 * r17675; double r17677 = c; double r17678 = r17676 * r17677; double r17679 = b; double r17680 = -r17679; double r17681 = r17679 * r17679; double r17682 = r17681 - r17678; double r17683 = sqrt(r17682); double r17684 = r17680 - r17683; double r17685 = r17678 / r17684; double r17686 = 2.0; double r17687 = r17686 * r17675; double r17688 = r17685 / r17687; return r17688; } 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 r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); mpfr_init_set_str(r17692, "4", 10, MPFR_RNDN); mpfr_init(r17693); mpfr_init(r17694); mpfr_init(r17695); mpfr_init(r17696); mpfr_init(r17697); mpfr_init(r17698); mpfr_init(r17699); mpfr_init_set_str(r17700, "2", 10, MPFR_RNDN); mpfr_init(r17701); mpfr_init(r17702); } double f_im(double a, double b, double c) { mpfr_set_d(r17689, b, MPFR_RNDN); mpfr_neg(r17690, r17689, MPFR_RNDN); mpfr_sqr(r17691, r17689, MPFR_RNDN); ; mpfr_set_d(r17693, a, MPFR_RNDN); mpfr_mul(r17694, r17692, r17693, MPFR_RNDN); mpfr_set_d(r17695, c, MPFR_RNDN); mpfr_mul(r17696, r17694, r17695, MPFR_RNDN); mpfr_sub(r17697, r17691, r17696, MPFR_RNDN); mpfr_sqrt(r17698, r17697, MPFR_RNDN); mpfr_add(r17699, r17690, r17698, MPFR_RNDN); ; mpfr_mul(r17701, r17700, r17693, MPFR_RNDN); mpfr_div(r17702, r17699, r17701, MPFR_RNDN); return mpfr_get_d(r17702, MPFR_RNDN); } static mpfr_t r17703, r17704, r17705, r17706, r17707, r17708, r17709, r17710, r17711, r17712, r17713, r17714, r17715, r17716, r17717; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17703, "4", 10, MPFR_RNDN); mpfr_init(r17704); mpfr_init(r17705); mpfr_init(r17706); mpfr_init(r17707); mpfr_init(r17708); mpfr_init(r17709); mpfr_init(r17710); mpfr_init(r17711); mpfr_init(r17712); mpfr_init(r17713); mpfr_init(r17714); mpfr_init_set_str(r17715, "2", 10, MPFR_RNDN); mpfr_init(r17716); mpfr_init(r17717); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17704, a, MPFR_RNDN); mpfr_mul(r17705, r17703, r17704, MPFR_RNDN); mpfr_set_d(r17706, c, MPFR_RNDN); mpfr_mul(r17707, r17705, r17706, MPFR_RNDN); mpfr_set_d(r17708, b, MPFR_RNDN); mpfr_neg(r17709, r17708, MPFR_RNDN); mpfr_sqr(r17710, r17708, MPFR_RNDN); mpfr_sub(r17711, r17710, r17707, MPFR_RNDN); mpfr_sqrt(r17712, r17711, MPFR_RNDN); mpfr_sub(r17713, r17709, r17712, MPFR_RNDN); mpfr_div(r17714, r17707, r17713, MPFR_RNDN); ; mpfr_mul(r17716, r17715, r17704, MPFR_RNDN); mpfr_div(r17717, r17714, r17716, MPFR_RNDN); return mpfr_get_d(r17717, MPFR_RNDN); } static mpfr_t r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17718, "4", 10, MPFR_RNDN); mpfr_init(r17719); mpfr_init(r17720); mpfr_init(r17721); mpfr_init(r17722); mpfr_init(r17723); mpfr_init(r17724); mpfr_init(r17725); mpfr_init(r17726); mpfr_init(r17727); mpfr_init(r17728); mpfr_init(r17729); mpfr_init_set_str(r17730, "2", 10, MPFR_RNDN); mpfr_init(r17731); mpfr_init(r17732); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17719, a, MPFR_RNDN); mpfr_mul(r17720, r17718, r17719, MPFR_RNDN); mpfr_set_d(r17721, c, MPFR_RNDN); mpfr_mul(r17722, r17720, r17721, MPFR_RNDN); mpfr_set_d(r17723, b, MPFR_RNDN); mpfr_neg(r17724, r17723, MPFR_RNDN); mpfr_sqr(r17725, r17723, MPFR_RNDN); mpfr_sub(r17726, r17725, r17722, MPFR_RNDN); mpfr_sqrt(r17727, r17726, MPFR_RNDN); mpfr_sub(r17728, r17724, r17727, MPFR_RNDN); mpfr_div(r17729, r17722, r17728, MPFR_RNDN); ; mpfr_mul(r17731, r17730, r17719, MPFR_RNDN); mpfr_div(r17732, r17729, r17731, MPFR_RNDN); return mpfr_get_d(r17732, MPFR_RNDN); }