#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r27034 = 2; float r27035 = atan2(1.0, 0.0); float r27036 = r27034 * r27035; float r27037 = 3; float r27038 = r27036 / r27037; float r27039 = g; float r27040 = -r27039; float r27041 = h; float r27042 = r27040 / r27041; float r27043 = acos(r27042); float r27044 = r27043 / r27037; float r27045 = r27038 + r27044; float r27046 = cos(r27045); float r27047 = r27034 * r27046; return r27047; } double f_id(double g, double h) { double r27048 = 2; double r27049 = atan2(1.0, 0.0); double r27050 = r27048 * r27049; double r27051 = 3; double r27052 = r27050 / r27051; double r27053 = g; double r27054 = -r27053; double r27055 = h; double r27056 = r27054 / r27055; double r27057 = acos(r27056); double r27058 = r27057 / r27051; double r27059 = r27052 + r27058; double r27060 = cos(r27059); double r27061 = r27048 * r27060; return r27061; } double f_of(float g, float h) { float r27062 = 2; float r27063 = 3; float r27064 = r27062 / r27063; float r27065 = atan2(1.0, 0.0); float r27066 = g; float r27067 = h; float r27068 = r27066 / r27067; float r27069 = -r27068; float r27070 = acos(r27069); float r27071 = r27070 / r27063; float r27072 = fma(r27064, r27065, r27071); float r27073 = cos(r27072); float r27074 = r27073 * r27073; float r27075 = cbrt(r27074); float r27076 = cbrt(r27073); float r27077 = r27075 * r27076; float r27078 = r27062 * r27077; return r27078; } double f_od(double g, double h) { double r27079 = 2; double r27080 = 3; double r27081 = r27079 / r27080; double r27082 = atan2(1.0, 0.0); double r27083 = g; double r27084 = h; double r27085 = r27083 / r27084; double r27086 = -r27085; double r27087 = acos(r27086); double r27088 = r27087 / r27080; double r27089 = fma(r27081, r27082, r27088); double r27090 = cos(r27089); double r27091 = r27090 * r27090; double r27092 = cbrt(r27091); double r27093 = cbrt(r27090); double r27094 = r27092 * r27093; double r27095 = r27079 * r27094; return r27095; } 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 r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108, r27109; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27096, "2", 10, MPFR_RNDN); mpfr_init(r27097); mpfr_init(r27098); mpfr_init_set_str(r27099, "3", 10, MPFR_RNDN); mpfr_init(r27100); mpfr_init(r27101); mpfr_init(r27102); mpfr_init(r27103); mpfr_init(r27104); mpfr_init(r27105); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init(r27109); } double f_im(double g, double h) { ; mpfr_const_pi(r27097, MPFR_RNDN); mpfr_mul(r27098, r27096, r27097, MPFR_RNDN); ; mpfr_div(r27100, r27098, r27099, MPFR_RNDN); mpfr_set_d(r27101, g, MPFR_RNDN); mpfr_neg(r27102, r27101, MPFR_RNDN); mpfr_set_d(r27103, h, MPFR_RNDN); mpfr_div(r27104, r27102, r27103, MPFR_RNDN); mpfr_acos(r27105, r27104, MPFR_RNDN); mpfr_div(r27106, r27105, r27099, MPFR_RNDN); mpfr_add(r27107, r27100, r27106, MPFR_RNDN); mpfr_cos(r27108, r27107, MPFR_RNDN); mpfr_mul(r27109, r27096, r27108, MPFR_RNDN); return mpfr_get_d(r27109, MPFR_RNDN); } static mpfr_t r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120, r27121, r27122, r27123, r27124, r27125, r27126; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27110, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27111, "3", 10, MPFR_RNDN); mpfr_init(r27112); mpfr_init(r27113); mpfr_init(r27114); mpfr_init(r27115); mpfr_init(r27116); mpfr_init(r27117); mpfr_init(r27118); mpfr_init(r27119); mpfr_init(r27120); mpfr_init(r27121); mpfr_init(r27122); mpfr_init(r27123); mpfr_init(r27124); mpfr_init(r27125); mpfr_init(r27126); } double f_fm(double g, double h) { ; ; mpfr_div(r27112, r27110, r27111, MPFR_RNDN); mpfr_const_pi(r27113, MPFR_RNDN); mpfr_set_d(r27114, g, MPFR_RNDN); mpfr_set_d(r27115, h, MPFR_RNDN); mpfr_div(r27116, r27114, r27115, MPFR_RNDN); mpfr_neg(r27117, r27116, MPFR_RNDN); mpfr_acos(r27118, r27117, MPFR_RNDN); mpfr_div(r27119, r27118, r27111, MPFR_RNDN); mpfr_fma(r27120, r27112, r27113, r27119, MPFR_RNDN); mpfr_cos(r27121, r27120, MPFR_RNDN); mpfr_mul(r27122, r27121, r27121, MPFR_RNDN); mpfr_cbrt(r27123, r27122, MPFR_RNDN); mpfr_cbrt(r27124, r27121, MPFR_RNDN); mpfr_mul(r27125, r27123, r27124, MPFR_RNDN); mpfr_mul(r27126, r27110, r27125, MPFR_RNDN); return mpfr_get_d(r27126, MPFR_RNDN); } static mpfr_t r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142, r27143; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27127, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27128, "3", 10, MPFR_RNDN); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init(r27135); mpfr_init(r27136); mpfr_init(r27137); mpfr_init(r27138); mpfr_init(r27139); mpfr_init(r27140); mpfr_init(r27141); mpfr_init(r27142); mpfr_init(r27143); } double f_dm(double g, double h) { ; ; mpfr_div(r27129, r27127, r27128, MPFR_RNDN); mpfr_const_pi(r27130, MPFR_RNDN); mpfr_set_d(r27131, g, MPFR_RNDN); mpfr_set_d(r27132, h, MPFR_RNDN); mpfr_div(r27133, r27131, r27132, MPFR_RNDN); mpfr_neg(r27134, r27133, MPFR_RNDN); mpfr_acos(r27135, r27134, MPFR_RNDN); mpfr_div(r27136, r27135, r27128, MPFR_RNDN); mpfr_fma(r27137, r27129, r27130, r27136, MPFR_RNDN); mpfr_cos(r27138, r27137, MPFR_RNDN); mpfr_mul(r27139, r27138, r27138, MPFR_RNDN); mpfr_cbrt(r27140, r27139, MPFR_RNDN); mpfr_cbrt(r27141, r27138, MPFR_RNDN); mpfr_mul(r27142, r27140, r27141, MPFR_RNDN); mpfr_mul(r27143, r27127, r27142, MPFR_RNDN); return mpfr_get_d(r27143, MPFR_RNDN); }