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