#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r10503 = 1.0f; float r10504 = 2.0f; float r10505 = a; float r10506 = r10504 * r10505; float r10507 = r10503 / r10506; float r10508 = g; float r10509 = -r10508; float r10510 = r10508 * r10508; float r10511 = h; float r10512 = r10511 * r10511; float r10513 = r10510 - r10512; float r10514 = sqrt(r10513); float r10515 = r10509 + r10514; float r10516 = r10507 * r10515; float r10517 = cbrt(r10516); float r10518 = r10509 - r10514; float r10519 = r10507 * r10518; float r10520 = cbrt(r10519); float r10521 = r10517 + r10520; return r10521; } double f_id(double g, double h, double a) { double r10522 = 1.0; double r10523 = 2.0; double r10524 = a; double r10525 = r10523 * r10524; double r10526 = r10522 / r10525; double r10527 = g; double r10528 = -r10527; double r10529 = r10527 * r10527; double r10530 = h; double r10531 = r10530 * r10530; double r10532 = r10529 - r10531; double r10533 = sqrt(r10532); double r10534 = r10528 + r10533; double r10535 = r10526 * r10534; double r10536 = cbrt(r10535); double r10537 = r10528 - r10533; double r10538 = r10526 * r10537; double r10539 = cbrt(r10538); double r10540 = r10536 + r10539; return r10540; } double f_of(float g, float h, float a) { float r10541 = g; float r10542 = -r10541; float r10543 = h; float r10544 = r10541 - r10543; float r10545 = r10541 + r10543; float r10546 = r10544 * r10545; float r10547 = sqrt(r10546); float r10548 = r10542 - r10547; float r10549 = cbrt(r10548); float r10550 = a; float r10551 = 2.0f; float r10552 = r10550 * r10551; float r10553 = cbrt(r10552); float r10554 = r10549 / r10553; float r10555 = r10542 + r10541; float r10556 = cbrt(r10555); float r10557 = 1.0f; float r10558 = r10557 / r10552; float r10559 = cbrt(r10558); float r10560 = r10556 * r10559; float r10561 = r10554 + r10560; float r10562 = -4.644802488352878e-113f; bool r10563 = r10561 <= r10562; float r10564 = 3.1391489186073185e-117f; bool r10565 = r10561 <= r10564; float r10566 = r10541 * r10541; float r10567 = r10543 * r10543; float r10568 = r10566 - r10567; float r10569 = sqrt(r10568); float r10570 = r10542 - r10569; float r10571 = r10570 * r10558; float r10572 = cbrt(r10571); float r10573 = r10541 + r10541; float r10574 = -r10573; float r10575 = cbrt(r10574); float r10576 = r10575 * r10559; float r10577 = r10572 + r10576; float r10578 = cbrt(r10570); float r10579 = r10559 * r10578; float r10580 = r10569 + r10542; float r10581 = cbrt(r10580); float r10582 = r10559 * r10581; float r10583 = r10579 + r10582; float r10584 = r10565 ? r10577 : r10583; float r10585 = r10563 ? r10561 : r10584; return r10585; } double f_od(double g, double h, double a) { double r10586 = g; double r10587 = -r10586; double r10588 = h; double r10589 = r10586 - r10588; double r10590 = r10586 + r10588; double r10591 = r10589 * r10590; double r10592 = sqrt(r10591); double r10593 = r10587 - r10592; double r10594 = cbrt(r10593); double r10595 = a; double r10596 = 2.0; double r10597 = r10595 * r10596; double r10598 = cbrt(r10597); double r10599 = r10594 / r10598; double r10600 = r10587 + r10586; double r10601 = cbrt(r10600); double r10602 = 1.0; double r10603 = r10602 / r10597; double r10604 = cbrt(r10603); double r10605 = r10601 * r10604; double r10606 = r10599 + r10605; double r10607 = -4.644802488352878e-113; bool r10608 = r10606 <= r10607; double r10609 = 3.1391489186073185e-117; bool r10610 = r10606 <= r10609; double r10611 = r10586 * r10586; double r10612 = r10588 * r10588; double r10613 = r10611 - r10612; double r10614 = sqrt(r10613); double r10615 = r10587 - r10614; double r10616 = r10615 * r10603; double r10617 = cbrt(r10616); double r10618 = r10586 + r10586; double r10619 = -r10618; double r10620 = cbrt(r10619); double r10621 = r10620 * r10604; double r10622 = r10617 + r10621; double r10623 = cbrt(r10615); double r10624 = r10604 * r10623; double r10625 = r10614 + r10587; double r10626 = cbrt(r10625); double r10627 = r10604 * r10626; double r10628 = r10624 + r10627; double r10629 = r10610 ? r10622 : r10628; double r10630 = r10608 ? r10606 : r10629; return r10630; } 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 r10631, r10632, r10633, r10634, r10635, r10636, r10637, r10638, r10639, r10640, r10641, r10642, r10643, r10644, r10645, r10646, r10647, r10648, r10649; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r10631, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10632, "2", 10, MPFR_RNDN); mpfr_init(r10633); mpfr_init(r10634); mpfr_init(r10635); mpfr_init(r10636); mpfr_init(r10637); mpfr_init(r10638); mpfr_init(r10639); mpfr_init(r10640); mpfr_init(r10641); mpfr_init(r10642); mpfr_init(r10643); mpfr_init(r10644); mpfr_init(r10645); mpfr_init(r10646); mpfr_init(r10647); mpfr_init(r10648); mpfr_init(r10649); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r10633, a, MPFR_RNDN); mpfr_mul(r10634, r10632, r10633, MPFR_RNDN); mpfr_div(r10635, r10631, r10634, MPFR_RNDN); mpfr_set_d(r10636, g, MPFR_RNDN); mpfr_neg(r10637, r10636, MPFR_RNDN); mpfr_mul(r10638, r10636, r10636, MPFR_RNDN); mpfr_set_d(r10639, h, MPFR_RNDN); mpfr_mul(r10640, r10639, r10639, MPFR_RNDN); mpfr_sub(r10641, r10638, r10640, MPFR_RNDN); mpfr_sqrt(r10642, r10641, MPFR_RNDN); mpfr_add(r10643, r10637, r10642, MPFR_RNDN); mpfr_mul(r10644, r10635, r10643, MPFR_RNDN); mpfr_cbrt(r10645, r10644, MPFR_RNDN); mpfr_sub(r10646, r10637, r10642, MPFR_RNDN); mpfr_mul(r10647, r10635, r10646, MPFR_RNDN); mpfr_cbrt(r10648, r10647, MPFR_RNDN); mpfr_add(r10649, r10645, r10648, MPFR_RNDN); return mpfr_get_d(r10649, MPFR_RNDN); } static mpfr_t r10650, r10651, r10652, r10653, r10654, r10655, r10656, r10657, r10658, r10659, r10660, r10661, r10662, r10663, r10664, r10665, r10666, r10667, r10668, r10669, r10670, r10671, r10672, r10673, r10674, r10675, r10676, r10677, r10678, r10679, r10680, r10681, r10682, r10683, r10684, r10685, r10686, r10687, r10688, r10689, r10690, r10691, r10692, r10693, r10694; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10650); mpfr_init(r10651); mpfr_init(r10652); mpfr_init(r10653); mpfr_init(r10654); mpfr_init(r10655); mpfr_init(r10656); mpfr_init(r10657); mpfr_init(r10658); mpfr_init(r10659); mpfr_init_set_str(r10660, "2", 10, MPFR_RNDN); mpfr_init(r10661); mpfr_init(r10662); mpfr_init(r10663); mpfr_init(r10664); mpfr_init(r10665); mpfr_init_set_str(r10666, "1", 10, MPFR_RNDN); mpfr_init(r10667); mpfr_init(r10668); mpfr_init(r10669); mpfr_init(r10670); mpfr_init_set_str(r10671, "-4.644802488352878e-113", 10, MPFR_RNDN); mpfr_init(r10672); mpfr_init_set_str(r10673, "3.1391489186073185e-117", 10, MPFR_RNDN); mpfr_init(r10674); mpfr_init(r10675); mpfr_init(r10676); mpfr_init(r10677); mpfr_init(r10678); mpfr_init(r10679); mpfr_init(r10680); mpfr_init(r10681); mpfr_init(r10682); mpfr_init(r10683); mpfr_init(r10684); mpfr_init(r10685); mpfr_init(r10686); mpfr_init(r10687); mpfr_init(r10688); mpfr_init(r10689); mpfr_init(r10690); mpfr_init(r10691); mpfr_init(r10692); mpfr_init(r10693); mpfr_init(r10694); } double f_fm(double g, double h, double a) { mpfr_set_d(r10650, g, MPFR_RNDN); mpfr_neg(r10651, r10650, MPFR_RNDN); mpfr_set_d(r10652, h, MPFR_RNDN); mpfr_sub(r10653, r10650, r10652, MPFR_RNDN); mpfr_add(r10654, r10650, r10652, MPFR_RNDN); mpfr_mul(r10655, r10653, r10654, MPFR_RNDN); mpfr_sqrt(r10656, r10655, MPFR_RNDN); mpfr_sub(r10657, r10651, r10656, MPFR_RNDN); mpfr_cbrt(r10658, r10657, MPFR_RNDN); mpfr_set_d(r10659, a, MPFR_RNDN); ; mpfr_mul(r10661, r10659, r10660, MPFR_RNDN); mpfr_cbrt(r10662, r10661, MPFR_RNDN); mpfr_div(r10663, r10658, r10662, MPFR_RNDN); mpfr_add(r10664, r10651, r10650, MPFR_RNDN); mpfr_cbrt(r10665, r10664, MPFR_RNDN); ; mpfr_div(r10667, r10666, r10661, MPFR_RNDN); mpfr_cbrt(r10668, r10667, MPFR_RNDN); mpfr_mul(r10669, r10665, r10668, MPFR_RNDN); mpfr_add(r10670, r10663, r10669, MPFR_RNDN); ; mpfr_set_si(r10672, mpfr_cmp(r10670, r10671) <= 0, MPFR_RNDN); ; mpfr_set_si(r10674, mpfr_cmp(r10670, r10673) <= 0, MPFR_RNDN); mpfr_mul(r10675, r10650, r10650, MPFR_RNDN); mpfr_mul(r10676, r10652, r10652, MPFR_RNDN); mpfr_sub(r10677, r10675, r10676, MPFR_RNDN); mpfr_sqrt(r10678, r10677, MPFR_RNDN); mpfr_sub(r10679, r10651, r10678, MPFR_RNDN); mpfr_mul(r10680, r10679, r10667, MPFR_RNDN); mpfr_cbrt(r10681, r10680, MPFR_RNDN); mpfr_add(r10682, r10650, r10650, MPFR_RNDN); mpfr_neg(r10683, r10682, MPFR_RNDN); mpfr_cbrt(r10684, r10683, MPFR_RNDN); mpfr_mul(r10685, r10684, r10668, MPFR_RNDN); mpfr_add(r10686, r10681, r10685, MPFR_RNDN); mpfr_cbrt(r10687, r10679, MPFR_RNDN); mpfr_mul(r10688, r10668, r10687, MPFR_RNDN); mpfr_add(r10689, r10678, r10651, MPFR_RNDN); mpfr_cbrt(r10690, r10689, MPFR_RNDN); mpfr_mul(r10691, r10668, r10690, MPFR_RNDN); mpfr_add(r10692, r10688, r10691, MPFR_RNDN); if (mpfr_get_si(r10674, MPFR_RNDN)) { mpfr_set(r10693, r10686, MPFR_RNDN); } else { mpfr_set(r10693, r10692, MPFR_RNDN); }; if (mpfr_get_si(r10672, MPFR_RNDN)) { mpfr_set(r10694, r10670, MPFR_RNDN); } else { mpfr_set(r10694, r10693, MPFR_RNDN); }; return mpfr_get_d(r10694, MPFR_RNDN); } static mpfr_t r10695, r10696, r10697, r10698, r10699, r10700, r10701, r10702, r10703, r10704, r10705, r10706, r10707, r10708, r10709, r10710, r10711, r10712, r10713, r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721, r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729, r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10695); mpfr_init(r10696); mpfr_init(r10697); mpfr_init(r10698); mpfr_init(r10699); mpfr_init(r10700); mpfr_init(r10701); mpfr_init(r10702); mpfr_init(r10703); mpfr_init(r10704); mpfr_init_set_str(r10705, "2", 10, MPFR_RNDN); mpfr_init(r10706); mpfr_init(r10707); mpfr_init(r10708); mpfr_init(r10709); mpfr_init(r10710); mpfr_init_set_str(r10711, "1", 10, MPFR_RNDN); mpfr_init(r10712); mpfr_init(r10713); mpfr_init(r10714); mpfr_init(r10715); mpfr_init_set_str(r10716, "-4.644802488352878e-113", 10, MPFR_RNDN); mpfr_init(r10717); mpfr_init_set_str(r10718, "3.1391489186073185e-117", 10, MPFR_RNDN); mpfr_init(r10719); mpfr_init(r10720); mpfr_init(r10721); mpfr_init(r10722); mpfr_init(r10723); mpfr_init(r10724); mpfr_init(r10725); mpfr_init(r10726); mpfr_init(r10727); mpfr_init(r10728); mpfr_init(r10729); mpfr_init(r10730); mpfr_init(r10731); mpfr_init(r10732); mpfr_init(r10733); mpfr_init(r10734); mpfr_init(r10735); mpfr_init(r10736); mpfr_init(r10737); mpfr_init(r10738); mpfr_init(r10739); } double f_dm(double g, double h, double a) { mpfr_set_d(r10695, g, MPFR_RNDN); mpfr_neg(r10696, r10695, MPFR_RNDN); mpfr_set_d(r10697, h, MPFR_RNDN); mpfr_sub(r10698, r10695, r10697, MPFR_RNDN); mpfr_add(r10699, r10695, r10697, MPFR_RNDN); mpfr_mul(r10700, r10698, r10699, MPFR_RNDN); mpfr_sqrt(r10701, r10700, MPFR_RNDN); mpfr_sub(r10702, r10696, r10701, MPFR_RNDN); mpfr_cbrt(r10703, r10702, MPFR_RNDN); mpfr_set_d(r10704, a, MPFR_RNDN); ; mpfr_mul(r10706, r10704, r10705, MPFR_RNDN); mpfr_cbrt(r10707, r10706, MPFR_RNDN); mpfr_div(r10708, r10703, r10707, MPFR_RNDN); mpfr_add(r10709, r10696, r10695, MPFR_RNDN); mpfr_cbrt(r10710, r10709, MPFR_RNDN); ; mpfr_div(r10712, r10711, r10706, MPFR_RNDN); mpfr_cbrt(r10713, r10712, MPFR_RNDN); mpfr_mul(r10714, r10710, r10713, MPFR_RNDN); mpfr_add(r10715, r10708, r10714, MPFR_RNDN); ; mpfr_set_si(r10717, mpfr_cmp(r10715, r10716) <= 0, MPFR_RNDN); ; mpfr_set_si(r10719, mpfr_cmp(r10715, r10718) <= 0, MPFR_RNDN); mpfr_mul(r10720, r10695, r10695, MPFR_RNDN); mpfr_mul(r10721, r10697, r10697, MPFR_RNDN); mpfr_sub(r10722, r10720, r10721, MPFR_RNDN); mpfr_sqrt(r10723, r10722, MPFR_RNDN); mpfr_sub(r10724, r10696, r10723, MPFR_RNDN); mpfr_mul(r10725, r10724, r10712, MPFR_RNDN); mpfr_cbrt(r10726, r10725, MPFR_RNDN); mpfr_add(r10727, r10695, r10695, MPFR_RNDN); mpfr_neg(r10728, r10727, MPFR_RNDN); mpfr_cbrt(r10729, r10728, MPFR_RNDN); mpfr_mul(r10730, r10729, r10713, MPFR_RNDN); mpfr_add(r10731, r10726, r10730, MPFR_RNDN); mpfr_cbrt(r10732, r10724, MPFR_RNDN); mpfr_mul(r10733, r10713, r10732, MPFR_RNDN); mpfr_add(r10734, r10723, r10696, MPFR_RNDN); mpfr_cbrt(r10735, r10734, MPFR_RNDN); mpfr_mul(r10736, r10713, r10735, MPFR_RNDN); mpfr_add(r10737, r10733, r10736, MPFR_RNDN); if (mpfr_get_si(r10719, MPFR_RNDN)) { mpfr_set(r10738, r10731, MPFR_RNDN); } else { mpfr_set(r10738, r10737, MPFR_RNDN); }; if (mpfr_get_si(r10717, MPFR_RNDN)) { mpfr_set(r10739, r10715, MPFR_RNDN); } else { mpfr_set(r10739, r10738, MPFR_RNDN); }; return mpfr_get_d(r10739, MPFR_RNDN); }