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