#include #include #include #include #include char *name = "The quadratic formula (r2)"; double f_if(float a, float b, float c) { float r25612 = b; float r25613 = -r25612; float r25614 = r25612 * r25612; float r25615 = 4; float r25616 = a; float r25617 = c; float r25618 = r25616 * r25617; float r25619 = r25615 * r25618; float r25620 = r25614 - r25619; float r25621 = sqrt(r25620); float r25622 = r25613 - r25621; float r25623 = 2; float r25624 = r25623 * r25616; float r25625 = r25622 / r25624; return r25625; } double f_id(double a, double b, double c) { double r25626 = b; double r25627 = -r25626; double r25628 = r25626 * r25626; double r25629 = 4; double r25630 = a; double r25631 = c; double r25632 = r25630 * r25631; double r25633 = r25629 * r25632; double r25634 = r25628 - r25633; double r25635 = sqrt(r25634); double r25636 = r25627 - r25635; double r25637 = 2; double r25638 = r25637 * r25630; double r25639 = r25636 / r25638; return r25639; } double f_of(float a, float b, float c) { float r25640 = b; float r25641 = -6.606286567414996e-28; bool r25642 = r25640 <= r25641; float r25643 = c; float r25644 = -r25643; float r25645 = r25644 / r25640; float r25646 = 9.633104463257503e+103; bool r25647 = r25640 <= r25646; float r25648 = 1; float r25649 = a; float r25650 = 2; float r25651 = r25649 * r25650; float r25652 = -r25640; float r25653 = r25640 * r25640; float r25654 = r25649 * r25643; float r25655 = 4; float r25656 = r25654 * r25655; float r25657 = r25653 - r25656; float r25658 = sqrt(r25657); float r25659 = r25652 - r25658; float r25660 = r25651 / r25659; float r25661 = r25648 / r25660; float r25662 = r25640 / r25649; float r25663 = -r25662; float r25664 = r25647 ? r25661 : r25663; float r25665 = r25642 ? r25645 : r25664; return r25665; } double f_od(double a, double b, double c) { double r25666 = b; double r25667 = -6.606286567414996e-28; bool r25668 = r25666 <= r25667; double r25669 = c; double r25670 = -r25669; double r25671 = r25670 / r25666; double r25672 = 9.633104463257503e+103; bool r25673 = r25666 <= r25672; double r25674 = 1; double r25675 = a; double r25676 = 2; double r25677 = r25675 * r25676; double r25678 = -r25666; double r25679 = r25666 * r25666; double r25680 = r25675 * r25669; double r25681 = 4; double r25682 = r25680 * r25681; double r25683 = r25679 - r25682; double r25684 = sqrt(r25683); double r25685 = r25678 - r25684; double r25686 = r25677 / r25685; double r25687 = r25674 / r25686; double r25688 = r25666 / r25675; double r25689 = -r25688; double r25690 = r25673 ? r25687 : r25689; double r25691 = r25668 ? r25671 : r25690; return r25691; } 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 r25692, r25693, r25694, r25695, r25696, r25697, r25698, r25699, r25700, r25701, r25702, r25703, r25704, r25705; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init(r25692); mpfr_init(r25693); mpfr_init(r25694); mpfr_init_set_str(r25695, "4", 10, MPFR_RNDN); mpfr_init(r25696); mpfr_init(r25697); mpfr_init(r25698); mpfr_init(r25699); mpfr_init(r25700); mpfr_init(r25701); mpfr_init(r25702); mpfr_init_set_str(r25703, "2", 10, MPFR_RNDN); mpfr_init(r25704); mpfr_init(r25705); } double f_im(double a, double b, double c) { mpfr_set_d(r25692, b, MPFR_RNDN); mpfr_neg(r25693, r25692, MPFR_RNDN); mpfr_mul(r25694, r25692, r25692, MPFR_RNDN); ; mpfr_set_d(r25696, a, MPFR_RNDN); mpfr_set_d(r25697, c, MPFR_RNDN); mpfr_mul(r25698, r25696, r25697, MPFR_RNDN); mpfr_mul(r25699, r25695, r25698, MPFR_RNDN); mpfr_sub(r25700, r25694, r25699, MPFR_RNDN); mpfr_sqrt(r25701, r25700, MPFR_RNDN); mpfr_sub(r25702, r25693, r25701, MPFR_RNDN); ; mpfr_mul(r25704, r25703, r25696, MPFR_RNDN); mpfr_div(r25705, r25702, r25704, MPFR_RNDN); return mpfr_get_d(r25705, MPFR_RNDN); } static mpfr_t r25706, r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729, r25730, r25731; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init(r25706); mpfr_init_set_str(r25707, "-6.606286567414996e-28", 10, MPFR_RNDN); mpfr_init(r25708); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init_set_str(r25712, "9.633104463257503e+103", 10, MPFR_RNDN); mpfr_init(r25713); mpfr_init_set_str(r25714, "1", 10, MPFR_RNDN); mpfr_init(r25715); mpfr_init_set_str(r25716, "2", 10, MPFR_RNDN); mpfr_init(r25717); mpfr_init(r25718); mpfr_init(r25719); mpfr_init(r25720); mpfr_init_set_str(r25721, "4", 10, MPFR_RNDN); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init(r25727); mpfr_init(r25728); mpfr_init(r25729); mpfr_init(r25730); mpfr_init(r25731); } double f_fm(double a, double b, double c) { mpfr_set_d(r25706, b, MPFR_RNDN); ; mpfr_set_si(r25708, mpfr_cmp(r25706, r25707) <= 0, MPFR_RNDN); mpfr_set_d(r25709, c, MPFR_RNDN); mpfr_neg(r25710, r25709, MPFR_RNDN); mpfr_div(r25711, r25710, r25706, MPFR_RNDN); ; mpfr_set_si(r25713, mpfr_cmp(r25706, r25712) <= 0, MPFR_RNDN); ; mpfr_set_d(r25715, a, MPFR_RNDN); ; mpfr_mul(r25717, r25715, r25716, MPFR_RNDN); mpfr_neg(r25718, r25706, MPFR_RNDN); mpfr_mul(r25719, r25706, r25706, MPFR_RNDN); mpfr_mul(r25720, r25715, r25709, MPFR_RNDN); ; mpfr_mul(r25722, r25720, r25721, MPFR_RNDN); mpfr_sub(r25723, r25719, r25722, MPFR_RNDN); mpfr_sqrt(r25724, r25723, MPFR_RNDN); mpfr_sub(r25725, r25718, r25724, MPFR_RNDN); mpfr_div(r25726, r25717, r25725, MPFR_RNDN); mpfr_div(r25727, r25714, r25726, MPFR_RNDN); mpfr_div(r25728, r25706, r25715, MPFR_RNDN); mpfr_neg(r25729, r25728, MPFR_RNDN); if (mpfr_get_si(r25713, MPFR_RNDN)) { mpfr_set(r25730, r25727, MPFR_RNDN); } else { mpfr_set(r25730, r25729, MPFR_RNDN); }; if (mpfr_get_si(r25708, MPFR_RNDN)) { mpfr_set(r25731, r25711, MPFR_RNDN); } else { mpfr_set(r25731, r25730, MPFR_RNDN); }; return mpfr_get_d(r25731, MPFR_RNDN); } static mpfr_t r25732, r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755, r25756, r25757; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init(r25732); mpfr_init_set_str(r25733, "-6.606286567414996e-28", 10, MPFR_RNDN); mpfr_init(r25734); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); mpfr_init_set_str(r25738, "9.633104463257503e+103", 10, MPFR_RNDN); mpfr_init(r25739); mpfr_init_set_str(r25740, "1", 10, MPFR_RNDN); mpfr_init(r25741); mpfr_init_set_str(r25742, "2", 10, MPFR_RNDN); mpfr_init(r25743); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init_set_str(r25747, "4", 10, MPFR_RNDN); mpfr_init(r25748); mpfr_init(r25749); mpfr_init(r25750); mpfr_init(r25751); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init(r25756); mpfr_init(r25757); } double f_dm(double a, double b, double c) { mpfr_set_d(r25732, b, MPFR_RNDN); ; mpfr_set_si(r25734, mpfr_cmp(r25732, r25733) <= 0, MPFR_RNDN); mpfr_set_d(r25735, c, MPFR_RNDN); mpfr_neg(r25736, r25735, MPFR_RNDN); mpfr_div(r25737, r25736, r25732, MPFR_RNDN); ; mpfr_set_si(r25739, mpfr_cmp(r25732, r25738) <= 0, MPFR_RNDN); ; mpfr_set_d(r25741, a, MPFR_RNDN); ; mpfr_mul(r25743, r25741, r25742, MPFR_RNDN); mpfr_neg(r25744, r25732, MPFR_RNDN); mpfr_mul(r25745, r25732, r25732, MPFR_RNDN); mpfr_mul(r25746, r25741, r25735, MPFR_RNDN); ; mpfr_mul(r25748, r25746, r25747, MPFR_RNDN); mpfr_sub(r25749, r25745, r25748, MPFR_RNDN); mpfr_sqrt(r25750, r25749, MPFR_RNDN); mpfr_sub(r25751, r25744, r25750, MPFR_RNDN); mpfr_div(r25752, r25743, r25751, MPFR_RNDN); mpfr_div(r25753, r25740, r25752, MPFR_RNDN); mpfr_div(r25754, r25732, r25741, MPFR_RNDN); mpfr_neg(r25755, r25754, MPFR_RNDN); if (mpfr_get_si(r25739, MPFR_RNDN)) { mpfr_set(r25756, r25753, MPFR_RNDN); } else { mpfr_set(r25756, r25755, MPFR_RNDN); }; if (mpfr_get_si(r25734, MPFR_RNDN)) { mpfr_set(r25757, r25737, MPFR_RNDN); } else { mpfr_set(r25757, r25756, MPFR_RNDN); }; return mpfr_get_d(r25757, MPFR_RNDN); }