#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r27237 = 2; float r27238 = atan2(1.0, 0.0); float r27239 = r27237 * r27238; float r27240 = 3; float r27241 = r27239 / r27240; float r27242 = g; float r27243 = -r27242; float r27244 = h; float r27245 = r27243 / r27244; float r27246 = acos(r27245); float r27247 = r27246 / r27240; float r27248 = r27241 + r27247; float r27249 = cos(r27248); float r27250 = r27237 * r27249; return r27250; } double f_id(double g, double h) { double r27251 = 2; double r27252 = atan2(1.0, 0.0); double r27253 = r27251 * r27252; double r27254 = 3; double r27255 = r27253 / r27254; double r27256 = g; double r27257 = -r27256; double r27258 = h; double r27259 = r27257 / r27258; double r27260 = acos(r27259); double r27261 = r27260 / r27254; double r27262 = r27255 + r27261; double r27263 = cos(r27262); double r27264 = r27251 * r27263; return r27264; } double f_of(float g, float h) { float r27265 = 2; float r27266 = 3; float r27267 = r27265 / r27266; float r27268 = atan2(1.0, 0.0); float r27269 = g; float r27270 = h; float r27271 = r27269 / r27270; float r27272 = -r27271; float r27273 = acos(r27272); float r27274 = r27273 / r27266; float r27275 = fma(r27267, r27268, r27274); float r27276 = cos(r27275); float r27277 = r27276 * r27276; float r27278 = cbrt(r27277); float r27279 = cbrt(r27276); float r27280 = r27278 * r27279; float r27281 = r27265 * r27280; return r27281; } double f_od(double g, double h) { double r27282 = 2; double r27283 = 3; double r27284 = r27282 / r27283; double r27285 = atan2(1.0, 0.0); double r27286 = g; double r27287 = h; double r27288 = r27286 / r27287; double r27289 = -r27288; double r27290 = acos(r27289); double r27291 = r27290 / r27283; double r27292 = fma(r27284, r27285, r27291); double r27293 = cos(r27292); double r27294 = r27293 * r27293; double r27295 = cbrt(r27294); double r27296 = cbrt(r27293); double r27297 = r27295 * r27296; double r27298 = r27282 * r27297; return r27298; } 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 r27299, r27300, r27301, r27302, r27303, r27304, r27305, r27306, r27307, r27308, r27309, r27310, r27311, r27312; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27299, "2", 10, MPFR_RNDN); mpfr_init(r27300); mpfr_init(r27301); mpfr_init_set_str(r27302, "3", 10, MPFR_RNDN); mpfr_init(r27303); mpfr_init(r27304); mpfr_init(r27305); mpfr_init(r27306); mpfr_init(r27307); mpfr_init(r27308); mpfr_init(r27309); mpfr_init(r27310); mpfr_init(r27311); mpfr_init(r27312); } double f_im(double g, double h) { ; mpfr_const_pi(r27300, MPFR_RNDN); mpfr_mul(r27301, r27299, r27300, MPFR_RNDN); ; mpfr_div(r27303, r27301, r27302, MPFR_RNDN); mpfr_set_d(r27304, g, MPFR_RNDN); mpfr_neg(r27305, r27304, MPFR_RNDN); mpfr_set_d(r27306, h, MPFR_RNDN); mpfr_div(r27307, r27305, r27306, MPFR_RNDN); mpfr_acos(r27308, r27307, MPFR_RNDN); mpfr_div(r27309, r27308, r27302, MPFR_RNDN); mpfr_add(r27310, r27303, r27309, MPFR_RNDN); mpfr_cos(r27311, r27310, MPFR_RNDN); mpfr_mul(r27312, r27299, r27311, MPFR_RNDN); return mpfr_get_d(r27312, MPFR_RNDN); } static mpfr_t r27313, r27314, r27315, r27316, r27317, r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27313, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27314, "3", 10, MPFR_RNDN); mpfr_init(r27315); mpfr_init(r27316); mpfr_init(r27317); mpfr_init(r27318); mpfr_init(r27319); mpfr_init(r27320); mpfr_init(r27321); mpfr_init(r27322); mpfr_init(r27323); mpfr_init(r27324); mpfr_init(r27325); mpfr_init(r27326); mpfr_init(r27327); mpfr_init(r27328); mpfr_init(r27329); } double f_fm(double g, double h) { ; ; mpfr_div(r27315, r27313, r27314, MPFR_RNDN); mpfr_const_pi(r27316, MPFR_RNDN); mpfr_set_d(r27317, g, MPFR_RNDN); mpfr_set_d(r27318, h, MPFR_RNDN); mpfr_div(r27319, r27317, r27318, MPFR_RNDN); mpfr_neg(r27320, r27319, MPFR_RNDN); mpfr_acos(r27321, r27320, MPFR_RNDN); mpfr_div(r27322, r27321, r27314, MPFR_RNDN); mpfr_fma(r27323, r27315, r27316, r27322, MPFR_RNDN); mpfr_cos(r27324, r27323, MPFR_RNDN); mpfr_mul(r27325, r27324, r27324, MPFR_RNDN); mpfr_cbrt(r27326, r27325, MPFR_RNDN); mpfr_cbrt(r27327, r27324, MPFR_RNDN); mpfr_mul(r27328, r27326, r27327, MPFR_RNDN); mpfr_mul(r27329, r27313, r27328, MPFR_RNDN); return mpfr_get_d(r27329, MPFR_RNDN); } static mpfr_t r27330, r27331, r27332, r27333, r27334, r27335, r27336, r27337, r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27330, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27331, "3", 10, MPFR_RNDN); mpfr_init(r27332); mpfr_init(r27333); mpfr_init(r27334); mpfr_init(r27335); mpfr_init(r27336); mpfr_init(r27337); mpfr_init(r27338); mpfr_init(r27339); mpfr_init(r27340); mpfr_init(r27341); mpfr_init(r27342); mpfr_init(r27343); mpfr_init(r27344); mpfr_init(r27345); mpfr_init(r27346); } double f_dm(double g, double h) { ; ; mpfr_div(r27332, r27330, r27331, MPFR_RNDN); mpfr_const_pi(r27333, MPFR_RNDN); mpfr_set_d(r27334, g, MPFR_RNDN); mpfr_set_d(r27335, h, MPFR_RNDN); mpfr_div(r27336, r27334, r27335, MPFR_RNDN); mpfr_neg(r27337, r27336, MPFR_RNDN); mpfr_acos(r27338, r27337, MPFR_RNDN); mpfr_div(r27339, r27338, r27331, MPFR_RNDN); mpfr_fma(r27340, r27332, r27333, r27339, MPFR_RNDN); mpfr_cos(r27341, r27340, MPFR_RNDN); mpfr_mul(r27342, r27341, r27341, MPFR_RNDN); mpfr_cbrt(r27343, r27342, MPFR_RNDN); mpfr_cbrt(r27344, r27341, MPFR_RNDN); mpfr_mul(r27345, r27343, r27344, MPFR_RNDN); mpfr_mul(r27346, r27330, r27345, MPFR_RNDN); return mpfr_get_d(r27346, MPFR_RNDN); }