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