#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r21950 = 1; float r21951 = 6; float r21952 = r21950 / r21951; float r21953 = -2; float r21954 = u1; float r21955 = log(r21954); float r21956 = r21953 * r21955; float r21957 = 0.5; float r21958 = pow(r21956, r21957); float r21959 = r21952 * r21958; float r21960 = 2; float r21961 = atan2(1.0, 0.0); float r21962 = r21960 * r21961; float r21963 = u2; float r21964 = r21962 * r21963; float r21965 = cos(r21964); float r21966 = r21959 * r21965; float r21967 = r21966 + r21957; return r21967; } double f_id(double u1, double u2) { double r21968 = 1; double r21969 = 6; double r21970 = r21968 / r21969; double r21971 = -2; double r21972 = u1; double r21973 = log(r21972); double r21974 = r21971 * r21973; double r21975 = 0.5; double r21976 = pow(r21974, r21975); double r21977 = r21970 * r21976; double r21978 = 2; double r21979 = atan2(1.0, 0.0); double r21980 = r21978 * r21979; double r21981 = u2; double r21982 = r21980 * r21981; double r21983 = cos(r21982); double r21984 = r21977 * r21983; double r21985 = r21984 + r21975; return r21985; } double f_of(float u1, float u2) { float r21986 = 0.5; float r21987 = atan2(1.0, 0.0); float r21988 = u2; float r21989 = r21987 * r21988; float r21990 = 2; float r21991 = r21989 * r21990; float r21992 = cos(r21991); float r21993 = -2; float r21994 = u1; float r21995 = log(r21994); float r21996 = r21993 * r21995; float r21997 = pow(r21996, r21986); float r21998 = 6; float r21999 = r21997 / r21998; float r22000 = r21992 * r21999; float r22001 = r21986 + r22000; return r22001; } double f_od(double u1, double u2) { double r22002 = 0.5; double r22003 = atan2(1.0, 0.0); double r22004 = u2; double r22005 = r22003 * r22004; double r22006 = 2; double r22007 = r22005 * r22006; double r22008 = cos(r22007); double r22009 = -2; double r22010 = u1; double r22011 = log(r22010); double r22012 = r22009 * r22011; double r22013 = pow(r22012, r22002); double r22014 = 6; double r22015 = r22013 / r22014; double r22016 = r22008 * r22015; double r22017 = r22002 + r22016; return r22017; } 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 r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22018, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22019, "6", 10, MPFR_RNDN); mpfr_init(r22020); mpfr_init_set_str(r22021, "-2", 10, MPFR_RNDN); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init_set_str(r22025, "0.5", 10, MPFR_RNDN); mpfr_init(r22026); mpfr_init(r22027); mpfr_init_set_str(r22028, "2", 10, MPFR_RNDN); mpfr_init(r22029); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); } double f_im(double u1, double u2) { ; ; mpfr_div(r22020, r22018, r22019, MPFR_RNDN); ; mpfr_set_d(r22022, u1, MPFR_RNDN); mpfr_log(r22023, r22022, MPFR_RNDN); mpfr_mul(r22024, r22021, r22023, MPFR_RNDN); ; mpfr_pow(r22026, r22024, r22025, MPFR_RNDN); mpfr_mul(r22027, r22020, r22026, MPFR_RNDN); ; mpfr_const_pi(r22029, MPFR_RNDN); mpfr_mul(r22030, r22028, r22029, MPFR_RNDN); mpfr_set_d(r22031, u2, MPFR_RNDN); mpfr_mul(r22032, r22030, r22031, MPFR_RNDN); mpfr_cos(r22033, r22032, MPFR_RNDN); mpfr_mul(r22034, r22027, r22033, MPFR_RNDN); mpfr_add(r22035, r22034, r22025, MPFR_RNDN); return mpfr_get_d(r22035, MPFR_RNDN); } static mpfr_t r22036, r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22036, "0.5", 10, MPFR_RNDN); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); mpfr_init_set_str(r22040, "2", 10, MPFR_RNDN); mpfr_init(r22041); mpfr_init(r22042); mpfr_init_set_str(r22043, "-2", 10, MPFR_RNDN); mpfr_init(r22044); mpfr_init(r22045); mpfr_init(r22046); mpfr_init(r22047); mpfr_init_set_str(r22048, "6", 10, MPFR_RNDN); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); } double f_fm(double u1, double u2) { ; mpfr_const_pi(r22037, MPFR_RNDN); mpfr_set_d(r22038, u2, MPFR_RNDN); mpfr_mul(r22039, r22037, r22038, MPFR_RNDN); ; mpfr_mul(r22041, r22039, r22040, MPFR_RNDN); mpfr_cos(r22042, r22041, MPFR_RNDN); ; mpfr_set_d(r22044, u1, MPFR_RNDN); mpfr_log(r22045, r22044, MPFR_RNDN); mpfr_mul(r22046, r22043, r22045, MPFR_RNDN); mpfr_pow(r22047, r22046, r22036, MPFR_RNDN); ; mpfr_div(r22049, r22047, r22048, MPFR_RNDN); mpfr_mul(r22050, r22042, r22049, MPFR_RNDN); mpfr_add(r22051, r22036, r22050, MPFR_RNDN); return mpfr_get_d(r22051, MPFR_RNDN); } static mpfr_t r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22052, "0.5", 10, MPFR_RNDN); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init_set_str(r22056, "2", 10, MPFR_RNDN); mpfr_init(r22057); mpfr_init(r22058); mpfr_init_set_str(r22059, "-2", 10, MPFR_RNDN); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init_set_str(r22064, "6", 10, MPFR_RNDN); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); } double f_dm(double u1, double u2) { ; mpfr_const_pi(r22053, MPFR_RNDN); mpfr_set_d(r22054, u2, MPFR_RNDN); mpfr_mul(r22055, r22053, r22054, MPFR_RNDN); ; mpfr_mul(r22057, r22055, r22056, MPFR_RNDN); mpfr_cos(r22058, r22057, MPFR_RNDN); ; mpfr_set_d(r22060, u1, MPFR_RNDN); mpfr_log(r22061, r22060, MPFR_RNDN); mpfr_mul(r22062, r22059, r22061, MPFR_RNDN); mpfr_pow(r22063, r22062, r22052, MPFR_RNDN); ; mpfr_div(r22065, r22063, r22064, MPFR_RNDN); mpfr_mul(r22066, r22058, r22065, MPFR_RNDN); mpfr_add(r22067, r22052, r22066, MPFR_RNDN); return mpfr_get_d(r22067, MPFR_RNDN); }