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