#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r23512 = 1; float r23513 = 2; float r23514 = a; float r23515 = r23513 * r23514; float r23516 = r23512 / r23515; float r23517 = g; float r23518 = -r23517; float r23519 = r23517 * r23517; float r23520 = h; float r23521 = r23520 * r23520; float r23522 = r23519 - r23521; float r23523 = sqrt(r23522); float r23524 = r23518 + r23523; float r23525 = r23516 * r23524; float r23526 = cbrt(r23525); float r23527 = r23518 - r23523; float r23528 = r23516 * r23527; float r23529 = cbrt(r23528); float r23530 = r23526 + r23529; return r23530; } double f_id(double g, double h, double a) { double r23531 = 1; double r23532 = 2; double r23533 = a; double r23534 = r23532 * r23533; double r23535 = r23531 / r23534; double r23536 = g; double r23537 = -r23536; double r23538 = r23536 * r23536; double r23539 = h; double r23540 = r23539 * r23539; double r23541 = r23538 - r23540; double r23542 = sqrt(r23541); double r23543 = r23537 + r23542; double r23544 = r23535 * r23543; double r23545 = cbrt(r23544); double r23546 = r23537 - r23542; double r23547 = r23535 * r23546; double r23548 = cbrt(r23547); double r23549 = r23545 + r23548; return r23549; } double f_of(float g, float h, float a) { float r23550 = g; float r23551 = -2.0103451662450893e-162; bool r23552 = r23550 <= r23551; float r23553 = 1; float r23554 = 2; float r23555 = a; float r23556 = r23554 * r23555; float r23557 = r23553 / r23556; float r23558 = -r23550; float r23559 = r23550 * r23550; float r23560 = h; float r23561 = r23560 * r23560; float r23562 = r23559 - r23561; float r23563 = sqrt(r23562); float r23564 = r23558 + r23563; float r23565 = r23557 * r23564; float r23566 = cbrt(r23565); float r23567 = r23561 / r23564; float r23568 = r23557 * r23567; float r23569 = cbrt(r23568); float r23570 = r23566 + r23569; float r23571 = r23558 - r23563; float r23572 = r23561 / r23571; float r23573 = r23557 * r23572; float r23574 = cbrt(r23573); float r23575 = r23557 * r23571; float r23576 = cbrt(r23575); float r23577 = r23574 + r23576; float r23578 = r23552 ? r23570 : r23577; return r23578; } double f_od(double g, double h, double a) { double r23579 = g; double r23580 = -2.0103451662450893e-162; bool r23581 = r23579 <= r23580; double r23582 = 1; double r23583 = 2; double r23584 = a; double r23585 = r23583 * r23584; double r23586 = r23582 / r23585; double r23587 = -r23579; double r23588 = r23579 * r23579; double r23589 = h; double r23590 = r23589 * r23589; double r23591 = r23588 - r23590; double r23592 = sqrt(r23591); double r23593 = r23587 + r23592; double r23594 = r23586 * r23593; double r23595 = cbrt(r23594); double r23596 = r23590 / r23593; double r23597 = r23586 * r23596; double r23598 = cbrt(r23597); double r23599 = r23595 + r23598; double r23600 = r23587 - r23592; double r23601 = r23590 / r23600; double r23602 = r23586 * r23601; double r23603 = cbrt(r23602); double r23604 = r23586 * r23600; double r23605 = cbrt(r23604); double r23606 = r23603 + r23605; double r23607 = r23581 ? r23599 : r23606; return r23607; } 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 r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init_set_str(r23608, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23609, "2", 10, MPFR_RNDN); 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); mpfr_init(r23626); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r23610, a, MPFR_RNDN); mpfr_mul(r23611, r23609, r23610, MPFR_RNDN); mpfr_div(r23612, r23608, r23611, MPFR_RNDN); mpfr_set_d(r23613, g, MPFR_RNDN); mpfr_neg(r23614, r23613, MPFR_RNDN); mpfr_mul(r23615, r23613, r23613, MPFR_RNDN); mpfr_set_d(r23616, h, MPFR_RNDN); mpfr_mul(r23617, r23616, r23616, MPFR_RNDN); mpfr_sub(r23618, r23615, r23617, MPFR_RNDN); mpfr_sqrt(r23619, r23618, MPFR_RNDN); mpfr_add(r23620, r23614, r23619, MPFR_RNDN); mpfr_mul(r23621, r23612, r23620, MPFR_RNDN); mpfr_cbrt(r23622, r23621, MPFR_RNDN); mpfr_sub(r23623, r23614, r23619, MPFR_RNDN); mpfr_mul(r23624, r23612, r23623, MPFR_RNDN); mpfr_cbrt(r23625, r23624, MPFR_RNDN); mpfr_add(r23626, r23622, r23625, MPFR_RNDN); return mpfr_get_d(r23626, MPFR_RNDN); } static mpfr_t 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, r23655; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23627); mpfr_init_set_str(r23628, "-2.0103451662450893e-162", 10, MPFR_RNDN); mpfr_init(r23629); mpfr_init_set_str(r23630, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23631, "2", 10, MPFR_RNDN); 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); mpfr_init(r23655); } double f_fm(double g, double h, double a) { mpfr_set_d(r23627, g, MPFR_RNDN); ; mpfr_set_si(r23629, mpfr_cmp(r23627, r23628) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23632, a, MPFR_RNDN); mpfr_mul(r23633, r23631, r23632, MPFR_RNDN); mpfr_div(r23634, r23630, r23633, MPFR_RNDN); mpfr_neg(r23635, r23627, MPFR_RNDN); mpfr_mul(r23636, r23627, r23627, MPFR_RNDN); mpfr_set_d(r23637, h, MPFR_RNDN); mpfr_mul(r23638, r23637, r23637, MPFR_RNDN); mpfr_sub(r23639, r23636, r23638, MPFR_RNDN); mpfr_sqrt(r23640, r23639, MPFR_RNDN); mpfr_add(r23641, r23635, r23640, MPFR_RNDN); mpfr_mul(r23642, r23634, r23641, MPFR_RNDN); mpfr_cbrt(r23643, r23642, MPFR_RNDN); mpfr_div(r23644, r23638, r23641, MPFR_RNDN); mpfr_mul(r23645, r23634, r23644, MPFR_RNDN); mpfr_cbrt(r23646, r23645, MPFR_RNDN); mpfr_add(r23647, r23643, r23646, MPFR_RNDN); mpfr_sub(r23648, r23635, r23640, MPFR_RNDN); mpfr_div(r23649, r23638, r23648, MPFR_RNDN); mpfr_mul(r23650, r23634, r23649, MPFR_RNDN); mpfr_cbrt(r23651, r23650, MPFR_RNDN); mpfr_mul(r23652, r23634, r23648, MPFR_RNDN); mpfr_cbrt(r23653, r23652, MPFR_RNDN); mpfr_add(r23654, r23651, r23653, MPFR_RNDN); if (mpfr_get_si(r23629, MPFR_RNDN)) { mpfr_set(r23655, r23647, MPFR_RNDN); } else { mpfr_set(r23655, r23654, MPFR_RNDN); }; return mpfr_get_d(r23655, MPFR_RNDN); } static mpfr_t r23656, r23657, r23658, r23659, r23660, r23661, r23662, r23663, r23664, r23665, r23666, r23667, r23668, r23669, r23670, r23671, r23672, r23673, r23674, r23675, r23676, r23677, r23678, r23679, r23680, r23681, r23682, r23683, r23684; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23656); mpfr_init_set_str(r23657, "-2.0103451662450893e-162", 10, MPFR_RNDN); mpfr_init(r23658); mpfr_init_set_str(r23659, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23660, "2", 10, MPFR_RNDN); mpfr_init(r23661); mpfr_init(r23662); mpfr_init(r23663); mpfr_init(r23664); 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); mpfr_init(r23675); mpfr_init(r23676); mpfr_init(r23677); mpfr_init(r23678); mpfr_init(r23679); mpfr_init(r23680); mpfr_init(r23681); mpfr_init(r23682); mpfr_init(r23683); mpfr_init(r23684); } double f_dm(double g, double h, double a) { mpfr_set_d(r23656, g, MPFR_RNDN); ; mpfr_set_si(r23658, mpfr_cmp(r23656, r23657) <= 0, MPFR_RNDN); ; ; mpfr_set_d(r23661, a, MPFR_RNDN); mpfr_mul(r23662, r23660, r23661, MPFR_RNDN); mpfr_div(r23663, r23659, r23662, MPFR_RNDN); mpfr_neg(r23664, r23656, MPFR_RNDN); mpfr_mul(r23665, r23656, r23656, MPFR_RNDN); mpfr_set_d(r23666, h, MPFR_RNDN); mpfr_mul(r23667, r23666, r23666, MPFR_RNDN); mpfr_sub(r23668, r23665, r23667, MPFR_RNDN); mpfr_sqrt(r23669, r23668, MPFR_RNDN); mpfr_add(r23670, r23664, r23669, MPFR_RNDN); mpfr_mul(r23671, r23663, r23670, MPFR_RNDN); mpfr_cbrt(r23672, r23671, MPFR_RNDN); mpfr_div(r23673, r23667, r23670, MPFR_RNDN); mpfr_mul(r23674, r23663, r23673, MPFR_RNDN); mpfr_cbrt(r23675, r23674, MPFR_RNDN); mpfr_add(r23676, r23672, r23675, MPFR_RNDN); mpfr_sub(r23677, r23664, r23669, MPFR_RNDN); mpfr_div(r23678, r23667, r23677, MPFR_RNDN); mpfr_mul(r23679, r23663, r23678, MPFR_RNDN); mpfr_cbrt(r23680, r23679, MPFR_RNDN); mpfr_mul(r23681, r23663, r23677, MPFR_RNDN); mpfr_cbrt(r23682, r23681, MPFR_RNDN); mpfr_add(r23683, r23680, r23682, MPFR_RNDN); if (mpfr_get_si(r23658, MPFR_RNDN)) { mpfr_set(r23684, r23676, MPFR_RNDN); } else { mpfr_set(r23684, r23683, MPFR_RNDN); }; return mpfr_get_d(r23684, MPFR_RNDN); }