#include #include #include #include #include char *name = "2-ancestry mixing, positive discriminant"; double f_if(float g, float h, float a) { float r20950 = 1; float r20951 = 2; float r20952 = a; float r20953 = r20951 * r20952; float r20954 = r20950 / r20953; float r20955 = g; float r20956 = -r20955; float r20957 = r20955 * r20955; float r20958 = h; float r20959 = r20958 * r20958; float r20960 = r20957 - r20959; float r20961 = sqrt(r20960); float r20962 = r20956 + r20961; float r20963 = r20954 * r20962; float r20964 = cbrt(r20963); float r20965 = r20956 - r20961; float r20966 = r20954 * r20965; float r20967 = cbrt(r20966); float r20968 = r20964 + r20967; return r20968; } double f_id(double g, double h, double a) { double r20969 = 1; double r20970 = 2; double r20971 = a; double r20972 = r20970 * r20971; double r20973 = r20969 / r20972; double r20974 = g; double r20975 = -r20974; double r20976 = r20974 * r20974; double r20977 = h; double r20978 = r20977 * r20977; double r20979 = r20976 - r20978; double r20980 = sqrt(r20979); double r20981 = r20975 + r20980; double r20982 = r20973 * r20981; double r20983 = cbrt(r20982); double r20984 = r20975 - r20980; double r20985 = r20973 * r20984; double r20986 = cbrt(r20985); double r20987 = r20983 + r20986; return r20987; } double f_of(float g, float h, float a) { float r20988 = g; float r20989 = 1.5437032229238084e-162; bool r20990 = r20988 <= r20989; float r20991 = 0; float r20992 = cbrt(r20991); float r20993 = -r20988; float r20994 = h; float r20995 = r20988 + r20994; float r20996 = r20988 - r20994; float r20997 = r20995 * r20996; float r20998 = sqrt(r20997); float r20999 = r20993 + r20998; float r21000 = cbrt(r20999); float r21001 = a; float r21002 = r21001 + r21001; float r21003 = cbrt(r21002); float r21004 = r21000 / r21003; float r21005 = r20992 + r21004; float r21006 = r20999 / r21002; float r21007 = cbrt(r21006); float r21008 = r20993 - r20998; float r21009 = cbrt(r21008); float r21010 = 1; float r21011 = r21010 / r21002; float r21012 = cbrt(r21011); float r21013 = r21009 * r21012; float r21014 = r21007 + r21013; float r21015 = r20990 ? r21005 : r21014; return r21015; } double f_od(double g, double h, double a) { double r21016 = g; double r21017 = 1.5437032229238084e-162; bool r21018 = r21016 <= r21017; double r21019 = 0; double r21020 = cbrt(r21019); double r21021 = -r21016; double r21022 = h; double r21023 = r21016 + r21022; double r21024 = r21016 - r21022; double r21025 = r21023 * r21024; double r21026 = sqrt(r21025); double r21027 = r21021 + r21026; double r21028 = cbrt(r21027); double r21029 = a; double r21030 = r21029 + r21029; double r21031 = cbrt(r21030); double r21032 = r21028 / r21031; double r21033 = r21020 + r21032; double r21034 = r21027 / r21030; double r21035 = cbrt(r21034); double r21036 = r21021 - r21026; double r21037 = cbrt(r21036); double r21038 = 1; double r21039 = r21038 / r21030; double r21040 = cbrt(r21039); double r21041 = r21037 * r21040; double r21042 = r21035 + r21041; double r21043 = r21018 ? r21033 : r21042; return r21043; } 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 r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055, r21056, r21057, r21058, r21059, r21060, r21061, r21062; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init_set_str(r21044, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21045, "2", 10, MPFR_RNDN); mpfr_init(r21046); mpfr_init(r21047); mpfr_init(r21048); mpfr_init(r21049); mpfr_init(r21050); mpfr_init(r21051); mpfr_init(r21052); mpfr_init(r21053); mpfr_init(r21054); mpfr_init(r21055); mpfr_init(r21056); mpfr_init(r21057); mpfr_init(r21058); mpfr_init(r21059); mpfr_init(r21060); mpfr_init(r21061); mpfr_init(r21062); } double f_im(double g, double h, double a) { ; ; mpfr_set_d(r21046, a, MPFR_RNDN); mpfr_mul(r21047, r21045, r21046, MPFR_RNDN); mpfr_div(r21048, r21044, r21047, MPFR_RNDN); mpfr_set_d(r21049, g, MPFR_RNDN); mpfr_neg(r21050, r21049, MPFR_RNDN); mpfr_mul(r21051, r21049, r21049, MPFR_RNDN); mpfr_set_d(r21052, h, MPFR_RNDN); mpfr_mul(r21053, r21052, r21052, MPFR_RNDN); mpfr_sub(r21054, r21051, r21053, MPFR_RNDN); mpfr_sqrt(r21055, r21054, MPFR_RNDN); mpfr_add(r21056, r21050, r21055, MPFR_RNDN); mpfr_mul(r21057, r21048, r21056, MPFR_RNDN); mpfr_cbrt(r21058, r21057, MPFR_RNDN); mpfr_sub(r21059, r21050, r21055, MPFR_RNDN); mpfr_mul(r21060, r21048, r21059, MPFR_RNDN); mpfr_cbrt(r21061, r21060, MPFR_RNDN); mpfr_add(r21062, r21058, r21061, MPFR_RNDN); return mpfr_get_d(r21062, MPFR_RNDN); } static mpfr_t r21063, r21064, r21065, r21066, r21067, r21068, r21069, r21070, r21071, r21072, r21073, r21074, r21075, r21076, r21077, r21078, r21079, r21080, r21081, r21082, r21083, r21084, r21085, r21086, r21087, r21088, r21089, r21090; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r21063); mpfr_init_set_str(r21064, "1.5437032229238084e-162", 10, MPFR_RNDN); mpfr_init(r21065); mpfr_init_set_str(r21066, "0", 10, MPFR_RNDN); mpfr_init(r21067); mpfr_init(r21068); mpfr_init(r21069); mpfr_init(r21070); mpfr_init(r21071); mpfr_init(r21072); mpfr_init(r21073); mpfr_init(r21074); mpfr_init(r21075); mpfr_init(r21076); mpfr_init(r21077); mpfr_init(r21078); mpfr_init(r21079); mpfr_init(r21080); mpfr_init(r21081); mpfr_init(r21082); mpfr_init(r21083); mpfr_init(r21084); mpfr_init_set_str(r21085, "1", 10, MPFR_RNDN); mpfr_init(r21086); mpfr_init(r21087); mpfr_init(r21088); mpfr_init(r21089); mpfr_init(r21090); } double f_fm(double g, double h, double a) { mpfr_set_d(r21063, g, MPFR_RNDN); ; mpfr_set_si(r21065, mpfr_cmp(r21063, r21064) <= 0, MPFR_RNDN); ; mpfr_cbrt(r21067, r21066, MPFR_RNDN); mpfr_neg(r21068, r21063, MPFR_RNDN); mpfr_set_d(r21069, h, MPFR_RNDN); mpfr_add(r21070, r21063, r21069, MPFR_RNDN); mpfr_sub(r21071, r21063, r21069, MPFR_RNDN); mpfr_mul(r21072, r21070, r21071, MPFR_RNDN); mpfr_sqrt(r21073, r21072, MPFR_RNDN); mpfr_add(r21074, r21068, r21073, MPFR_RNDN); mpfr_cbrt(r21075, r21074, MPFR_RNDN); mpfr_set_d(r21076, a, MPFR_RNDN); mpfr_add(r21077, r21076, r21076, MPFR_RNDN); mpfr_cbrt(r21078, r21077, MPFR_RNDN); mpfr_div(r21079, r21075, r21078, MPFR_RNDN); mpfr_add(r21080, r21067, r21079, MPFR_RNDN); mpfr_div(r21081, r21074, r21077, MPFR_RNDN); mpfr_cbrt(r21082, r21081, MPFR_RNDN); mpfr_sub(r21083, r21068, r21073, MPFR_RNDN); mpfr_cbrt(r21084, r21083, MPFR_RNDN); ; mpfr_div(r21086, r21085, r21077, MPFR_RNDN); mpfr_cbrt(r21087, r21086, MPFR_RNDN); mpfr_mul(r21088, r21084, r21087, MPFR_RNDN); mpfr_add(r21089, r21082, r21088, MPFR_RNDN); if (mpfr_get_si(r21065, MPFR_RNDN)) { mpfr_set(r21090, r21080, MPFR_RNDN); } else { mpfr_set(r21090, r21089, MPFR_RNDN); }; return mpfr_get_d(r21090, MPFR_RNDN); } static mpfr_t r21091, r21092, r21093, r21094, r21095, r21096, r21097, r21098, r21099, r21100, r21101, r21102, r21103, r21104, r21105, r21106, r21107, r21108, r21109, r21110, r21111, r21112, r21113, r21114, r21115, r21116, r21117, r21118; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r21091); mpfr_init_set_str(r21092, "1.5437032229238084e-162", 10, MPFR_RNDN); mpfr_init(r21093); mpfr_init_set_str(r21094, "0", 10, MPFR_RNDN); mpfr_init(r21095); mpfr_init(r21096); mpfr_init(r21097); mpfr_init(r21098); mpfr_init(r21099); mpfr_init(r21100); mpfr_init(r21101); mpfr_init(r21102); mpfr_init(r21103); mpfr_init(r21104); mpfr_init(r21105); mpfr_init(r21106); mpfr_init(r21107); mpfr_init(r21108); mpfr_init(r21109); mpfr_init(r21110); mpfr_init(r21111); mpfr_init(r21112); mpfr_init_set_str(r21113, "1", 10, MPFR_RNDN); mpfr_init(r21114); mpfr_init(r21115); mpfr_init(r21116); mpfr_init(r21117); mpfr_init(r21118); } double f_dm(double g, double h, double a) { mpfr_set_d(r21091, g, MPFR_RNDN); ; mpfr_set_si(r21093, mpfr_cmp(r21091, r21092) <= 0, MPFR_RNDN); ; mpfr_cbrt(r21095, r21094, MPFR_RNDN); mpfr_neg(r21096, r21091, MPFR_RNDN); mpfr_set_d(r21097, h, MPFR_RNDN); mpfr_add(r21098, r21091, r21097, MPFR_RNDN); mpfr_sub(r21099, r21091, r21097, MPFR_RNDN); mpfr_mul(r21100, r21098, r21099, MPFR_RNDN); mpfr_sqrt(r21101, r21100, MPFR_RNDN); mpfr_add(r21102, r21096, r21101, MPFR_RNDN); mpfr_cbrt(r21103, r21102, MPFR_RNDN); mpfr_set_d(r21104, a, MPFR_RNDN); mpfr_add(r21105, r21104, r21104, MPFR_RNDN); mpfr_cbrt(r21106, r21105, MPFR_RNDN); mpfr_div(r21107, r21103, r21106, MPFR_RNDN); mpfr_add(r21108, r21095, r21107, MPFR_RNDN); mpfr_div(r21109, r21102, r21105, MPFR_RNDN); mpfr_cbrt(r21110, r21109, MPFR_RNDN); mpfr_sub(r21111, r21096, r21101, MPFR_RNDN); mpfr_cbrt(r21112, r21111, MPFR_RNDN); ; mpfr_div(r21114, r21113, r21105, MPFR_RNDN); mpfr_cbrt(r21115, r21114, MPFR_RNDN); mpfr_mul(r21116, r21112, r21115, MPFR_RNDN); mpfr_add(r21117, r21110, r21116, MPFR_RNDN); if (mpfr_get_si(r21093, MPFR_RNDN)) { mpfr_set(r21118, r21108, MPFR_RNDN); } else { mpfr_set(r21118, r21117, MPFR_RNDN); }; return mpfr_get_d(r21118, MPFR_RNDN); }