#include #include #include #include #include char *name = "2-ancestry mixing, negative discriminant"; double f_if(float g, float h) { float r26909 = 2; float r26910 = atan2(1.0, 0.0); float r26911 = r26909 * r26910; float r26912 = 3; float r26913 = r26911 / r26912; float r26914 = g; float r26915 = -r26914; float r26916 = h; float r26917 = r26915 / r26916; float r26918 = acos(r26917); float r26919 = r26918 / r26912; float r26920 = r26913 + r26919; float r26921 = cos(r26920); float r26922 = r26909 * r26921; return r26922; } double f_id(double g, double h) { double r26923 = 2; double r26924 = atan2(1.0, 0.0); double r26925 = r26923 * r26924; double r26926 = 3; double r26927 = r26925 / r26926; double r26928 = g; double r26929 = -r26928; double r26930 = h; double r26931 = r26929 / r26930; double r26932 = acos(r26931); double r26933 = r26932 / r26926; double r26934 = r26927 + r26933; double r26935 = cos(r26934); double r26936 = r26923 * r26935; return r26936; } double f_of(float g, float h) { float r26937 = 2; float r26938 = 3; float r26939 = r26937 / r26938; float r26940 = atan2(1.0, 0.0); float r26941 = g; float r26942 = h; float r26943 = r26941 / r26942; float r26944 = -r26943; float r26945 = acos(r26944); float r26946 = sqrt(r26945); float r26947 = r26938 / r26946; float r26948 = r26946 / r26947; float r26949 = fma(r26939, r26940, r26948); float r26950 = cos(r26949); float r26951 = r26937 * r26950; return r26951; } double f_od(double g, double h) { double r26952 = 2; double r26953 = 3; double r26954 = r26952 / r26953; double r26955 = atan2(1.0, 0.0); double r26956 = g; double r26957 = h; double r26958 = r26956 / r26957; double r26959 = -r26958; double r26960 = acos(r26959); double r26961 = sqrt(r26960); double r26962 = r26953 / r26961; double r26963 = r26961 / r26962; double r26964 = fma(r26954, r26955, r26963); double r26965 = cos(r26964); double r26966 = r26952 * r26965; return r26966; } 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 r26967, r26968, r26969, r26970, r26971, r26972, r26973, r26974, r26975, r26976, r26977, r26978, r26979, r26980; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r26967, "2", 10, MPFR_RNDN); mpfr_init(r26968); mpfr_init(r26969); mpfr_init_set_str(r26970, "3", 10, MPFR_RNDN); mpfr_init(r26971); mpfr_init(r26972); mpfr_init(r26973); mpfr_init(r26974); mpfr_init(r26975); mpfr_init(r26976); mpfr_init(r26977); mpfr_init(r26978); mpfr_init(r26979); mpfr_init(r26980); } double f_im(double g, double h) { ; mpfr_const_pi(r26968, MPFR_RNDN); mpfr_mul(r26969, r26967, r26968, MPFR_RNDN); ; mpfr_div(r26971, r26969, r26970, MPFR_RNDN); mpfr_set_d(r26972, g, MPFR_RNDN); mpfr_neg(r26973, r26972, MPFR_RNDN); mpfr_set_d(r26974, h, MPFR_RNDN); mpfr_div(r26975, r26973, r26974, MPFR_RNDN); mpfr_acos(r26976, r26975, MPFR_RNDN); mpfr_div(r26977, r26976, r26970, MPFR_RNDN); mpfr_add(r26978, r26971, r26977, MPFR_RNDN); mpfr_cos(r26979, r26978, MPFR_RNDN); mpfr_mul(r26980, r26967, r26979, MPFR_RNDN); return mpfr_get_d(r26980, MPFR_RNDN); } static mpfr_t r26981, r26982, r26983, r26984, r26985, r26986, r26987, r26988, r26989, r26990, r26991, r26992, r26993, r26994, r26995; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26981, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26982, "3", 10, MPFR_RNDN); mpfr_init(r26983); mpfr_init(r26984); mpfr_init(r26985); mpfr_init(r26986); mpfr_init(r26987); mpfr_init(r26988); mpfr_init(r26989); mpfr_init(r26990); mpfr_init(r26991); mpfr_init(r26992); mpfr_init(r26993); mpfr_init(r26994); mpfr_init(r26995); } double f_fm(double g, double h) { ; ; mpfr_div(r26983, r26981, r26982, MPFR_RNDN); mpfr_const_pi(r26984, MPFR_RNDN); mpfr_set_d(r26985, g, MPFR_RNDN); mpfr_set_d(r26986, h, MPFR_RNDN); mpfr_div(r26987, r26985, r26986, MPFR_RNDN); mpfr_neg(r26988, r26987, MPFR_RNDN); mpfr_acos(r26989, r26988, MPFR_RNDN); mpfr_sqrt(r26990, r26989, MPFR_RNDN); mpfr_div(r26991, r26982, r26990, MPFR_RNDN); mpfr_div(r26992, r26990, r26991, MPFR_RNDN); mpfr_fma(r26993, r26983, r26984, r26992, MPFR_RNDN); mpfr_cos(r26994, r26993, MPFR_RNDN); mpfr_mul(r26995, r26981, r26994, MPFR_RNDN); return mpfr_get_d(r26995, MPFR_RNDN); } static mpfr_t r26996, r26997, r26998, r26999, r27000, r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26996, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26997, "3", 10, MPFR_RNDN); mpfr_init(r26998); mpfr_init(r26999); mpfr_init(r27000); mpfr_init(r27001); mpfr_init(r27002); mpfr_init(r27003); mpfr_init(r27004); mpfr_init(r27005); mpfr_init(r27006); mpfr_init(r27007); mpfr_init(r27008); mpfr_init(r27009); mpfr_init(r27010); } double f_dm(double g, double h) { ; ; mpfr_div(r26998, r26996, r26997, MPFR_RNDN); mpfr_const_pi(r26999, MPFR_RNDN); mpfr_set_d(r27000, g, MPFR_RNDN); mpfr_set_d(r27001, h, MPFR_RNDN); mpfr_div(r27002, r27000, r27001, MPFR_RNDN); mpfr_neg(r27003, r27002, MPFR_RNDN); mpfr_acos(r27004, r27003, MPFR_RNDN); mpfr_sqrt(r27005, r27004, MPFR_RNDN); mpfr_div(r27006, r26997, r27005, MPFR_RNDN); mpfr_div(r27007, r27005, r27006, MPFR_RNDN); mpfr_fma(r27008, r26998, r26999, r27007, MPFR_RNDN); mpfr_cos(r27009, r27008, MPFR_RNDN); mpfr_mul(r27010, r26996, r27009, MPFR_RNDN); return mpfr_get_d(r27010, MPFR_RNDN); }