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