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