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