#include #include #include #include #include char *name = "The quadratic formula (r1)"; double f_if(float a, float b, float c) { float r18548 = b; float r18549 = -r18548; float r18550 = r18548 * r18548; float r18551 = 4.0f; float r18552 = a; float r18553 = r18551 * r18552; float r18554 = c; float r18555 = r18553 * r18554; float r18556 = r18550 - r18555; float r18557 = sqrt(r18556); float r18558 = r18549 + r18557; float r18559 = 2.0f; float r18560 = r18559 * r18552; float r18561 = r18558 / r18560; return r18561; } double f_id(double a, double b, double c) { double r18562 = b; double r18563 = -r18562; double r18564 = r18562 * r18562; double r18565 = 4.0; double r18566 = a; double r18567 = r18565 * r18566; double r18568 = c; double r18569 = r18567 * r18568; double r18570 = r18564 - r18569; double r18571 = sqrt(r18570); double r18572 = r18563 + r18571; double r18573 = 2.0; double r18574 = r18573 * r18566; double r18575 = r18572 / r18574; return r18575; } double f_of(float a, float b, float c) { float r18576 = b; float r18577 = -2.2913476789857995e+35f; bool r18578 = r18576 <= r18577; float r18579 = c; float r18580 = r18579 / r18576; float r18581 = a; float r18582 = r18576 / r18581; float r18583 = r18580 - r18582; float r18584 = 2.890760014197913e-73f; bool r18585 = r18576 <= r18584; float r18586 = -r18576; float r18587 = r18576 * r18576; float r18588 = 4.0f; float r18589 = r18588 * r18581; float r18590 = r18589 * r18579; float r18591 = r18587 - r18590; float r18592 = sqrt(r18591); float r18593 = r18586 + r18592; float r18594 = 2.0f; float r18595 = r18594 * r18581; float r18596 = r18593 / r18595; float r18597 = -2.0f; float r18598 = r18597 / r18594; float r18599 = r18580 * r18598; float r18600 = r18585 ? r18596 : r18599; float r18601 = r18578 ? r18583 : r18600; return r18601; } double f_od(double a, double b, double c) { double r18602 = b; double r18603 = -2.2913476789857995e+35; bool r18604 = r18602 <= r18603; double r18605 = c; double r18606 = r18605 / r18602; double r18607 = a; double r18608 = r18602 / r18607; double r18609 = r18606 - r18608; double r18610 = 2.890760014197913e-73; bool r18611 = r18602 <= r18610; double r18612 = -r18602; double r18613 = r18602 * r18602; double r18614 = 4.0; double r18615 = r18614 * r18607; double r18616 = r18615 * r18605; double r18617 = r18613 - r18616; double r18618 = sqrt(r18617); double r18619 = r18612 + r18618; double r18620 = 2.0; double r18621 = r18620 * r18607; double r18622 = r18619 / r18621; double r18623 = -2.0; double r18624 = r18623 / r18620; double r18625 = r18606 * r18624; double r18626 = r18611 ? r18622 : r18625; double r18627 = r18604 ? r18609 : r18626; return r18627; } 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 r18628, r18629, r18630, r18631, r18632, r18633, r18634, r18635, r18636, r18637, r18638, r18639, r18640, r18641; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18628); mpfr_init(r18629); mpfr_init(r18630); mpfr_init_set_str(r18631, "4", 10, MPFR_RNDN); mpfr_init(r18632); mpfr_init(r18633); mpfr_init(r18634); mpfr_init(r18635); mpfr_init(r18636); mpfr_init(r18637); mpfr_init(r18638); mpfr_init_set_str(r18639, "2", 10, MPFR_RNDN); mpfr_init(r18640); mpfr_init(r18641); } double f_im(double a, double b, double c) { mpfr_set_d(r18628, b, MPFR_RNDN); mpfr_neg(r18629, r18628, MPFR_RNDN); mpfr_sqr(r18630, r18628, MPFR_RNDN); ; mpfr_set_d(r18632, a, MPFR_RNDN); mpfr_mul(r18633, r18631, r18632, MPFR_RNDN); mpfr_set_d(r18634, c, MPFR_RNDN); mpfr_mul(r18635, r18633, r18634, MPFR_RNDN); mpfr_sub(r18636, r18630, r18635, MPFR_RNDN); mpfr_sqrt(r18637, r18636, MPFR_RNDN); mpfr_add(r18638, r18629, r18637, MPFR_RNDN); ; mpfr_mul(r18640, r18639, r18632, MPFR_RNDN); mpfr_div(r18641, r18638, r18640, MPFR_RNDN); return mpfr_get_d(r18641, MPFR_RNDN); } static mpfr_t r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18642); mpfr_init_set_str(r18643, "-2.2913476789857995e+35", 10, MPFR_RNDN); mpfr_init(r18644); mpfr_init(r18645); mpfr_init(r18646); mpfr_init(r18647); mpfr_init(r18648); mpfr_init(r18649); mpfr_init_set_str(r18650, "2.890760014197913e-73", 10, MPFR_RNDN); mpfr_init(r18651); mpfr_init(r18652); mpfr_init(r18653); mpfr_init_set_str(r18654, "4", 10, MPFR_RNDN); mpfr_init(r18655); mpfr_init(r18656); mpfr_init(r18657); mpfr_init(r18658); mpfr_init(r18659); mpfr_init_set_str(r18660, "2", 10, MPFR_RNDN); mpfr_init(r18661); mpfr_init(r18662); mpfr_init_set_str(r18663, "-2", 10, MPFR_RNDN); mpfr_init(r18664); mpfr_init(r18665); mpfr_init(r18666); mpfr_init(r18667); } double f_fm(double a, double b, double c) { mpfr_set_d(r18642, b, MPFR_RNDN); ; mpfr_set_si(r18644, mpfr_cmp(r18642, r18643) <= 0, MPFR_RNDN); mpfr_set_d(r18645, c, MPFR_RNDN); mpfr_div(r18646, r18645, r18642, MPFR_RNDN); mpfr_set_d(r18647, a, MPFR_RNDN); mpfr_div(r18648, r18642, r18647, MPFR_RNDN); mpfr_sub(r18649, r18646, r18648, MPFR_RNDN); ; mpfr_set_si(r18651, mpfr_cmp(r18642, r18650) <= 0, MPFR_RNDN); mpfr_neg(r18652, r18642, MPFR_RNDN); mpfr_sqr(r18653, r18642, MPFR_RNDN); ; mpfr_mul(r18655, r18654, r18647, MPFR_RNDN); mpfr_mul(r18656, r18655, r18645, MPFR_RNDN); mpfr_sub(r18657, r18653, r18656, MPFR_RNDN); mpfr_sqrt(r18658, r18657, MPFR_RNDN); mpfr_add(r18659, r18652, r18658, MPFR_RNDN); ; mpfr_mul(r18661, r18660, r18647, MPFR_RNDN); mpfr_div(r18662, r18659, r18661, MPFR_RNDN); ; mpfr_div(r18664, r18663, r18660, MPFR_RNDN); mpfr_mul(r18665, r18646, r18664, MPFR_RNDN); if (mpfr_get_si(r18651, MPFR_RNDN)) { mpfr_set(r18666, r18662, MPFR_RNDN); } else { mpfr_set(r18666, r18665, MPFR_RNDN); }; if (mpfr_get_si(r18644, MPFR_RNDN)) { mpfr_set(r18667, r18649, MPFR_RNDN); } else { mpfr_set(r18667, r18666, MPFR_RNDN); }; return mpfr_get_d(r18667, MPFR_RNDN); } static mpfr_t r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681, r18682, r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18668); mpfr_init_set_str(r18669, "-2.2913476789857995e+35", 10, MPFR_RNDN); mpfr_init(r18670); mpfr_init(r18671); mpfr_init(r18672); mpfr_init(r18673); mpfr_init(r18674); mpfr_init(r18675); mpfr_init_set_str(r18676, "2.890760014197913e-73", 10, MPFR_RNDN); mpfr_init(r18677); mpfr_init(r18678); mpfr_init(r18679); mpfr_init_set_str(r18680, "4", 10, MPFR_RNDN); mpfr_init(r18681); mpfr_init(r18682); mpfr_init(r18683); mpfr_init(r18684); mpfr_init(r18685); mpfr_init_set_str(r18686, "2", 10, MPFR_RNDN); mpfr_init(r18687); mpfr_init(r18688); mpfr_init_set_str(r18689, "-2", 10, MPFR_RNDN); mpfr_init(r18690); mpfr_init(r18691); mpfr_init(r18692); mpfr_init(r18693); } double f_dm(double a, double b, double c) { mpfr_set_d(r18668, b, MPFR_RNDN); ; mpfr_set_si(r18670, mpfr_cmp(r18668, r18669) <= 0, MPFR_RNDN); mpfr_set_d(r18671, c, MPFR_RNDN); mpfr_div(r18672, r18671, r18668, MPFR_RNDN); mpfr_set_d(r18673, a, MPFR_RNDN); mpfr_div(r18674, r18668, r18673, MPFR_RNDN); mpfr_sub(r18675, r18672, r18674, MPFR_RNDN); ; mpfr_set_si(r18677, mpfr_cmp(r18668, r18676) <= 0, MPFR_RNDN); mpfr_neg(r18678, r18668, MPFR_RNDN); mpfr_sqr(r18679, r18668, MPFR_RNDN); ; mpfr_mul(r18681, r18680, r18673, MPFR_RNDN); mpfr_mul(r18682, r18681, r18671, MPFR_RNDN); mpfr_sub(r18683, r18679, r18682, MPFR_RNDN); mpfr_sqrt(r18684, r18683, MPFR_RNDN); mpfr_add(r18685, r18678, r18684, MPFR_RNDN); ; mpfr_mul(r18687, r18686, r18673, MPFR_RNDN); mpfr_div(r18688, r18685, r18687, MPFR_RNDN); ; mpfr_div(r18690, r18689, r18686, MPFR_RNDN); mpfr_mul(r18691, r18672, r18690, MPFR_RNDN); if (mpfr_get_si(r18677, MPFR_RNDN)) { mpfr_set(r18692, r18688, MPFR_RNDN); } else { mpfr_set(r18692, r18691, MPFR_RNDN); }; if (mpfr_get_si(r18670, MPFR_RNDN)) { mpfr_set(r18693, r18675, MPFR_RNDN); } else { mpfr_set(r18693, r18692, MPFR_RNDN); }; return mpfr_get_d(r18693, MPFR_RNDN); }