#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r23599 = 2; float r23600 = atan2(1.0, 0.0); float r23601 = r23599 * r23600; float r23602 = 3; float r23603 = r23601 / r23602; float r23604 = g; float r23605 = -r23604; float r23606 = h; float r23607 = r23605 / r23606; float r23608 = acos(r23607); float r23609 = r23608 / r23602; float r23610 = r23603 + r23609; float r23611 = cos(r23610); float r23612 = r23599 * r23611; return r23612; } double f_id(double g, double h) { double r23613 = 2; double r23614 = atan2(1.0, 0.0); double r23615 = r23613 * r23614; double r23616 = 3; double r23617 = r23615 / r23616; double r23618 = g; double r23619 = -r23618; double r23620 = h; double r23621 = r23619 / r23620; double r23622 = acos(r23621); double r23623 = r23622 / r23616; double r23624 = r23617 + r23623; double r23625 = cos(r23624); double r23626 = r23613 * r23625; return r23626; } double f_of(float g, float h) { float r23627 = 2; float r23628 = 3; float r23629 = r23627 / r23628; float r23630 = atan2(1.0, 0.0); float r23631 = g; float r23632 = h; float r23633 = r23631 / r23632; float r23634 = -r23633; float r23635 = acos(r23634); float r23636 = r23635 / r23628; float r23637 = fma(r23629, r23630, r23636); float r23638 = cos(r23637); float r23639 = r23638 * r23638; float r23640 = cbrt(r23639); float r23641 = cbrt(r23638); float r23642 = r23640 * r23641; float r23643 = r23627 * r23642; return r23643; } double f_od(double g, double h) { double r23644 = 2; double r23645 = 3; double r23646 = r23644 / r23645; double r23647 = atan2(1.0, 0.0); double r23648 = g; double r23649 = h; double r23650 = r23648 / r23649; double r23651 = -r23650; double r23652 = acos(r23651); double r23653 = r23652 / r23645; double r23654 = fma(r23646, r23647, r23653); double r23655 = cos(r23654); double r23656 = r23655 * r23655; double r23657 = cbrt(r23656); double r23658 = cbrt(r23655); double r23659 = r23657 * r23658; double r23660 = r23644 * r23659; return r23660; } 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 r23661, r23662, r23663, r23664, r23665, r23666, r23667, r23668, r23669, r23670, r23671, r23672, r23673, r23674; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23661, "2", 10, MPFR_RNDN); mpfr_init(r23662); mpfr_init(r23663); mpfr_init_set_str(r23664, "3", 10, MPFR_RNDN); mpfr_init(r23665); mpfr_init(r23666); mpfr_init(r23667); mpfr_init(r23668); mpfr_init(r23669); mpfr_init(r23670); mpfr_init(r23671); mpfr_init(r23672); mpfr_init(r23673); mpfr_init(r23674); } double f_im(double g, double h) { ; mpfr_const_pi(r23662, MPFR_RNDN); mpfr_mul(r23663, r23661, r23662, MPFR_RNDN); ; mpfr_div(r23665, r23663, r23664, MPFR_RNDN); mpfr_set_d(r23666, g, MPFR_RNDN); mpfr_neg(r23667, r23666, MPFR_RNDN); mpfr_set_d(r23668, h, MPFR_RNDN); mpfr_div(r23669, r23667, r23668, MPFR_RNDN); mpfr_acos(r23670, r23669, MPFR_RNDN); mpfr_div(r23671, r23670, r23664, MPFR_RNDN); mpfr_add(r23672, r23665, r23671, MPFR_RNDN); mpfr_cos(r23673, r23672, MPFR_RNDN); mpfr_mul(r23674, r23661, r23673, MPFR_RNDN); return mpfr_get_d(r23674, MPFR_RNDN); } static mpfr_t r23675, r23676, r23677, r23678, r23679, r23680, r23681, r23682, r23683, r23684, r23685, r23686, r23687, r23688, r23689, r23690, r23691; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23675, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23676, "3", 10, MPFR_RNDN); mpfr_init(r23677); mpfr_init(r23678); mpfr_init(r23679); mpfr_init(r23680); mpfr_init(r23681); mpfr_init(r23682); mpfr_init(r23683); mpfr_init(r23684); mpfr_init(r23685); mpfr_init(r23686); mpfr_init(r23687); mpfr_init(r23688); mpfr_init(r23689); mpfr_init(r23690); mpfr_init(r23691); } double f_fm(double g, double h) { ; ; mpfr_div(r23677, r23675, r23676, MPFR_RNDN); mpfr_const_pi(r23678, MPFR_RNDN); mpfr_set_d(r23679, g, MPFR_RNDN); mpfr_set_d(r23680, h, MPFR_RNDN); mpfr_div(r23681, r23679, r23680, MPFR_RNDN); mpfr_neg(r23682, r23681, MPFR_RNDN); mpfr_acos(r23683, r23682, MPFR_RNDN); mpfr_div(r23684, r23683, r23676, MPFR_RNDN); mpfr_fma(r23685, r23677, r23678, r23684, MPFR_RNDN); mpfr_cos(r23686, r23685, MPFR_RNDN); mpfr_mul(r23687, r23686, r23686, MPFR_RNDN); mpfr_cbrt(r23688, r23687, MPFR_RNDN); mpfr_cbrt(r23689, r23686, MPFR_RNDN); mpfr_mul(r23690, r23688, r23689, MPFR_RNDN); mpfr_mul(r23691, r23675, r23690, MPFR_RNDN); return mpfr_get_d(r23691, MPFR_RNDN); } static mpfr_t r23692, r23693, r23694, r23695, r23696, r23697, r23698, r23699, r23700, r23701, r23702, r23703, r23704, r23705, r23706, r23707, r23708; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23692, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23693, "3", 10, MPFR_RNDN); mpfr_init(r23694); mpfr_init(r23695); mpfr_init(r23696); mpfr_init(r23697); mpfr_init(r23698); mpfr_init(r23699); mpfr_init(r23700); mpfr_init(r23701); mpfr_init(r23702); mpfr_init(r23703); mpfr_init(r23704); mpfr_init(r23705); mpfr_init(r23706); mpfr_init(r23707); mpfr_init(r23708); } double f_dm(double g, double h) { ; ; mpfr_div(r23694, r23692, r23693, MPFR_RNDN); mpfr_const_pi(r23695, MPFR_RNDN); mpfr_set_d(r23696, g, MPFR_RNDN); mpfr_set_d(r23697, h, MPFR_RNDN); mpfr_div(r23698, r23696, r23697, MPFR_RNDN); mpfr_neg(r23699, r23698, MPFR_RNDN); mpfr_acos(r23700, r23699, MPFR_RNDN); mpfr_div(r23701, r23700, r23693, MPFR_RNDN); mpfr_fma(r23702, r23694, r23695, r23701, MPFR_RNDN); mpfr_cos(r23703, r23702, MPFR_RNDN); mpfr_mul(r23704, r23703, r23703, MPFR_RNDN); mpfr_cbrt(r23705, r23704, MPFR_RNDN); mpfr_cbrt(r23706, r23703, MPFR_RNDN); mpfr_mul(r23707, r23705, r23706, MPFR_RNDN); mpfr_mul(r23708, r23692, r23707, MPFR_RNDN); return mpfr_get_d(r23708, MPFR_RNDN); }