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