#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r10612 = 2.0f; float r10613 = atan2(1.0, 0.0); float r10614 = r10612 * r10613; float r10615 = 3.0f; float r10616 = r10614 / r10615; float r10617 = g; float r10618 = -r10617; float r10619 = h; float r10620 = r10618 / r10619; float r10621 = acos(r10620); float r10622 = r10621 / r10615; float r10623 = r10616 + r10622; float r10624 = cos(r10623); float r10625 = r10612 * r10624; return r10625; } double f_id(double g, double h) { double r10626 = 2.0; double r10627 = atan2(1.0, 0.0); double r10628 = r10626 * r10627; double r10629 = 3.0; double r10630 = r10628 / r10629; double r10631 = g; double r10632 = -r10631; double r10633 = h; double r10634 = r10632 / r10633; double r10635 = acos(r10634); double r10636 = r10635 / r10629; double r10637 = r10630 + r10636; double r10638 = cos(r10637); double r10639 = r10626 * r10638; return r10639; } double f_of(float g, float h) { float r10640 = 1.0f; float r10641 = 2.0f; float r10642 = r10640 + r10641; float r10643 = atan2(1.0, 0.0); float r10644 = 3.0f; float r10645 = r10643 / r10644; float r10646 = r10642 * r10645; float r10647 = cos(r10646); float r10648 = g; float r10649 = h; float r10650 = r10648 / r10649; float r10651 = acos(r10650); float r10652 = r10651 / r10644; float r10653 = cos(r10652); float r10654 = r10647 * r10653; float r10655 = r10643 * r10641; float r10656 = r10655 / r10644; float r10657 = r10656 + r10645; float r10658 = sin(r10657); float r10659 = sin(r10652); float r10660 = r10658 * r10659; float r10661 = r10654 + r10660; float r10662 = r10661 * r10641; return r10662; } double f_od(double g, double h) { double r10663 = 1.0; double r10664 = 2.0; double r10665 = r10663 + r10664; double r10666 = atan2(1.0, 0.0); double r10667 = 3.0; double r10668 = r10666 / r10667; double r10669 = r10665 * r10668; double r10670 = cos(r10669); double r10671 = g; double r10672 = h; double r10673 = r10671 / r10672; double r10674 = acos(r10673); double r10675 = r10674 / r10667; double r10676 = cos(r10675); double r10677 = r10670 * r10676; double r10678 = r10666 * r10664; double r10679 = r10678 / r10667; double r10680 = r10679 + r10668; double r10681 = sin(r10680); double r10682 = sin(r10675); double r10683 = r10681 * r10682; double r10684 = r10677 + r10683; double r10685 = r10684 * r10664; return r10685; } 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 r10686, r10687, r10688, r10689, r10690, r10691, r10692, r10693, r10694, r10695, r10696, r10697, r10698, r10699; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r10686, "2", 10, MPFR_RNDN); mpfr_init(r10687); mpfr_init(r10688); mpfr_init_set_str(r10689, "3", 10, MPFR_RNDN); mpfr_init(r10690); mpfr_init(r10691); mpfr_init(r10692); mpfr_init(r10693); mpfr_init(r10694); mpfr_init(r10695); mpfr_init(r10696); mpfr_init(r10697); mpfr_init(r10698); mpfr_init(r10699); } double f_im(double g, double h) { ; mpfr_const_pi(r10687, MPFR_RNDN); mpfr_mul(r10688, r10686, r10687, MPFR_RNDN); ; mpfr_div(r10690, r10688, r10689, MPFR_RNDN); mpfr_set_d(r10691, g, MPFR_RNDN); mpfr_neg(r10692, r10691, MPFR_RNDN); mpfr_set_d(r10693, h, MPFR_RNDN); mpfr_div(r10694, r10692, r10693, MPFR_RNDN); mpfr_acos(r10695, r10694, MPFR_RNDN); mpfr_div(r10696, r10695, r10689, MPFR_RNDN); mpfr_add(r10697, r10690, r10696, MPFR_RNDN); mpfr_cos(r10698, r10697, MPFR_RNDN); mpfr_mul(r10699, r10686, r10698, MPFR_RNDN); return mpfr_get_d(r10699, MPFR_RNDN); } static mpfr_t r10700, r10701, r10702, r10703, r10704, r10705, r10706, r10707, r10708, r10709, r10710, r10711, r10712, r10713, r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721, r10722; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r10700, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10701, "2", 10, MPFR_RNDN); mpfr_init(r10702); mpfr_init(r10703); mpfr_init_set_str(r10704, "3", 10, MPFR_RNDN); mpfr_init(r10705); mpfr_init(r10706); mpfr_init(r10707); mpfr_init(r10708); mpfr_init(r10709); mpfr_init(r10710); mpfr_init(r10711); mpfr_init(r10712); mpfr_init(r10713); mpfr_init(r10714); mpfr_init(r10715); mpfr_init(r10716); mpfr_init(r10717); mpfr_init(r10718); mpfr_init(r10719); mpfr_init(r10720); mpfr_init(r10721); mpfr_init(r10722); } double f_fm(double g, double h) { ; ; mpfr_add(r10702, r10700, r10701, MPFR_RNDN); mpfr_const_pi(r10703, MPFR_RNDN); ; mpfr_div(r10705, r10703, r10704, MPFR_RNDN); mpfr_mul(r10706, r10702, r10705, MPFR_RNDN); mpfr_cos(r10707, r10706, MPFR_RNDN); mpfr_set_d(r10708, g, MPFR_RNDN); mpfr_set_d(r10709, h, MPFR_RNDN); mpfr_div(r10710, r10708, r10709, MPFR_RNDN); mpfr_acos(r10711, r10710, MPFR_RNDN); mpfr_div(r10712, r10711, r10704, MPFR_RNDN); mpfr_cos(r10713, r10712, MPFR_RNDN); mpfr_mul(r10714, r10707, r10713, MPFR_RNDN); mpfr_mul(r10715, r10703, r10701, MPFR_RNDN); mpfr_div(r10716, r10715, r10704, MPFR_RNDN); mpfr_add(r10717, r10716, r10705, MPFR_RNDN); mpfr_sin(r10718, r10717, MPFR_RNDN); mpfr_sin(r10719, r10712, MPFR_RNDN); mpfr_mul(r10720, r10718, r10719, MPFR_RNDN); mpfr_add(r10721, r10714, r10720, MPFR_RNDN); mpfr_mul(r10722, r10721, r10701, MPFR_RNDN); return mpfr_get_d(r10722, MPFR_RNDN); } static mpfr_t r10723, r10724, r10725, r10726, r10727, r10728, r10729, r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739, r10740, r10741, r10742, r10743, r10744, r10745; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r10723, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10724, "2", 10, MPFR_RNDN); mpfr_init(r10725); mpfr_init(r10726); mpfr_init_set_str(r10727, "3", 10, MPFR_RNDN); 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); mpfr_init(r10740); mpfr_init(r10741); mpfr_init(r10742); mpfr_init(r10743); mpfr_init(r10744); mpfr_init(r10745); } double f_dm(double g, double h) { ; ; mpfr_add(r10725, r10723, r10724, MPFR_RNDN); mpfr_const_pi(r10726, MPFR_RNDN); ; mpfr_div(r10728, r10726, r10727, MPFR_RNDN); mpfr_mul(r10729, r10725, r10728, MPFR_RNDN); mpfr_cos(r10730, r10729, MPFR_RNDN); mpfr_set_d(r10731, g, MPFR_RNDN); mpfr_set_d(r10732, h, MPFR_RNDN); mpfr_div(r10733, r10731, r10732, MPFR_RNDN); mpfr_acos(r10734, r10733, MPFR_RNDN); mpfr_div(r10735, r10734, r10727, MPFR_RNDN); mpfr_cos(r10736, r10735, MPFR_RNDN); mpfr_mul(r10737, r10730, r10736, MPFR_RNDN); mpfr_mul(r10738, r10726, r10724, MPFR_RNDN); mpfr_div(r10739, r10738, r10727, MPFR_RNDN); mpfr_add(r10740, r10739, r10728, MPFR_RNDN); mpfr_sin(r10741, r10740, MPFR_RNDN); mpfr_sin(r10742, r10735, MPFR_RNDN); mpfr_mul(r10743, r10741, r10742, MPFR_RNDN); mpfr_add(r10744, r10737, r10743, MPFR_RNDN); mpfr_mul(r10745, r10744, r10724, MPFR_RNDN); return mpfr_get_d(r10745, MPFR_RNDN); }