#include #include #include #include #include char *name = "The quadratic formula (r2)"; double f_if(float a, float b, float c) { float r22580 = b; float r22581 = -r22580; float r22582 = r22580 * r22580; float r22583 = 4; float r22584 = a; float r22585 = c; float r22586 = r22584 * r22585; float r22587 = r22583 * r22586; float r22588 = r22582 - r22587; float r22589 = sqrt(r22588); float r22590 = r22581 - r22589; float r22591 = 2; float r22592 = r22591 * r22584; float r22593 = r22590 / r22592; return r22593; } double f_id(double a, double b, double c) { double r22594 = b; double r22595 = -r22594; double r22596 = r22594 * r22594; double r22597 = 4; double r22598 = a; double r22599 = c; double r22600 = r22598 * r22599; double r22601 = r22597 * r22600; double r22602 = r22596 - r22601; double r22603 = sqrt(r22602); double r22604 = r22595 - r22603; double r22605 = 2; double r22606 = r22605 * r22598; double r22607 = r22604 / r22606; return r22607; } double f_of(float a, float b, float c) { float r22608 = b; float r22609 = -1.3775702482140898e-117; bool r22610 = r22608 <= r22609; float r22611 = c; float r22612 = -r22611; float r22613 = 1; float r22614 = r22608 / r22613; float r22615 = r22612 / r22614; float r22616 = 3.2254415904266855e+123; bool r22617 = r22608 <= r22616; float r22618 = 2; float r22619 = a; float r22620 = r22618 * r22619; float r22621 = -r22608; float r22622 = r22608 * r22608; float r22623 = 4; float r22624 = r22619 * r22611; float r22625 = r22623 * r22624; float r22626 = r22622 - r22625; float r22627 = sqrt(r22626); float r22628 = r22621 - r22627; float r22629 = r22620 / r22628; float r22630 = r22613 / r22629; float r22631 = r22621 / r22619; float r22632 = r22617 ? r22630 : r22631; float r22633 = r22610 ? r22615 : r22632; return r22633; } double f_od(double a, double b, double c) { double r22634 = b; double r22635 = -1.3775702482140898e-117; bool r22636 = r22634 <= r22635; double r22637 = c; double r22638 = -r22637; double r22639 = 1; double r22640 = r22634 / r22639; double r22641 = r22638 / r22640; double r22642 = 3.2254415904266855e+123; bool r22643 = r22634 <= r22642; double r22644 = 2; double r22645 = a; double r22646 = r22644 * r22645; double r22647 = -r22634; double r22648 = r22634 * r22634; double r22649 = 4; double r22650 = r22645 * r22637; double r22651 = r22649 * r22650; double r22652 = r22648 - r22651; double r22653 = sqrt(r22652); double r22654 = r22647 - r22653; double r22655 = r22646 / r22654; double r22656 = r22639 / r22655; double r22657 = r22647 / r22645; double r22658 = r22643 ? r22656 : r22657; double r22659 = r22636 ? r22641 : r22658; return r22659; } 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 r22660, r22661, r22662, r22663, r22664, r22665, r22666, r22667, r22668, r22669, r22670, r22671, r22672, r22673; void setup_mpfr_f_im() { mpfr_set_default_prec(3216); mpfr_init(r22660); mpfr_init(r22661); mpfr_init(r22662); mpfr_init_set_str(r22663, "4", 10, MPFR_RNDN); mpfr_init(r22664); mpfr_init(r22665); mpfr_init(r22666); mpfr_init(r22667); mpfr_init(r22668); mpfr_init(r22669); mpfr_init(r22670); mpfr_init_set_str(r22671, "2", 10, MPFR_RNDN); mpfr_init(r22672); mpfr_init(r22673); } double f_im(double a, double b, double c) { mpfr_set_d(r22660, b, MPFR_RNDN); mpfr_neg(r22661, r22660, MPFR_RNDN); mpfr_mul(r22662, r22660, r22660, MPFR_RNDN); ; mpfr_set_d(r22664, a, MPFR_RNDN); mpfr_set_d(r22665, c, MPFR_RNDN); mpfr_mul(r22666, r22664, r22665, MPFR_RNDN); mpfr_mul(r22667, r22663, r22666, MPFR_RNDN); mpfr_sub(r22668, r22662, r22667, MPFR_RNDN); mpfr_sqrt(r22669, r22668, MPFR_RNDN); mpfr_sub(r22670, r22661, r22669, MPFR_RNDN); ; mpfr_mul(r22672, r22671, r22664, MPFR_RNDN); mpfr_div(r22673, r22670, r22672, MPFR_RNDN); return mpfr_get_d(r22673, MPFR_RNDN); } static mpfr_t r22674, r22675, r22676, r22677, r22678, r22679, r22680, r22681, r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689, r22690, r22691, r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699; void setup_mpfr_f_fm() { mpfr_set_default_prec(3216); mpfr_init(r22674); mpfr_init_set_str(r22675, "-1.3775702482140898e-117", 10, MPFR_RNDN); mpfr_init(r22676); mpfr_init(r22677); mpfr_init(r22678); mpfr_init_set_str(r22679, "1", 10, MPFR_RNDN); mpfr_init(r22680); mpfr_init(r22681); mpfr_init_set_str(r22682, "3.2254415904266855e+123", 10, MPFR_RNDN); mpfr_init(r22683); mpfr_init_set_str(r22684, "2", 10, MPFR_RNDN); mpfr_init(r22685); mpfr_init(r22686); mpfr_init(r22687); mpfr_init(r22688); mpfr_init_set_str(r22689, "4", 10, MPFR_RNDN); mpfr_init(r22690); mpfr_init(r22691); mpfr_init(r22692); mpfr_init(r22693); mpfr_init(r22694); mpfr_init(r22695); mpfr_init(r22696); mpfr_init(r22697); mpfr_init(r22698); mpfr_init(r22699); } double f_fm(double a, double b, double c) { mpfr_set_d(r22674, b, MPFR_RNDN); ; mpfr_set_si(r22676, mpfr_cmp(r22674, r22675) <= 0, MPFR_RNDN); mpfr_set_d(r22677, c, MPFR_RNDN); mpfr_neg(r22678, r22677, MPFR_RNDN); ; mpfr_div(r22680, r22674, r22679, MPFR_RNDN); mpfr_div(r22681, r22678, r22680, MPFR_RNDN); ; mpfr_set_si(r22683, mpfr_cmp(r22674, r22682) <= 0, MPFR_RNDN); ; mpfr_set_d(r22685, a, MPFR_RNDN); mpfr_mul(r22686, r22684, r22685, MPFR_RNDN); mpfr_neg(r22687, r22674, MPFR_RNDN); mpfr_mul(r22688, r22674, r22674, MPFR_RNDN); ; mpfr_mul(r22690, r22685, r22677, MPFR_RNDN); mpfr_mul(r22691, r22689, r22690, MPFR_RNDN); mpfr_sub(r22692, r22688, r22691, MPFR_RNDN); mpfr_sqrt(r22693, r22692, MPFR_RNDN); mpfr_sub(r22694, r22687, r22693, MPFR_RNDN); mpfr_div(r22695, r22686, r22694, MPFR_RNDN); mpfr_div(r22696, r22679, r22695, MPFR_RNDN); mpfr_div(r22697, r22687, r22685, MPFR_RNDN); if (mpfr_get_si(r22683, MPFR_RNDN)) { mpfr_set(r22698, r22696, MPFR_RNDN); } else { mpfr_set(r22698, r22697, MPFR_RNDN); }; if (mpfr_get_si(r22676, MPFR_RNDN)) { mpfr_set(r22699, r22681, MPFR_RNDN); } else { mpfr_set(r22699, r22698, MPFR_RNDN); }; return mpfr_get_d(r22699, MPFR_RNDN); } static mpfr_t r22700, r22701, r22702, r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725; void setup_mpfr_f_dm() { mpfr_set_default_prec(3216); mpfr_init(r22700); mpfr_init_set_str(r22701, "-1.3775702482140898e-117", 10, MPFR_RNDN); mpfr_init(r22702); mpfr_init(r22703); mpfr_init(r22704); mpfr_init_set_str(r22705, "1", 10, MPFR_RNDN); mpfr_init(r22706); mpfr_init(r22707); mpfr_init_set_str(r22708, "3.2254415904266855e+123", 10, MPFR_RNDN); mpfr_init(r22709); mpfr_init_set_str(r22710, "2", 10, MPFR_RNDN); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); mpfr_init_set_str(r22715, "4", 10, MPFR_RNDN); mpfr_init(r22716); mpfr_init(r22717); mpfr_init(r22718); mpfr_init(r22719); mpfr_init(r22720); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); mpfr_init(r22725); } double f_dm(double a, double b, double c) { mpfr_set_d(r22700, b, MPFR_RNDN); ; mpfr_set_si(r22702, mpfr_cmp(r22700, r22701) <= 0, MPFR_RNDN); mpfr_set_d(r22703, c, MPFR_RNDN); mpfr_neg(r22704, r22703, MPFR_RNDN); ; mpfr_div(r22706, r22700, r22705, MPFR_RNDN); mpfr_div(r22707, r22704, r22706, MPFR_RNDN); ; mpfr_set_si(r22709, mpfr_cmp(r22700, r22708) <= 0, MPFR_RNDN); ; mpfr_set_d(r22711, a, MPFR_RNDN); mpfr_mul(r22712, r22710, r22711, MPFR_RNDN); mpfr_neg(r22713, r22700, MPFR_RNDN); mpfr_mul(r22714, r22700, r22700, MPFR_RNDN); ; mpfr_mul(r22716, r22711, r22703, MPFR_RNDN); mpfr_mul(r22717, r22715, r22716, MPFR_RNDN); mpfr_sub(r22718, r22714, r22717, MPFR_RNDN); mpfr_sqrt(r22719, r22718, MPFR_RNDN); mpfr_sub(r22720, r22713, r22719, MPFR_RNDN); mpfr_div(r22721, r22712, r22720, MPFR_RNDN); mpfr_div(r22722, r22705, r22721, MPFR_RNDN); mpfr_div(r22723, r22713, r22711, MPFR_RNDN); if (mpfr_get_si(r22709, MPFR_RNDN)) { mpfr_set(r22724, r22722, MPFR_RNDN); } else { mpfr_set(r22724, r22723, MPFR_RNDN); }; if (mpfr_get_si(r22702, MPFR_RNDN)) { mpfr_set(r22725, r22707, MPFR_RNDN); } else { mpfr_set(r22725, r22724, MPFR_RNDN); }; return mpfr_get_d(r22725, MPFR_RNDN); }