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