#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r23482 = 1; float r23483 = 2; float r23484 = a; float r23485 = r23483 * r23484; float r23486 = r23482 / r23485; float r23487 = g; float r23488 = -r23487; float r23489 = r23487 * r23487; float r23490 = h; float r23491 = r23490 * r23490; float r23492 = r23489 - r23491; float r23493 = sqrt(r23492); float r23494 = r23488 + r23493; float r23495 = r23486 * r23494; float r23496 = cbrt(r23495); float r23497 = r23488 - r23493; float r23498 = r23486 * r23497; float r23499 = cbrt(r23498); float r23500 = r23496 + r23499; return r23500; } double f_id(double g, double h, double a) { double r23501 = 1; double r23502 = 2; double r23503 = a; double r23504 = r23502 * r23503; double r23505 = r23501 / r23504; double r23506 = g; double r23507 = -r23506; double r23508 = r23506 * r23506; double r23509 = h; double r23510 = r23509 * r23509; double r23511 = r23508 - r23510; double r23512 = sqrt(r23511); double r23513 = r23507 + r23512; double r23514 = r23505 * r23513; double r23515 = cbrt(r23514); double r23516 = r23507 - r23512; double r23517 = r23505 * r23516; double r23518 = cbrt(r23517); double r23519 = r23515 + r23518; return r23519; } double f_of(float g, float h, float a) { float r23520 = g; float r23521 = -1.683025889411105e-162; bool r23522 = r23520 <= r23521; float r23523 = 1; float r23524 = 2; float r23525 = a; float r23526 = r23524 * r23525; float r23527 = r23523 / r23526; float r23528 = -r23520; float r23529 = r23520 * r23520; float r23530 = h; float r23531 = r23530 * r23530; float r23532 = r23529 - r23531; float r23533 = sqrt(r23532); float r23534 = r23528 + r23533; float r23535 = r23527 * r23534; float r23536 = cbrt(r23535); float r23537 = r23531 / r23534; float r23538 = r23527 * r23537; float r23539 = cbrt(r23538); float r23540 = r23536 + r23539; float r23541 = r23528 - r23533; float r23542 = r23531 / r23541; float r23543 = r23527 * r23542; float r23544 = cbrt(r23543); float r23545 = r23527 * r23541; float r23546 = cbrt(r23545); float r23547 = r23544 + r23546; float r23548 = r23522 ? r23540 : r23547; return r23548; } double f_od(double g, double h, double a) { double r23549 = g; double r23550 = -1.683025889411105e-162; bool r23551 = r23549 <= r23550; double r23552 = 1; double r23553 = 2; double r23554 = a; double r23555 = r23553 * r23554; double r23556 = r23552 / r23555; double r23557 = -r23549; double r23558 = r23549 * r23549; double r23559 = h; double r23560 = r23559 * r23559; double r23561 = r23558 - r23560; double r23562 = sqrt(r23561); double r23563 = r23557 + r23562; double r23564 = r23556 * r23563; double r23565 = cbrt(r23564); double r23566 = r23560 / r23563; double r23567 = r23556 * r23566; double r23568 = cbrt(r23567); double r23569 = r23565 + r23568; double r23570 = r23557 - r23562; double r23571 = r23560 / r23570; double r23572 = r23556 * r23571; double r23573 = cbrt(r23572); double r23574 = r23556 * r23570; double r23575 = cbrt(r23574); double r23576 = r23573 + r23575; double r23577 = r23551 ? r23569 : r23576; return r23577; } 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 r23578, r23579, r23580, r23581, r23582, r23583, r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593, r23594, r23595, r23596; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init_set_str(r23578, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23579, "2", 10, MPFR_RNDN); mpfr_init(r23580); mpfr_init(r23581); mpfr_init(r23582); mpfr_init(r23583); mpfr_init(r23584); mpfr_init(r23585); mpfr_init(r23586); mpfr_init(r23587); mpfr_init(r23588); mpfr_init(r23589); mpfr_init(r23590); mpfr_init(r23591); mpfr_init(r23592); mpfr_init(r23593); mpfr_init(r23594); mpfr_init(r23595); mpfr_init(r23596); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r23580, a, MPFR_RNDN); mpfr_mul(r23581, r23579, r23580, MPFR_RNDN); mpfr_div(r23582, r23578, r23581, MPFR_RNDN); mpfr_set_d(r23583, g, MPFR_RNDN); mpfr_neg(r23584, r23583, MPFR_RNDN); mpfr_mul(r23585, r23583, r23583, MPFR_RNDN); mpfr_set_d(r23586, h, MPFR_RNDN); mpfr_mul(r23587, r23586, r23586, MPFR_RNDN); mpfr_sub(r23588, r23585, r23587, MPFR_RNDN); mpfr_sqrt(r23589, r23588, MPFR_RNDN); mpfr_add(r23590, r23584, r23589, MPFR_RNDN); mpfr_mul(r23591, r23582, r23590, MPFR_RNDN); mpfr_cbrt(r23592, r23591, MPFR_RNDN); mpfr_sub(r23593, r23584, r23589, MPFR_RNDN); mpfr_mul(r23594, r23582, r23593, MPFR_RNDN); mpfr_cbrt(r23595, r23594, MPFR_RNDN); mpfr_add(r23596, r23592, r23595, MPFR_RNDN); return mpfr_get_d(r23596, MPFR_RNDN); } static mpfr_t r23597, r23598, r23599, r23600, r23601, r23602, r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23597); mpfr_init_set_str(r23598, "-1.683025889411105e-162", 10, MPFR_RNDN); mpfr_init(r23599); mpfr_init_set_str(r23600, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23601, "2", 10, MPFR_RNDN); mpfr_init(r23602); mpfr_init(r23603); mpfr_init(r23604); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init(r23611); mpfr_init(r23612); mpfr_init(r23613); mpfr_init(r23614); mpfr_init(r23615); mpfr_init(r23616); mpfr_init(r23617); mpfr_init(r23618); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); } double f_fm(double g, double h, double a) { mpfr_set_d(r23597, g, MPFR_RNDN); ; mpfr_set_si(r23599, mpfr_cmp(r23597, r23598) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23602, a, MPFR_RNDN); mpfr_mul(r23603, r23601, r23602, MPFR_RNDN); mpfr_div(r23604, r23600, r23603, MPFR_RNDN); mpfr_neg(r23605, r23597, MPFR_RNDN); mpfr_mul(r23606, r23597, r23597, MPFR_RNDN); mpfr_set_d(r23607, h, MPFR_RNDN); mpfr_mul(r23608, r23607, r23607, MPFR_RNDN); mpfr_sub(r23609, r23606, r23608, MPFR_RNDN); mpfr_sqrt(r23610, r23609, MPFR_RNDN); mpfr_add(r23611, r23605, r23610, MPFR_RNDN); mpfr_mul(r23612, r23604, r23611, MPFR_RNDN); mpfr_cbrt(r23613, r23612, MPFR_RNDN); mpfr_div(r23614, r23608, r23611, MPFR_RNDN); mpfr_mul(r23615, r23604, r23614, MPFR_RNDN); mpfr_cbrt(r23616, r23615, MPFR_RNDN); mpfr_add(r23617, r23613, r23616, MPFR_RNDN); mpfr_sub(r23618, r23605, r23610, MPFR_RNDN); mpfr_div(r23619, r23608, r23618, MPFR_RNDN); mpfr_mul(r23620, r23604, r23619, MPFR_RNDN); mpfr_cbrt(r23621, r23620, MPFR_RNDN); mpfr_mul(r23622, r23604, r23618, MPFR_RNDN); mpfr_cbrt(r23623, r23622, MPFR_RNDN); mpfr_add(r23624, r23621, r23623, MPFR_RNDN); if (mpfr_get_si(r23599, MPFR_RNDN)) { mpfr_set(r23625, r23617, MPFR_RNDN); } else { mpfr_set(r23625, r23624, MPFR_RNDN); }; return mpfr_get_d(r23625, MPFR_RNDN); } static mpfr_t r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648, r23649, r23650, r23651, r23652, r23653, r23654; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23626); mpfr_init_set_str(r23627, "-1.683025889411105e-162", 10, MPFR_RNDN); mpfr_init(r23628); mpfr_init_set_str(r23629, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23630, "2", 10, MPFR_RNDN); mpfr_init(r23631); mpfr_init(r23632); mpfr_init(r23633); mpfr_init(r23634); mpfr_init(r23635); mpfr_init(r23636); mpfr_init(r23637); mpfr_init(r23638); mpfr_init(r23639); mpfr_init(r23640); mpfr_init(r23641); mpfr_init(r23642); mpfr_init(r23643); mpfr_init(r23644); mpfr_init(r23645); mpfr_init(r23646); mpfr_init(r23647); mpfr_init(r23648); mpfr_init(r23649); mpfr_init(r23650); mpfr_init(r23651); mpfr_init(r23652); mpfr_init(r23653); mpfr_init(r23654); } double f_dm(double g, double h, double a) { mpfr_set_d(r23626, g, MPFR_RNDN); ; mpfr_set_si(r23628, mpfr_cmp(r23626, r23627) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23631, a, MPFR_RNDN); mpfr_mul(r23632, r23630, r23631, MPFR_RNDN); mpfr_div(r23633, r23629, r23632, MPFR_RNDN); mpfr_neg(r23634, r23626, MPFR_RNDN); mpfr_mul(r23635, r23626, r23626, MPFR_RNDN); mpfr_set_d(r23636, h, MPFR_RNDN); mpfr_mul(r23637, r23636, r23636, MPFR_RNDN); mpfr_sub(r23638, r23635, r23637, MPFR_RNDN); mpfr_sqrt(r23639, r23638, MPFR_RNDN); mpfr_add(r23640, r23634, r23639, MPFR_RNDN); mpfr_mul(r23641, r23633, r23640, MPFR_RNDN); mpfr_cbrt(r23642, r23641, MPFR_RNDN); mpfr_div(r23643, r23637, r23640, MPFR_RNDN); mpfr_mul(r23644, r23633, r23643, MPFR_RNDN); mpfr_cbrt(r23645, r23644, MPFR_RNDN); mpfr_add(r23646, r23642, r23645, MPFR_RNDN); mpfr_sub(r23647, r23634, r23639, MPFR_RNDN); mpfr_div(r23648, r23637, r23647, MPFR_RNDN); mpfr_mul(r23649, r23633, r23648, MPFR_RNDN); mpfr_cbrt(r23650, r23649, MPFR_RNDN); mpfr_mul(r23651, r23633, r23647, MPFR_RNDN); mpfr_cbrt(r23652, r23651, MPFR_RNDN); mpfr_add(r23653, r23650, r23652, MPFR_RNDN); if (mpfr_get_si(r23628, MPFR_RNDN)) { mpfr_set(r23654, r23646, MPFR_RNDN); } else { mpfr_set(r23654, r23653, MPFR_RNDN); }; return mpfr_get_d(r23654, MPFR_RNDN); }