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