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