#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17655 = b; float r17656 = -r17655; float r17657 = r17655 * r17655; float r17658 = 4.0f; float r17659 = a; float r17660 = r17658 * r17659; float r17661 = c; float r17662 = r17660 * r17661; float r17663 = r17657 - r17662; float r17664 = sqrt(r17663); float r17665 = r17656 + r17664; float r17666 = 2.0f; float r17667 = r17666 * r17659; float r17668 = r17665 / r17667; return r17668; } double f_id(double a, double b, double c) { double r17669 = b; double r17670 = -r17669; double r17671 = r17669 * r17669; double r17672 = 4.0; double r17673 = a; double r17674 = r17672 * r17673; double r17675 = c; double r17676 = r17674 * r17675; double r17677 = r17671 - r17676; double r17678 = sqrt(r17677); double r17679 = r17670 + r17678; double r17680 = 2.0; double r17681 = r17680 * r17673; double r17682 = r17679 / r17681; return r17682; } double f_of(float a, float b, float c) { float r17683 = 4.0f; float r17684 = a; float r17685 = r17683 * r17684; float r17686 = 1.0f; float r17687 = r17685 / r17686; float r17688 = c; float r17689 = b; float r17690 = -r17689; float r17691 = r17689 * r17689; float r17692 = r17685 * r17688; float r17693 = r17691 - r17692; float r17694 = sqrt(r17693); float r17695 = r17690 - r17694; float r17696 = r17688 / r17695; float r17697 = r17687 * r17696; float r17698 = 2.0f; float r17699 = r17698 * r17684; float r17700 = r17697 / r17699; return r17700; } double f_od(double a, double b, double c) { double r17701 = 4.0; double r17702 = a; double r17703 = r17701 * r17702; double r17704 = 1.0; double r17705 = r17703 / r17704; double r17706 = c; double r17707 = b; double r17708 = -r17707; double r17709 = r17707 * r17707; double r17710 = r17703 * r17706; double r17711 = r17709 - r17710; double r17712 = sqrt(r17711); double r17713 = r17708 - r17712; double r17714 = r17706 / r17713; double r17715 = r17705 * r17714; double r17716 = 2.0; double r17717 = r17716 * r17702; double r17718 = r17715 / r17717; return r17718; } 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 r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17719); mpfr_init(r17720); mpfr_init(r17721); mpfr_init_set_str(r17722, "4", 10, MPFR_RNDN); 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_im(double a, double b, double c) { mpfr_set_d(r17719, b, MPFR_RNDN); mpfr_neg(r17720, r17719, MPFR_RNDN); mpfr_sqr(r17721, r17719, MPFR_RNDN); ; mpfr_set_d(r17723, a, MPFR_RNDN); mpfr_mul(r17724, r17722, r17723, MPFR_RNDN); mpfr_set_d(r17725, c, MPFR_RNDN); mpfr_mul(r17726, r17724, r17725, MPFR_RNDN); mpfr_sub(r17727, r17721, r17726, MPFR_RNDN); mpfr_sqrt(r17728, r17727, MPFR_RNDN); mpfr_add(r17729, r17720, r17728, MPFR_RNDN); ; mpfr_mul(r17731, r17730, r17723, MPFR_RNDN); mpfr_div(r17732, r17729, r17731, MPFR_RNDN); return mpfr_get_d(r17732, MPFR_RNDN); } static mpfr_t r17733, r17734, r17735, r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743, r17744, r17745, r17746, r17747, r17748, r17749, r17750; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17733, "4", 10, MPFR_RNDN); mpfr_init(r17734); mpfr_init(r17735); mpfr_init_set_str(r17736, "1", 10, MPFR_RNDN); mpfr_init(r17737); mpfr_init(r17738); mpfr_init(r17739); mpfr_init(r17740); mpfr_init(r17741); mpfr_init(r17742); mpfr_init(r17743); mpfr_init(r17744); mpfr_init(r17745); mpfr_init(r17746); mpfr_init(r17747); mpfr_init_set_str(r17748, "2", 10, MPFR_RNDN); mpfr_init(r17749); mpfr_init(r17750); } double f_fm(double a, double b, double c) { ; mpfr_set_d(r17734, a, MPFR_RNDN); mpfr_mul(r17735, r17733, r17734, MPFR_RNDN); ; mpfr_div(r17737, r17735, r17736, MPFR_RNDN); mpfr_set_d(r17738, c, MPFR_RNDN); mpfr_set_d(r17739, b, MPFR_RNDN); mpfr_neg(r17740, r17739, MPFR_RNDN); mpfr_sqr(r17741, r17739, MPFR_RNDN); mpfr_mul(r17742, r17735, r17738, MPFR_RNDN); mpfr_sub(r17743, r17741, r17742, MPFR_RNDN); mpfr_sqrt(r17744, r17743, MPFR_RNDN); mpfr_sub(r17745, r17740, r17744, MPFR_RNDN); mpfr_div(r17746, r17738, r17745, MPFR_RNDN); mpfr_mul(r17747, r17737, r17746, MPFR_RNDN); ; mpfr_mul(r17749, r17748, r17734, MPFR_RNDN); mpfr_div(r17750, r17747, r17749, MPFR_RNDN); return mpfr_get_d(r17750, MPFR_RNDN); } static mpfr_t r17751, r17752, r17753, r17754, r17755, r17756, r17757, r17758, r17759, r17760, r17761, r17762, r17763, r17764, r17765, r17766, r17767, r17768; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17751, "4", 10, MPFR_RNDN); mpfr_init(r17752); mpfr_init(r17753); mpfr_init_set_str(r17754, "1", 10, MPFR_RNDN); mpfr_init(r17755); mpfr_init(r17756); mpfr_init(r17757); mpfr_init(r17758); mpfr_init(r17759); mpfr_init(r17760); mpfr_init(r17761); mpfr_init(r17762); mpfr_init(r17763); mpfr_init(r17764); mpfr_init(r17765); mpfr_init_set_str(r17766, "2", 10, MPFR_RNDN); mpfr_init(r17767); mpfr_init(r17768); } double f_dm(double a, double b, double c) { ; mpfr_set_d(r17752, a, MPFR_RNDN); mpfr_mul(r17753, r17751, r17752, MPFR_RNDN); ; mpfr_div(r17755, r17753, r17754, MPFR_RNDN); mpfr_set_d(r17756, c, MPFR_RNDN); mpfr_set_d(r17757, b, MPFR_RNDN); mpfr_neg(r17758, r17757, MPFR_RNDN); mpfr_sqr(r17759, r17757, MPFR_RNDN); mpfr_mul(r17760, r17753, r17756, MPFR_RNDN); mpfr_sub(r17761, r17759, r17760, MPFR_RNDN); mpfr_sqrt(r17762, r17761, MPFR_RNDN); mpfr_sub(r17763, r17758, r17762, MPFR_RNDN); mpfr_div(r17764, r17756, r17763, MPFR_RNDN); mpfr_mul(r17765, r17755, r17764, MPFR_RNDN); ; mpfr_mul(r17767, r17766, r17752, MPFR_RNDN); mpfr_div(r17768, r17765, r17767, MPFR_RNDN); return mpfr_get_d(r17768, MPFR_RNDN); }