#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r17558 = b; float r17559 = -r17558; float r17560 = r17558 * r17558; float r17561 = 4.0f; float r17562 = a; float r17563 = r17561 * r17562; float r17564 = c; float r17565 = r17563 * r17564; float r17566 = r17560 - r17565; float r17567 = sqrt(r17566); float r17568 = r17559 + r17567; float r17569 = 2.0f; float r17570 = r17569 * r17562; float r17571 = r17568 / r17570; return r17571; } double f_id(double a, double b, double c) { double r17572 = b; double r17573 = -r17572; double r17574 = r17572 * r17572; double r17575 = 4.0; double r17576 = a; double r17577 = r17575 * r17576; double r17578 = c; double r17579 = r17577 * r17578; double r17580 = r17574 - r17579; double r17581 = sqrt(r17580); double r17582 = r17573 + r17581; double r17583 = 2.0; double r17584 = r17583 * r17576; double r17585 = r17582 / r17584; return r17585; } double f_of(float a, float b, float c) { float r17586 = b; float r17587 = -1.5511731923780108e+19f; bool r17588 = r17586 <= r17587; float r17589 = c; float r17590 = r17589 / r17586; float r17591 = a; float r17592 = r17586 / r17591; float r17593 = r17590 - r17592; float r17594 = 4.954544863267074e-08f; bool r17595 = r17586 <= r17594; float r17596 = -r17586; float r17597 = r17586 * r17586; float r17598 = 4.0f; float r17599 = r17598 * r17591; float r17600 = r17599 * r17589; float r17601 = r17597 - r17600; float r17602 = sqrt(r17601); float r17603 = r17596 + r17602; float r17604 = 2.0f; float r17605 = r17604 * r17591; float r17606 = r17603 / r17605; float r17607 = -2.0f; float r17608 = r17607 / r17604; float r17609 = r17590 * r17608; float r17610 = r17595 ? r17606 : r17609; float r17611 = r17588 ? r17593 : r17610; return r17611; } double f_od(double a, double b, double c) { double r17612 = b; double r17613 = -1.5511731923780108e+19; bool r17614 = r17612 <= r17613; double r17615 = c; double r17616 = r17615 / r17612; double r17617 = a; double r17618 = r17612 / r17617; double r17619 = r17616 - r17618; double r17620 = 4.954544863267074e-08; bool r17621 = r17612 <= r17620; double r17622 = -r17612; double r17623 = r17612 * r17612; double r17624 = 4.0; double r17625 = r17624 * r17617; double r17626 = r17625 * r17615; double r17627 = r17623 - r17626; double r17628 = sqrt(r17627); double r17629 = r17622 + r17628; double r17630 = 2.0; double r17631 = r17630 * r17617; double r17632 = r17629 / r17631; double r17633 = -2.0; double r17634 = r17633 / r17630; double r17635 = r17616 * r17634; double r17636 = r17621 ? r17632 : r17635; double r17637 = r17614 ? r17619 : r17636; return r17637; } 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 r17638, r17639, r17640, r17641, r17642, r17643, r17644, r17645, r17646, r17647, r17648, r17649, r17650, r17651; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17638); mpfr_init(r17639); mpfr_init(r17640); mpfr_init_set_str(r17641, "4", 10, MPFR_RNDN); mpfr_init(r17642); mpfr_init(r17643); mpfr_init(r17644); mpfr_init(r17645); mpfr_init(r17646); mpfr_init(r17647); mpfr_init(r17648); mpfr_init_set_str(r17649, "2", 10, MPFR_RNDN); mpfr_init(r17650); mpfr_init(r17651); } double f_im(double a, double b, double c) { mpfr_set_d(r17638, b, MPFR_RNDN); mpfr_neg(r17639, r17638, MPFR_RNDN); mpfr_sqr(r17640, r17638, MPFR_RNDN); ; mpfr_set_d(r17642, a, MPFR_RNDN); mpfr_mul(r17643, r17641, r17642, MPFR_RNDN); mpfr_set_d(r17644, c, MPFR_RNDN); mpfr_mul(r17645, r17643, r17644, MPFR_RNDN); mpfr_sub(r17646, r17640, r17645, MPFR_RNDN); mpfr_sqrt(r17647, r17646, MPFR_RNDN); mpfr_add(r17648, r17639, r17647, MPFR_RNDN); ; mpfr_mul(r17650, r17649, r17642, MPFR_RNDN); mpfr_div(r17651, r17648, r17650, MPFR_RNDN); return mpfr_get_d(r17651, MPFR_RNDN); } static mpfr_t r17652, r17653, r17654, r17655, r17656, r17657, r17658, r17659, r17660, r17661, r17662, r17663, r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17652); mpfr_init_set_str(r17653, "-1.5511732f+19", 10, MPFR_RNDN); mpfr_init(r17654); mpfr_init(r17655); mpfr_init(r17656); mpfr_init(r17657); mpfr_init(r17658); mpfr_init(r17659); mpfr_init_set_str(r17660, "4.954545f-08", 10, MPFR_RNDN); mpfr_init(r17661); mpfr_init(r17662); mpfr_init(r17663); mpfr_init_set_str(r17664, "4", 10, MPFR_RNDN); 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); mpfr_init_set_str(r17673, "-2", 10, MPFR_RNDN); mpfr_init(r17674); mpfr_init(r17675); mpfr_init(r17676); mpfr_init(r17677); } double f_fm(double a, double b, double c) { mpfr_set_d(r17652, b, MPFR_RNDN); ; mpfr_set_si(r17654, mpfr_cmp(r17652, r17653) <= 0, MPFR_RNDN); mpfr_set_d(r17655, c, MPFR_RNDN); mpfr_div(r17656, r17655, r17652, MPFR_RNDN); mpfr_set_d(r17657, a, MPFR_RNDN); mpfr_div(r17658, r17652, r17657, MPFR_RNDN); mpfr_sub(r17659, r17656, r17658, MPFR_RNDN); ; mpfr_set_si(r17661, mpfr_cmp(r17652, r17660) <= 0, MPFR_RNDN); mpfr_neg(r17662, r17652, MPFR_RNDN); mpfr_sqr(r17663, r17652, MPFR_RNDN); ; mpfr_mul(r17665, r17664, r17657, MPFR_RNDN); mpfr_mul(r17666, r17665, r17655, MPFR_RNDN); mpfr_sub(r17667, r17663, r17666, MPFR_RNDN); mpfr_sqrt(r17668, r17667, MPFR_RNDN); mpfr_add(r17669, r17662, r17668, MPFR_RNDN); ; mpfr_mul(r17671, r17670, r17657, MPFR_RNDN); mpfr_div(r17672, r17669, r17671, MPFR_RNDN); ; mpfr_div(r17674, r17673, r17670, MPFR_RNDN); mpfr_mul(r17675, r17656, r17674, MPFR_RNDN); if (mpfr_get_si(r17661, MPFR_RNDN)) { mpfr_set(r17676, r17672, MPFR_RNDN); } else { mpfr_set(r17676, r17675, MPFR_RNDN); }; if (mpfr_get_si(r17654, MPFR_RNDN)) { mpfr_set(r17677, r17659, MPFR_RNDN); } else { mpfr_set(r17677, r17676, MPFR_RNDN); }; return mpfr_get_d(r17677, MPFR_RNDN); } static mpfr_t r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17678); mpfr_init_set_str(r17679, "-1.5511732f+19", 10, MPFR_RNDN); mpfr_init(r17680); mpfr_init(r17681); mpfr_init(r17682); mpfr_init(r17683); mpfr_init(r17684); mpfr_init(r17685); mpfr_init_set_str(r17686, "4.954545f-08", 10, MPFR_RNDN); mpfr_init(r17687); mpfr_init(r17688); mpfr_init(r17689); mpfr_init_set_str(r17690, "4", 10, MPFR_RNDN); mpfr_init(r17691); mpfr_init(r17692); mpfr_init(r17693); mpfr_init(r17694); mpfr_init(r17695); mpfr_init_set_str(r17696, "2", 10, MPFR_RNDN); mpfr_init(r17697); mpfr_init(r17698); mpfr_init_set_str(r17699, "-2", 10, MPFR_RNDN); mpfr_init(r17700); mpfr_init(r17701); mpfr_init(r17702); mpfr_init(r17703); } double f_dm(double a, double b, double c) { mpfr_set_d(r17678, b, MPFR_RNDN); ; mpfr_set_si(r17680, mpfr_cmp(r17678, r17679) <= 0, MPFR_RNDN); mpfr_set_d(r17681, c, MPFR_RNDN); mpfr_div(r17682, r17681, r17678, MPFR_RNDN); mpfr_set_d(r17683, a, MPFR_RNDN); mpfr_div(r17684, r17678, r17683, MPFR_RNDN); mpfr_sub(r17685, r17682, r17684, MPFR_RNDN); ; mpfr_set_si(r17687, mpfr_cmp(r17678, r17686) <= 0, MPFR_RNDN); mpfr_neg(r17688, r17678, MPFR_RNDN); mpfr_sqr(r17689, r17678, MPFR_RNDN); ; mpfr_mul(r17691, r17690, r17683, MPFR_RNDN); mpfr_mul(r17692, r17691, r17681, MPFR_RNDN); mpfr_sub(r17693, r17689, r17692, MPFR_RNDN); mpfr_sqrt(r17694, r17693, MPFR_RNDN); mpfr_add(r17695, r17688, r17694, MPFR_RNDN); ; mpfr_mul(r17697, r17696, r17683, MPFR_RNDN); mpfr_div(r17698, r17695, r17697, MPFR_RNDN); ; mpfr_div(r17700, r17699, r17696, MPFR_RNDN); mpfr_mul(r17701, r17682, r17700, MPFR_RNDN); if (mpfr_get_si(r17687, MPFR_RNDN)) { mpfr_set(r17702, r17698, MPFR_RNDN); } else { mpfr_set(r17702, r17701, MPFR_RNDN); }; if (mpfr_get_si(r17680, MPFR_RNDN)) { mpfr_set(r17703, r17685, MPFR_RNDN); } else { mpfr_set(r17703, r17702, MPFR_RNDN); }; return mpfr_get_d(r17703, MPFR_RNDN); }