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