#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r24953 = 2; float r24954 = atan2(1.0, 0.0); float r24955 = r24953 * r24954; float r24956 = 3; float r24957 = r24955 / r24956; float r24958 = g; float r24959 = -r24958; float r24960 = h; float r24961 = r24959 / r24960; float r24962 = acos(r24961); float r24963 = r24962 / r24956; float r24964 = r24957 + r24963; float r24965 = cos(r24964); float r24966 = r24953 * r24965; return r24966; } double f_id(double g, double h) { double r24967 = 2; double r24968 = atan2(1.0, 0.0); double r24969 = r24967 * r24968; double r24970 = 3; double r24971 = r24969 / r24970; double r24972 = g; double r24973 = -r24972; double r24974 = h; double r24975 = r24973 / r24974; double r24976 = acos(r24975); double r24977 = r24976 / r24970; double r24978 = r24971 + r24977; double r24979 = cos(r24978); double r24980 = r24967 * r24979; return r24980; } double f_of(float g, float h) { float r24981 = 2; float r24982 = 3; float r24983 = r24981 / r24982; float r24984 = atan2(1.0, 0.0); float r24985 = g; float r24986 = h; float r24987 = r24985 / r24986; float r24988 = -r24987; float r24989 = acos(r24988); float r24990 = r24989 / r24982; float r24991 = fma(r24983, r24984, r24990); float r24992 = cos(r24991); float r24993 = r24992 * r24992; float r24994 = cbrt(r24993); float r24995 = cbrt(r24992); float r24996 = r24994 * r24995; float r24997 = r24981 * r24996; return r24997; } double f_od(double g, double h) { double r24998 = 2; double r24999 = 3; double r25000 = r24998 / r24999; double r25001 = atan2(1.0, 0.0); double r25002 = g; double r25003 = h; double r25004 = r25002 / r25003; double r25005 = -r25004; double r25006 = acos(r25005); double r25007 = r25006 / r24999; double r25008 = fma(r25000, r25001, r25007); double r25009 = cos(r25008); double r25010 = r25009 * r25009; double r25011 = cbrt(r25010); double r25012 = cbrt(r25009); double r25013 = r25011 * r25012; double r25014 = r24998 * r25013; return r25014; } 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 r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025, r25026, r25027, r25028; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r25015, "2", 10, MPFR_RNDN); mpfr_init(r25016); mpfr_init(r25017); mpfr_init_set_str(r25018, "3", 10, MPFR_RNDN); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); mpfr_init(r25023); mpfr_init(r25024); mpfr_init(r25025); mpfr_init(r25026); mpfr_init(r25027); mpfr_init(r25028); } double f_im(double g, double h) { ; mpfr_const_pi(r25016, MPFR_RNDN); mpfr_mul(r25017, r25015, r25016, MPFR_RNDN); ; mpfr_div(r25019, r25017, r25018, MPFR_RNDN); mpfr_set_d(r25020, g, MPFR_RNDN); mpfr_neg(r25021, r25020, MPFR_RNDN); mpfr_set_d(r25022, h, MPFR_RNDN); mpfr_div(r25023, r25021, r25022, MPFR_RNDN); mpfr_acos(r25024, r25023, MPFR_RNDN); mpfr_div(r25025, r25024, r25018, MPFR_RNDN); mpfr_add(r25026, r25019, r25025, MPFR_RNDN); mpfr_cos(r25027, r25026, MPFR_RNDN); mpfr_mul(r25028, r25015, r25027, MPFR_RNDN); return mpfr_get_d(r25028, MPFR_RNDN); } static mpfr_t r25029, r25030, r25031, r25032, r25033, r25034, r25035, r25036, r25037, r25038, r25039, r25040, r25041, r25042, r25043, r25044, r25045; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r25029, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25030, "3", 10, MPFR_RNDN); mpfr_init(r25031); mpfr_init(r25032); mpfr_init(r25033); mpfr_init(r25034); mpfr_init(r25035); mpfr_init(r25036); mpfr_init(r25037); mpfr_init(r25038); mpfr_init(r25039); mpfr_init(r25040); mpfr_init(r25041); mpfr_init(r25042); mpfr_init(r25043); mpfr_init(r25044); mpfr_init(r25045); } double f_fm(double g, double h) { ; ; mpfr_div(r25031, r25029, r25030, MPFR_RNDN); mpfr_const_pi(r25032, MPFR_RNDN); mpfr_set_d(r25033, g, MPFR_RNDN); mpfr_set_d(r25034, h, MPFR_RNDN); mpfr_div(r25035, r25033, r25034, MPFR_RNDN); mpfr_neg(r25036, r25035, MPFR_RNDN); mpfr_acos(r25037, r25036, MPFR_RNDN); mpfr_div(r25038, r25037, r25030, MPFR_RNDN); mpfr_fma(r25039, r25031, r25032, r25038, MPFR_RNDN); mpfr_cos(r25040, r25039, MPFR_RNDN); mpfr_mul(r25041, r25040, r25040, MPFR_RNDN); mpfr_cbrt(r25042, r25041, MPFR_RNDN); mpfr_cbrt(r25043, r25040, MPFR_RNDN); mpfr_mul(r25044, r25042, r25043, MPFR_RNDN); mpfr_mul(r25045, r25029, r25044, MPFR_RNDN); return mpfr_get_d(r25045, MPFR_RNDN); } static mpfr_t r25046, r25047, r25048, r25049, r25050, r25051, r25052, r25053, r25054, r25055, r25056, r25057, r25058, r25059, r25060, r25061, r25062; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r25046, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25047, "3", 10, MPFR_RNDN); mpfr_init(r25048); mpfr_init(r25049); mpfr_init(r25050); mpfr_init(r25051); mpfr_init(r25052); mpfr_init(r25053); mpfr_init(r25054); mpfr_init(r25055); mpfr_init(r25056); mpfr_init(r25057); mpfr_init(r25058); mpfr_init(r25059); mpfr_init(r25060); mpfr_init(r25061); mpfr_init(r25062); } double f_dm(double g, double h) { ; ; mpfr_div(r25048, r25046, r25047, MPFR_RNDN); mpfr_const_pi(r25049, MPFR_RNDN); mpfr_set_d(r25050, g, MPFR_RNDN); mpfr_set_d(r25051, h, MPFR_RNDN); mpfr_div(r25052, r25050, r25051, MPFR_RNDN); mpfr_neg(r25053, r25052, MPFR_RNDN); mpfr_acos(r25054, r25053, MPFR_RNDN); mpfr_div(r25055, r25054, r25047, MPFR_RNDN); mpfr_fma(r25056, r25048, r25049, r25055, MPFR_RNDN); mpfr_cos(r25057, r25056, MPFR_RNDN); mpfr_mul(r25058, r25057, r25057, MPFR_RNDN); mpfr_cbrt(r25059, r25058, MPFR_RNDN); mpfr_cbrt(r25060, r25057, MPFR_RNDN); mpfr_mul(r25061, r25059, r25060, MPFR_RNDN); mpfr_mul(r25062, r25046, r25061, MPFR_RNDN); return mpfr_get_d(r25062, MPFR_RNDN); }