#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r27026 = 1; float r27027 = 2; float r27028 = a; float r27029 = r27027 * r27028; float r27030 = r27026 / r27029; float r27031 = g; float r27032 = -r27031; float r27033 = r27031 * r27031; float r27034 = h; float r27035 = r27034 * r27034; float r27036 = r27033 - r27035; float r27037 = sqrt(r27036); float r27038 = r27032 + r27037; float r27039 = r27030 * r27038; float r27040 = cbrt(r27039); float r27041 = r27032 - r27037; float r27042 = r27030 * r27041; float r27043 = cbrt(r27042); float r27044 = r27040 + r27043; return r27044; } double f_id(double g, double h, double a) { double r27045 = 1; double r27046 = 2; double r27047 = a; double r27048 = r27046 * r27047; double r27049 = r27045 / r27048; double r27050 = g; double r27051 = -r27050; double r27052 = r27050 * r27050; double r27053 = h; double r27054 = r27053 * r27053; double r27055 = r27052 - r27054; double r27056 = sqrt(r27055); double r27057 = r27051 + r27056; double r27058 = r27049 * r27057; double r27059 = cbrt(r27058); double r27060 = r27051 - r27056; double r27061 = r27049 * r27060; double r27062 = cbrt(r27061); double r27063 = r27059 + r27062; return r27063; } double f_of(float g, float h, float a) { float r27064 = g; float r27065 = h; float r27066 = r27064 + r27065; float r27067 = r27064 - r27065; float r27068 = r27066 * r27067; float r27069 = sqrt(r27068); float r27070 = r27069 - r27064; float r27071 = cbrt(r27070); float r27072 = r27071 * r27071; float r27073 = r27072 * r27071; float r27074 = a; float r27075 = 2; float r27076 = r27074 * r27075; float r27077 = r27073 / r27076; float r27078 = cbrt(r27077); float r27079 = -r27064; float r27080 = r27079 - r27069; float r27081 = r27080 / r27076; float r27082 = cbrt(r27081); float r27083 = r27078 + r27082; return r27083; } double f_od(double g, double h, double a) { double r27084 = g; double r27085 = h; double r27086 = r27084 + r27085; double r27087 = r27084 - r27085; double r27088 = r27086 * r27087; double r27089 = sqrt(r27088); double r27090 = r27089 - r27084; double r27091 = cbrt(r27090); double r27092 = r27091 * r27091; double r27093 = r27092 * r27091; double r27094 = a; double r27095 = 2; double r27096 = r27094 * r27095; double r27097 = r27093 / r27096; double r27098 = cbrt(r27097); double r27099 = -r27084; double r27100 = r27099 - r27089; double r27101 = r27100 / r27096; double r27102 = cbrt(r27101); double r27103 = r27098 + r27102; return r27103; } 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 r27104, r27105, r27106, r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120, r27121, r27122; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27104, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27105, "2", 10, MPFR_RNDN); mpfr_init(r27106); mpfr_init(r27107); mpfr_init(r27108); mpfr_init(r27109); mpfr_init(r27110); mpfr_init(r27111); 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); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r27106, a, MPFR_RNDN); mpfr_mul(r27107, r27105, r27106, MPFR_RNDN); mpfr_div(r27108, r27104, r27107, MPFR_RNDN); mpfr_set_d(r27109, g, MPFR_RNDN); mpfr_neg(r27110, r27109, MPFR_RNDN); mpfr_mul(r27111, r27109, r27109, MPFR_RNDN); mpfr_set_d(r27112, h, MPFR_RNDN); mpfr_mul(r27113, r27112, r27112, MPFR_RNDN); mpfr_sub(r27114, r27111, r27113, MPFR_RNDN); mpfr_sqrt(r27115, r27114, MPFR_RNDN); mpfr_add(r27116, r27110, r27115, MPFR_RNDN); mpfr_mul(r27117, r27108, r27116, MPFR_RNDN); mpfr_cbrt(r27118, r27117, MPFR_RNDN); mpfr_sub(r27119, r27110, r27115, MPFR_RNDN); mpfr_mul(r27120, r27108, r27119, MPFR_RNDN); mpfr_cbrt(r27121, r27120, MPFR_RNDN); mpfr_add(r27122, r27118, r27121, MPFR_RNDN); return mpfr_get_d(r27122, MPFR_RNDN); } static mpfr_t r27123, r27124, r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135, r27136, r27137, r27138, r27139, r27140, r27141, r27142; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27123); mpfr_init(r27124); mpfr_init(r27125); mpfr_init(r27126); mpfr_init(r27127); mpfr_init(r27128); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init_set_str(r27134, "2", 10, MPFR_RNDN); mpfr_init(r27135); mpfr_init(r27136); mpfr_init(r27137); mpfr_init(r27138); mpfr_init(r27139); mpfr_init(r27140); mpfr_init(r27141); mpfr_init(r27142); } double f_fm(double g, double h, double a) { mpfr_set_d(r27123, g, MPFR_RNDN); mpfr_set_d(r27124, h, MPFR_RNDN); mpfr_add(r27125, r27123, r27124, MPFR_RNDN); mpfr_sub(r27126, r27123, r27124, MPFR_RNDN); mpfr_mul(r27127, r27125, r27126, MPFR_RNDN); mpfr_sqrt(r27128, r27127, MPFR_RNDN); mpfr_sub(r27129, r27128, r27123, MPFR_RNDN); mpfr_cbrt(r27130, r27129, MPFR_RNDN); mpfr_mul(r27131, r27130, r27130, MPFR_RNDN); mpfr_mul(r27132, r27131, r27130, MPFR_RNDN); mpfr_set_d(r27133, a, MPFR_RNDN); ; mpfr_mul(r27135, r27133, r27134, MPFR_RNDN); mpfr_div(r27136, r27132, r27135, MPFR_RNDN); mpfr_cbrt(r27137, r27136, MPFR_RNDN); mpfr_neg(r27138, r27123, MPFR_RNDN); mpfr_sub(r27139, r27138, r27128, MPFR_RNDN); mpfr_div(r27140, r27139, r27135, MPFR_RNDN); mpfr_cbrt(r27141, r27140, MPFR_RNDN); mpfr_add(r27142, r27137, r27141, MPFR_RNDN); return mpfr_get_d(r27142, MPFR_RNDN); } static mpfr_t r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161, r27162; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27143); mpfr_init(r27144); mpfr_init(r27145); mpfr_init(r27146); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init(r27151); mpfr_init(r27152); mpfr_init(r27153); mpfr_init_set_str(r27154, "2", 10, MPFR_RNDN); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); mpfr_init(r27158); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); mpfr_init(r27162); } double f_dm(double g, double h, double a) { mpfr_set_d(r27143, g, MPFR_RNDN); mpfr_set_d(r27144, h, MPFR_RNDN); mpfr_add(r27145, r27143, r27144, MPFR_RNDN); mpfr_sub(r27146, r27143, r27144, MPFR_RNDN); mpfr_mul(r27147, r27145, r27146, MPFR_RNDN); mpfr_sqrt(r27148, r27147, MPFR_RNDN); mpfr_sub(r27149, r27148, r27143, MPFR_RNDN); mpfr_cbrt(r27150, r27149, MPFR_RNDN); mpfr_mul(r27151, r27150, r27150, MPFR_RNDN); mpfr_mul(r27152, r27151, r27150, MPFR_RNDN); mpfr_set_d(r27153, a, MPFR_RNDN); ; mpfr_mul(r27155, r27153, r27154, MPFR_RNDN); mpfr_div(r27156, r27152, r27155, MPFR_RNDN); mpfr_cbrt(r27157, r27156, MPFR_RNDN); mpfr_neg(r27158, r27143, MPFR_RNDN); mpfr_sub(r27159, r27158, r27148, MPFR_RNDN); mpfr_div(r27160, r27159, r27155, MPFR_RNDN); mpfr_cbrt(r27161, r27160, MPFR_RNDN); mpfr_add(r27162, r27157, r27161, MPFR_RNDN); return mpfr_get_d(r27162, MPFR_RNDN); }