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