#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r20933 = 1; float r20934 = 6; float r20935 = r20933 / r20934; float r20936 = -2; float r20937 = u1; float r20938 = log(r20937); float r20939 = r20936 * r20938; float r20940 = 0.5; float r20941 = pow(r20939, r20940); float r20942 = r20935 * r20941; float r20943 = 2; float r20944 = atan2(1.0, 0.0); float r20945 = r20943 * r20944; float r20946 = u2; float r20947 = r20945 * r20946; float r20948 = cos(r20947); float r20949 = r20942 * r20948; float r20950 = r20949 + r20940; return r20950; } double f_id(double u1, double u2) { double r20951 = 1; double r20952 = 6; double r20953 = r20951 / r20952; double r20954 = -2; double r20955 = u1; double r20956 = log(r20955); double r20957 = r20954 * r20956; double r20958 = 0.5; double r20959 = pow(r20957, r20958); double r20960 = r20953 * r20959; double r20961 = 2; double r20962 = atan2(1.0, 0.0); double r20963 = r20961 * r20962; double r20964 = u2; double r20965 = r20963 * r20964; double r20966 = cos(r20965); double r20967 = r20960 * r20966; double r20968 = r20967 + r20958; return r20968; } double f_of(float u1, float u2) { float r20969 = 1; float r20970 = 6; float r20971 = -2; float r20972 = u1; float r20973 = log(r20972); float r20974 = r20971 * r20973; float r20975 = 0.5; float r20976 = pow(r20974, r20975); float r20977 = r20970 / r20976; float r20978 = r20969 / r20977; float r20979 = u2; float r20980 = atan2(1.0, 0.0); float r20981 = r20980 + r20980; float r20982 = r20979 * r20981; float r20983 = cos(r20982); float r20984 = fma(r20978, r20983, r20975); return r20984; } double f_od(double u1, double u2) { double r20985 = 1; double r20986 = 6; double r20987 = -2; double r20988 = u1; double r20989 = log(r20988); double r20990 = r20987 * r20989; double r20991 = 0.5; double r20992 = pow(r20990, r20991); double r20993 = r20986 / r20992; double r20994 = r20985 / r20993; double r20995 = u2; double r20996 = atan2(1.0, 0.0); double r20997 = r20996 + r20996; double r20998 = r20995 * r20997; double r20999 = cos(r20998); double r21000 = fma(r20994, r20999, r20991); return r21000; } 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 r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21001, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21002, "6", 10, MPFR_RNDN); mpfr_init(r21003); mpfr_init_set_str(r21004, "-2", 10, MPFR_RNDN); mpfr_init(r21005); mpfr_init(r21006); mpfr_init(r21007); mpfr_init_set_str(r21008, "0.5", 10, MPFR_RNDN); mpfr_init(r21009); mpfr_init(r21010); mpfr_init_set_str(r21011, "2", 10, MPFR_RNDN); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init(r21015); mpfr_init(r21016); mpfr_init(r21017); mpfr_init(r21018); } double f_im(double u1, double u2) { ; ; mpfr_div(r21003, r21001, r21002, MPFR_RNDN); ; mpfr_set_d(r21005, u1, MPFR_RNDN); mpfr_log(r21006, r21005, MPFR_RNDN); mpfr_mul(r21007, r21004, r21006, MPFR_RNDN); ; mpfr_pow(r21009, r21007, r21008, MPFR_RNDN); mpfr_mul(r21010, r21003, r21009, MPFR_RNDN); ; mpfr_const_pi(r21012, MPFR_RNDN); mpfr_mul(r21013, r21011, r21012, MPFR_RNDN); mpfr_set_d(r21014, u2, MPFR_RNDN); mpfr_mul(r21015, r21013, r21014, MPFR_RNDN); mpfr_cos(r21016, r21015, MPFR_RNDN); mpfr_mul(r21017, r21010, r21016, MPFR_RNDN); mpfr_add(r21018, r21017, r21008, MPFR_RNDN); return mpfr_get_d(r21018, MPFR_RNDN); } static mpfr_t r21019, r21020, r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21019, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21020, "6", 10, MPFR_RNDN); mpfr_init_set_str(r21021, "-2", 10, MPFR_RNDN); mpfr_init(r21022); mpfr_init(r21023); mpfr_init(r21024); mpfr_init_set_str(r21025, "0.5", 10, MPFR_RNDN); mpfr_init(r21026); mpfr_init(r21027); mpfr_init(r21028); mpfr_init(r21029); mpfr_init(r21030); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); } double f_fm(double u1, double u2) { ; ; ; mpfr_set_d(r21022, u1, MPFR_RNDN); mpfr_log(r21023, r21022, MPFR_RNDN); mpfr_mul(r21024, r21021, r21023, MPFR_RNDN); ; mpfr_pow(r21026, r21024, r21025, MPFR_RNDN); mpfr_div(r21027, r21020, r21026, MPFR_RNDN); mpfr_div(r21028, r21019, r21027, MPFR_RNDN); mpfr_set_d(r21029, u2, MPFR_RNDN); mpfr_const_pi(r21030, MPFR_RNDN); mpfr_add(r21031, r21030, r21030, MPFR_RNDN); mpfr_mul(r21032, r21029, r21031, MPFR_RNDN); mpfr_cos(r21033, r21032, MPFR_RNDN); mpfr_fma(r21034, r21028, r21033, r21025, MPFR_RNDN); return mpfr_get_d(r21034, MPFR_RNDN); } static mpfr_t r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21035, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21036, "6", 10, MPFR_RNDN); mpfr_init_set_str(r21037, "-2", 10, MPFR_RNDN); mpfr_init(r21038); mpfr_init(r21039); mpfr_init(r21040); mpfr_init_set_str(r21041, "0.5", 10, MPFR_RNDN); mpfr_init(r21042); mpfr_init(r21043); mpfr_init(r21044); mpfr_init(r21045); mpfr_init(r21046); mpfr_init(r21047); mpfr_init(r21048); mpfr_init(r21049); mpfr_init(r21050); } double f_dm(double u1, double u2) { ; ; ; mpfr_set_d(r21038, u1, MPFR_RNDN); mpfr_log(r21039, r21038, MPFR_RNDN); mpfr_mul(r21040, r21037, r21039, MPFR_RNDN); ; mpfr_pow(r21042, r21040, r21041, MPFR_RNDN); mpfr_div(r21043, r21036, r21042, MPFR_RNDN); mpfr_div(r21044, r21035, r21043, MPFR_RNDN); mpfr_set_d(r21045, u2, MPFR_RNDN); mpfr_const_pi(r21046, MPFR_RNDN); mpfr_add(r21047, r21046, r21046, MPFR_RNDN); mpfr_mul(r21048, r21045, r21047, MPFR_RNDN); mpfr_cos(r21049, r21048, MPFR_RNDN); mpfr_fma(r21050, r21044, r21049, r21041, MPFR_RNDN); return mpfr_get_d(r21050, MPFR_RNDN); }