#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r24433 = 2; float r24434 = atan2(1.0, 0.0); float r24435 = r24433 * r24434; float r24436 = 3; float r24437 = r24435 / r24436; float r24438 = g; float r24439 = -r24438; float r24440 = h; float r24441 = r24439 / r24440; float r24442 = acos(r24441); float r24443 = r24442 / r24436; float r24444 = r24437 + r24443; float r24445 = cos(r24444); float r24446 = r24433 * r24445; return r24446; } double f_id(double g, double h) { double r24447 = 2; double r24448 = atan2(1.0, 0.0); double r24449 = r24447 * r24448; double r24450 = 3; double r24451 = r24449 / r24450; double r24452 = g; double r24453 = -r24452; double r24454 = h; double r24455 = r24453 / r24454; double r24456 = acos(r24455); double r24457 = r24456 / r24450; double r24458 = r24451 + r24457; double r24459 = cos(r24458); double r24460 = r24447 * r24459; return r24460; } double f_of(float g, float h) { float r24461 = 2; float r24462 = atan2(1.0, 0.0); float r24463 = 3; float r24464 = r24463 / r24461; float r24465 = r24462 / r24464; float r24466 = g; float r24467 = -r24466; float r24468 = h; float r24469 = r24467 / r24468; float r24470 = acos(r24469); float r24471 = r24470 / r24463; float r24472 = r24465 + r24471; float r24473 = cos(r24472); float r24474 = r24473 * r24473; float r24475 = cbrt(r24474); float r24476 = r24462 / r24463; float r24477 = r24461 * r24476; float r24478 = r24477 + r24471; float r24479 = cos(r24478); float r24480 = cbrt(r24479); float r24481 = r24475 * r24480; float r24482 = r24461 * r24481; return r24482; } double f_od(double g, double h) { double r24483 = 2; double r24484 = atan2(1.0, 0.0); double r24485 = 3; double r24486 = r24485 / r24483; double r24487 = r24484 / r24486; double r24488 = g; double r24489 = -r24488; double r24490 = h; double r24491 = r24489 / r24490; double r24492 = acos(r24491); double r24493 = r24492 / r24485; double r24494 = r24487 + r24493; double r24495 = cos(r24494); double r24496 = r24495 * r24495; double r24497 = cbrt(r24496); double r24498 = r24484 / r24485; double r24499 = r24483 * r24498; double r24500 = r24499 + r24493; double r24501 = cos(r24500); double r24502 = cbrt(r24501); double r24503 = r24497 * r24502; double r24504 = r24483 * r24503; return r24504; } 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 r24505, r24506, r24507, r24508, r24509, r24510, r24511, r24512, r24513, r24514, r24515, r24516, r24517, r24518; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24505, "2", 10, MPFR_RNDN); mpfr_init(r24506); mpfr_init(r24507); mpfr_init_set_str(r24508, "3", 10, MPFR_RNDN); mpfr_init(r24509); mpfr_init(r24510); mpfr_init(r24511); mpfr_init(r24512); mpfr_init(r24513); mpfr_init(r24514); mpfr_init(r24515); mpfr_init(r24516); mpfr_init(r24517); mpfr_init(r24518); } double f_im(double g, double h) { ; mpfr_const_pi(r24506, MPFR_RNDN); mpfr_mul(r24507, r24505, r24506, MPFR_RNDN); ; mpfr_div(r24509, r24507, r24508, MPFR_RNDN); mpfr_set_d(r24510, g, MPFR_RNDN); mpfr_neg(r24511, r24510, MPFR_RNDN); mpfr_set_d(r24512, h, MPFR_RNDN); mpfr_div(r24513, r24511, r24512, MPFR_RNDN); mpfr_acos(r24514, r24513, MPFR_RNDN); mpfr_div(r24515, r24514, r24508, MPFR_RNDN); mpfr_add(r24516, r24509, r24515, MPFR_RNDN); mpfr_cos(r24517, r24516, MPFR_RNDN); mpfr_mul(r24518, r24505, r24517, MPFR_RNDN); return mpfr_get_d(r24518, MPFR_RNDN); } static mpfr_t r24519, r24520, r24521, r24522, r24523, r24524, r24525, r24526, r24527, r24528, r24529, r24530, r24531, r24532, r24533, r24534, r24535, r24536, r24537, r24538, r24539, r24540; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24519, "2", 10, MPFR_RNDN); mpfr_init(r24520); mpfr_init_set_str(r24521, "3", 10, MPFR_RNDN); mpfr_init(r24522); mpfr_init(r24523); mpfr_init(r24524); mpfr_init(r24525); mpfr_init(r24526); mpfr_init(r24527); mpfr_init(r24528); mpfr_init(r24529); mpfr_init(r24530); mpfr_init(r24531); mpfr_init(r24532); mpfr_init(r24533); mpfr_init(r24534); mpfr_init(r24535); mpfr_init(r24536); mpfr_init(r24537); mpfr_init(r24538); mpfr_init(r24539); mpfr_init(r24540); } double f_fm(double g, double h) { ; mpfr_const_pi(r24520, MPFR_RNDN); ; mpfr_div(r24522, r24521, r24519, MPFR_RNDN); mpfr_div(r24523, r24520, r24522, MPFR_RNDN); mpfr_set_d(r24524, g, MPFR_RNDN); mpfr_neg(r24525, r24524, MPFR_RNDN); mpfr_set_d(r24526, h, MPFR_RNDN); mpfr_div(r24527, r24525, r24526, MPFR_RNDN); mpfr_acos(r24528, r24527, MPFR_RNDN); mpfr_div(r24529, r24528, r24521, MPFR_RNDN); mpfr_add(r24530, r24523, r24529, MPFR_RNDN); mpfr_cos(r24531, r24530, MPFR_RNDN); mpfr_mul(r24532, r24531, r24531, MPFR_RNDN); mpfr_cbrt(r24533, r24532, MPFR_RNDN); mpfr_div(r24534, r24520, r24521, MPFR_RNDN); mpfr_mul(r24535, r24519, r24534, MPFR_RNDN); mpfr_add(r24536, r24535, r24529, MPFR_RNDN); mpfr_cos(r24537, r24536, MPFR_RNDN); mpfr_cbrt(r24538, r24537, MPFR_RNDN); mpfr_mul(r24539, r24533, r24538, MPFR_RNDN); mpfr_mul(r24540, r24519, r24539, MPFR_RNDN); return mpfr_get_d(r24540, MPFR_RNDN); } static mpfr_t r24541, r24542, r24543, r24544, r24545, r24546, r24547, r24548, r24549, r24550, r24551, r24552, r24553, r24554, r24555, r24556, r24557, r24558, r24559, r24560, r24561, r24562; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24541, "2", 10, MPFR_RNDN); mpfr_init(r24542); mpfr_init_set_str(r24543, "3", 10, MPFR_RNDN); mpfr_init(r24544); mpfr_init(r24545); mpfr_init(r24546); mpfr_init(r24547); mpfr_init(r24548); mpfr_init(r24549); mpfr_init(r24550); mpfr_init(r24551); mpfr_init(r24552); mpfr_init(r24553); mpfr_init(r24554); mpfr_init(r24555); mpfr_init(r24556); mpfr_init(r24557); mpfr_init(r24558); mpfr_init(r24559); mpfr_init(r24560); mpfr_init(r24561); mpfr_init(r24562); } double f_dm(double g, double h) { ; mpfr_const_pi(r24542, MPFR_RNDN); ; mpfr_div(r24544, r24543, r24541, MPFR_RNDN); mpfr_div(r24545, r24542, r24544, MPFR_RNDN); mpfr_set_d(r24546, g, MPFR_RNDN); mpfr_neg(r24547, r24546, MPFR_RNDN); mpfr_set_d(r24548, h, MPFR_RNDN); mpfr_div(r24549, r24547, r24548, MPFR_RNDN); mpfr_acos(r24550, r24549, MPFR_RNDN); mpfr_div(r24551, r24550, r24543, MPFR_RNDN); mpfr_add(r24552, r24545, r24551, MPFR_RNDN); mpfr_cos(r24553, r24552, MPFR_RNDN); mpfr_mul(r24554, r24553, r24553, MPFR_RNDN); mpfr_cbrt(r24555, r24554, MPFR_RNDN); mpfr_div(r24556, r24542, r24543, MPFR_RNDN); mpfr_mul(r24557, r24541, r24556, MPFR_RNDN); mpfr_add(r24558, r24557, r24551, MPFR_RNDN); mpfr_cos(r24559, r24558, MPFR_RNDN); mpfr_cbrt(r24560, r24559, MPFR_RNDN); mpfr_mul(r24561, r24555, r24560, MPFR_RNDN); mpfr_mul(r24562, r24541, r24561, MPFR_RNDN); return mpfr_get_d(r24562, MPFR_RNDN); }