#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r24376 = 2; float r24377 = atan2(1.0, 0.0); float r24378 = r24376 * r24377; float r24379 = 3; float r24380 = r24378 / r24379; float r24381 = g; float r24382 = -r24381; float r24383 = h; float r24384 = r24382 / r24383; float r24385 = acos(r24384); float r24386 = r24385 / r24379; float r24387 = r24380 + r24386; float r24388 = cos(r24387); float r24389 = r24376 * r24388; return r24389; } double f_id(double g, double h) { double r24390 = 2; double r24391 = atan2(1.0, 0.0); double r24392 = r24390 * r24391; double r24393 = 3; double r24394 = r24392 / r24393; double r24395 = g; double r24396 = -r24395; double r24397 = h; double r24398 = r24396 / r24397; double r24399 = acos(r24398); double r24400 = r24399 / r24393; double r24401 = r24394 + r24400; double r24402 = cos(r24401); double r24403 = r24390 * r24402; return r24403; } double f_of(float g, float h) { float r24404 = 2; float r24405 = 3; float r24406 = r24404 / r24405; float r24407 = atan2(1.0, 0.0); float r24408 = g; float r24409 = h; float r24410 = r24408 / r24409; float r24411 = -r24410; float r24412 = acos(r24411); float r24413 = r24412 / r24405; float r24414 = fma(r24406, r24407, r24413); float r24415 = cos(r24414); float r24416 = expm1(r24415); float r24417 = log1p(r24416); float r24418 = r24404 * r24417; return r24418; } double f_od(double g, double h) { double r24419 = 2; double r24420 = 3; double r24421 = r24419 / r24420; double r24422 = atan2(1.0, 0.0); double r24423 = g; double r24424 = h; double r24425 = r24423 / r24424; double r24426 = -r24425; double r24427 = acos(r24426); double r24428 = r24427 / r24420; double r24429 = fma(r24421, r24422, r24428); double r24430 = cos(r24429); double r24431 = expm1(r24430); double r24432 = log1p(r24431); double r24433 = r24419 * r24432; return r24433; } 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 r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24434, "2", 10, MPFR_RNDN); mpfr_init(r24435); mpfr_init(r24436); mpfr_init_set_str(r24437, "3", 10, MPFR_RNDN); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); } double f_im(double g, double h) { ; mpfr_const_pi(r24435, MPFR_RNDN); mpfr_mul(r24436, r24434, r24435, MPFR_RNDN); ; mpfr_div(r24438, r24436, r24437, MPFR_RNDN); mpfr_set_d(r24439, g, MPFR_RNDN); mpfr_neg(r24440, r24439, MPFR_RNDN); mpfr_set_d(r24441, h, MPFR_RNDN); mpfr_div(r24442, r24440, r24441, MPFR_RNDN); mpfr_acos(r24443, r24442, MPFR_RNDN); mpfr_div(r24444, r24443, r24437, MPFR_RNDN); mpfr_add(r24445, r24438, r24444, MPFR_RNDN); mpfr_cos(r24446, r24445, MPFR_RNDN); mpfr_mul(r24447, r24434, r24446, MPFR_RNDN); return mpfr_get_d(r24447, MPFR_RNDN); } static mpfr_t r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24448, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24449, "3", 10, MPFR_RNDN); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init(r24458); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); } double f_fm(double g, double h) { ; ; mpfr_div(r24450, r24448, r24449, MPFR_RNDN); mpfr_const_pi(r24451, MPFR_RNDN); mpfr_set_d(r24452, g, MPFR_RNDN); mpfr_set_d(r24453, h, MPFR_RNDN); mpfr_div(r24454, r24452, r24453, MPFR_RNDN); mpfr_neg(r24455, r24454, MPFR_RNDN); mpfr_acos(r24456, r24455, MPFR_RNDN); mpfr_div(r24457, r24456, r24449, MPFR_RNDN); mpfr_fma(r24458, r24450, r24451, r24457, MPFR_RNDN); mpfr_cos(r24459, r24458, MPFR_RNDN); mpfr_expm1(r24460, r24459, MPFR_RNDN); mpfr_log1p(r24461, r24460, MPFR_RNDN); mpfr_mul(r24462, r24448, r24461, MPFR_RNDN); return mpfr_get_d(r24462, MPFR_RNDN); } static mpfr_t r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24463, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24464, "3", 10, MPFR_RNDN); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); mpfr_init(r24469); mpfr_init(r24470); mpfr_init(r24471); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init(r24477); } double f_dm(double g, double h) { ; ; mpfr_div(r24465, r24463, r24464, MPFR_RNDN); mpfr_const_pi(r24466, MPFR_RNDN); mpfr_set_d(r24467, g, MPFR_RNDN); mpfr_set_d(r24468, h, MPFR_RNDN); mpfr_div(r24469, r24467, r24468, MPFR_RNDN); mpfr_neg(r24470, r24469, MPFR_RNDN); mpfr_acos(r24471, r24470, MPFR_RNDN); mpfr_div(r24472, r24471, r24464, MPFR_RNDN); mpfr_fma(r24473, r24465, r24466, r24472, MPFR_RNDN); mpfr_cos(r24474, r24473, MPFR_RNDN); mpfr_expm1(r24475, r24474, MPFR_RNDN); mpfr_log1p(r24476, r24475, MPFR_RNDN); mpfr_mul(r24477, r24463, r24476, MPFR_RNDN); return mpfr_get_d(r24477, MPFR_RNDN); }