#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r23635 = 1; float r23636 = 2; float r23637 = a; float r23638 = r23636 * r23637; float r23639 = r23635 / r23638; float r23640 = g; float r23641 = -r23640; float r23642 = r23640 * r23640; float r23643 = h; float r23644 = r23643 * r23643; float r23645 = r23642 - r23644; float r23646 = sqrt(r23645); float r23647 = r23641 + r23646; float r23648 = r23639 * r23647; float r23649 = cbrt(r23648); float r23650 = r23641 - r23646; float r23651 = r23639 * r23650; float r23652 = cbrt(r23651); float r23653 = r23649 + r23652; return r23653; } double f_id(double g, double h, double a) { double r23654 = 1; double r23655 = 2; double r23656 = a; double r23657 = r23655 * r23656; double r23658 = r23654 / r23657; double r23659 = g; double r23660 = -r23659; double r23661 = r23659 * r23659; double r23662 = h; double r23663 = r23662 * r23662; double r23664 = r23661 - r23663; double r23665 = sqrt(r23664); double r23666 = r23660 + r23665; double r23667 = r23658 * r23666; double r23668 = cbrt(r23667); double r23669 = r23660 - r23665; double r23670 = r23658 * r23669; double r23671 = cbrt(r23670); double r23672 = r23668 + r23671; return r23672; } double f_of(float g, float h, float a) { float r23673 = g; float r23674 = -1.6784217127088642e-162; bool r23675 = r23673 <= r23674; float r23676 = 1; float r23677 = 2; float r23678 = a; float r23679 = r23677 * r23678; float r23680 = r23676 / r23679; float r23681 = -r23673; float r23682 = r23673 * r23673; float r23683 = h; float r23684 = r23683 * r23683; float r23685 = r23682 - r23684; float r23686 = sqrt(r23685); float r23687 = r23681 + r23686; float r23688 = r23680 * r23687; float r23689 = cbrt(r23688); float r23690 = r23684 / r23687; float r23691 = r23680 * r23690; float r23692 = cbrt(r23691); float r23693 = r23689 + r23692; float r23694 = r23681 - r23686; float r23695 = r23684 / r23694; float r23696 = r23680 * r23695; float r23697 = cbrt(r23696); float r23698 = r23680 * r23694; float r23699 = cbrt(r23698); float r23700 = r23697 + r23699; float r23701 = r23675 ? r23693 : r23700; return r23701; } double f_od(double g, double h, double a) { double r23702 = g; double r23703 = -1.6784217127088642e-162; bool r23704 = r23702 <= r23703; double r23705 = 1; double r23706 = 2; double r23707 = a; double r23708 = r23706 * r23707; double r23709 = r23705 / r23708; double r23710 = -r23702; double r23711 = r23702 * r23702; double r23712 = h; double r23713 = r23712 * r23712; double r23714 = r23711 - r23713; double r23715 = sqrt(r23714); double r23716 = r23710 + r23715; double r23717 = r23709 * r23716; double r23718 = cbrt(r23717); double r23719 = r23713 / r23716; double r23720 = r23709 * r23719; double r23721 = cbrt(r23720); double r23722 = r23718 + r23721; double r23723 = r23710 - r23715; double r23724 = r23713 / r23723; double r23725 = r23709 * r23724; double r23726 = cbrt(r23725); double r23727 = r23709 * r23723; double r23728 = cbrt(r23727); double r23729 = r23726 + r23728; double r23730 = r23704 ? r23722 : r23729; return r23730; } 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 r23731, r23732, r23733, r23734, r23735, r23736, r23737, r23738, r23739, r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747, r23748, r23749; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init_set_str(r23731, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23732, "2", 10, MPFR_RNDN); mpfr_init(r23733); mpfr_init(r23734); mpfr_init(r23735); mpfr_init(r23736); mpfr_init(r23737); mpfr_init(r23738); mpfr_init(r23739); mpfr_init(r23740); mpfr_init(r23741); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init(r23745); mpfr_init(r23746); mpfr_init(r23747); mpfr_init(r23748); mpfr_init(r23749); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r23733, a, MPFR_RNDN); mpfr_mul(r23734, r23732, r23733, MPFR_RNDN); mpfr_div(r23735, r23731, r23734, MPFR_RNDN); mpfr_set_d(r23736, g, MPFR_RNDN); mpfr_neg(r23737, r23736, MPFR_RNDN); mpfr_mul(r23738, r23736, r23736, MPFR_RNDN); mpfr_set_d(r23739, h, MPFR_RNDN); mpfr_mul(r23740, r23739, r23739, MPFR_RNDN); mpfr_sub(r23741, r23738, r23740, MPFR_RNDN); mpfr_sqrt(r23742, r23741, MPFR_RNDN); mpfr_add(r23743, r23737, r23742, MPFR_RNDN); mpfr_mul(r23744, r23735, r23743, MPFR_RNDN); mpfr_cbrt(r23745, r23744, MPFR_RNDN); mpfr_sub(r23746, r23737, r23742, MPFR_RNDN); mpfr_mul(r23747, r23735, r23746, MPFR_RNDN); mpfr_cbrt(r23748, r23747, MPFR_RNDN); mpfr_add(r23749, r23745, r23748, MPFR_RNDN); return mpfr_get_d(r23749, MPFR_RNDN); } static mpfr_t r23750, r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23750); mpfr_init_set_str(r23751, "-1.6784217127088642e-162", 10, MPFR_RNDN); mpfr_init(r23752); mpfr_init_set_str(r23753, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23754, "2", 10, MPFR_RNDN); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init(r23768); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); } double f_fm(double g, double h, double a) { mpfr_set_d(r23750, g, MPFR_RNDN); ; mpfr_set_si(r23752, mpfr_cmp(r23750, r23751) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23755, a, MPFR_RNDN); mpfr_mul(r23756, r23754, r23755, MPFR_RNDN); mpfr_div(r23757, r23753, r23756, MPFR_RNDN); mpfr_neg(r23758, r23750, MPFR_RNDN); mpfr_mul(r23759, r23750, r23750, MPFR_RNDN); mpfr_set_d(r23760, h, MPFR_RNDN); mpfr_mul(r23761, r23760, r23760, MPFR_RNDN); mpfr_sub(r23762, r23759, r23761, MPFR_RNDN); mpfr_sqrt(r23763, r23762, MPFR_RNDN); mpfr_add(r23764, r23758, r23763, MPFR_RNDN); mpfr_mul(r23765, r23757, r23764, MPFR_RNDN); mpfr_cbrt(r23766, r23765, MPFR_RNDN); mpfr_div(r23767, r23761, r23764, MPFR_RNDN); mpfr_mul(r23768, r23757, r23767, MPFR_RNDN); mpfr_cbrt(r23769, r23768, MPFR_RNDN); mpfr_add(r23770, r23766, r23769, MPFR_RNDN); mpfr_sub(r23771, r23758, r23763, MPFR_RNDN); mpfr_div(r23772, r23761, r23771, MPFR_RNDN); mpfr_mul(r23773, r23757, r23772, MPFR_RNDN); mpfr_cbrt(r23774, r23773, MPFR_RNDN); mpfr_mul(r23775, r23757, r23771, MPFR_RNDN); mpfr_cbrt(r23776, r23775, MPFR_RNDN); mpfr_add(r23777, r23774, r23776, MPFR_RNDN); if (mpfr_get_si(r23752, MPFR_RNDN)) { mpfr_set(r23778, r23770, MPFR_RNDN); } else { mpfr_set(r23778, r23777, MPFR_RNDN); }; return mpfr_get_d(r23778, MPFR_RNDN); } static mpfr_t r23779, r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23779); mpfr_init_set_str(r23780, "-1.6784217127088642e-162", 10, MPFR_RNDN); mpfr_init(r23781); mpfr_init_set_str(r23782, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23783, "2", 10, MPFR_RNDN); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init(r23791); mpfr_init(r23792); mpfr_init(r23793); mpfr_init(r23794); mpfr_init(r23795); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init(r23801); mpfr_init(r23802); mpfr_init(r23803); mpfr_init(r23804); mpfr_init(r23805); mpfr_init(r23806); mpfr_init(r23807); } double f_dm(double g, double h, double a) { mpfr_set_d(r23779, g, MPFR_RNDN); ; mpfr_set_si(r23781, mpfr_cmp(r23779, r23780) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23784, a, MPFR_RNDN); mpfr_mul(r23785, r23783, r23784, MPFR_RNDN); mpfr_div(r23786, r23782, r23785, MPFR_RNDN); mpfr_neg(r23787, r23779, MPFR_RNDN); mpfr_mul(r23788, r23779, r23779, MPFR_RNDN); mpfr_set_d(r23789, h, MPFR_RNDN); mpfr_mul(r23790, r23789, r23789, MPFR_RNDN); mpfr_sub(r23791, r23788, r23790, MPFR_RNDN); mpfr_sqrt(r23792, r23791, MPFR_RNDN); mpfr_add(r23793, r23787, r23792, MPFR_RNDN); mpfr_mul(r23794, r23786, r23793, MPFR_RNDN); mpfr_cbrt(r23795, r23794, MPFR_RNDN); mpfr_div(r23796, r23790, r23793, MPFR_RNDN); mpfr_mul(r23797, r23786, r23796, MPFR_RNDN); mpfr_cbrt(r23798, r23797, MPFR_RNDN); mpfr_add(r23799, r23795, r23798, MPFR_RNDN); mpfr_sub(r23800, r23787, r23792, MPFR_RNDN); mpfr_div(r23801, r23790, r23800, MPFR_RNDN); mpfr_mul(r23802, r23786, r23801, MPFR_RNDN); mpfr_cbrt(r23803, r23802, MPFR_RNDN); mpfr_mul(r23804, r23786, r23800, MPFR_RNDN); mpfr_cbrt(r23805, r23804, MPFR_RNDN); mpfr_add(r23806, r23803, r23805, MPFR_RNDN); if (mpfr_get_si(r23781, MPFR_RNDN)) { mpfr_set(r23807, r23799, MPFR_RNDN); } else { mpfr_set(r23807, r23806, MPFR_RNDN); }; return mpfr_get_d(r23807, MPFR_RNDN); }