#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r20703 = 1; float r20704 = 6; float r20705 = r20703 / r20704; float r20706 = -2; float r20707 = u1; float r20708 = log(r20707); float r20709 = r20706 * r20708; float r20710 = 0.5; float r20711 = pow(r20709, r20710); float r20712 = r20705 * r20711; float r20713 = 2; float r20714 = atan2(1.0, 0.0); float r20715 = r20713 * r20714; float r20716 = u2; float r20717 = r20715 * r20716; float r20718 = cos(r20717); float r20719 = r20712 * r20718; float r20720 = r20719 + r20710; return r20720; } double f_id(double u1, double u2) { double r20721 = 1; double r20722 = 6; double r20723 = r20721 / r20722; double r20724 = -2; double r20725 = u1; double r20726 = log(r20725); double r20727 = r20724 * r20726; double r20728 = 0.5; double r20729 = pow(r20727, r20728); double r20730 = r20723 * r20729; double r20731 = 2; double r20732 = atan2(1.0, 0.0); double r20733 = r20731 * r20732; double r20734 = u2; double r20735 = r20733 * r20734; double r20736 = cos(r20735); double r20737 = r20730 * r20736; double r20738 = r20737 + r20728; return r20738; } double f_of(float u1, float u2) { float r20739 = 1; float r20740 = 6; float r20741 = -2; float r20742 = u1; float r20743 = log(r20742); float r20744 = r20741 * r20743; float r20745 = 0.5; float r20746 = pow(r20744, r20745); float r20747 = r20740 / r20746; float r20748 = r20739 / r20747; float r20749 = u2; float r20750 = atan2(1.0, 0.0); float r20751 = r20750 + r20750; float r20752 = r20749 * r20751; float r20753 = cos(r20752); float r20754 = fma(r20748, r20753, r20745); return r20754; } double f_od(double u1, double u2) { double r20755 = 1; double r20756 = 6; double r20757 = -2; double r20758 = u1; double r20759 = log(r20758); double r20760 = r20757 * r20759; double r20761 = 0.5; double r20762 = pow(r20760, r20761); double r20763 = r20756 / r20762; double r20764 = r20755 / r20763; double r20765 = u2; double r20766 = atan2(1.0, 0.0); double r20767 = r20766 + r20766; double r20768 = r20765 * r20767; double r20769 = cos(r20768); double r20770 = fma(r20764, r20769, r20761); return r20770; } 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 r20771, r20772, r20773, r20774, r20775, r20776, r20777, r20778, r20779, r20780, r20781, r20782, r20783, r20784, r20785, r20786, r20787, r20788; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20771, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20772, "6", 10, MPFR_RNDN); mpfr_init(r20773); mpfr_init_set_str(r20774, "-2", 10, MPFR_RNDN); mpfr_init(r20775); mpfr_init(r20776); mpfr_init(r20777); mpfr_init_set_str(r20778, "0.5", 10, MPFR_RNDN); mpfr_init(r20779); mpfr_init(r20780); mpfr_init_set_str(r20781, "2", 10, MPFR_RNDN); mpfr_init(r20782); mpfr_init(r20783); mpfr_init(r20784); mpfr_init(r20785); mpfr_init(r20786); mpfr_init(r20787); mpfr_init(r20788); } double f_im(double u1, double u2) { ; ; mpfr_div(r20773, r20771, r20772, MPFR_RNDN); ; mpfr_set_d(r20775, u1, MPFR_RNDN); mpfr_log(r20776, r20775, MPFR_RNDN); mpfr_mul(r20777, r20774, r20776, MPFR_RNDN); ; mpfr_pow(r20779, r20777, r20778, MPFR_RNDN); mpfr_mul(r20780, r20773, r20779, MPFR_RNDN); ; mpfr_const_pi(r20782, MPFR_RNDN); mpfr_mul(r20783, r20781, r20782, MPFR_RNDN); mpfr_set_d(r20784, u2, MPFR_RNDN); mpfr_mul(r20785, r20783, r20784, MPFR_RNDN); mpfr_cos(r20786, r20785, MPFR_RNDN); mpfr_mul(r20787, r20780, r20786, MPFR_RNDN); mpfr_add(r20788, r20787, r20778, MPFR_RNDN); return mpfr_get_d(r20788, MPFR_RNDN); } static mpfr_t r20789, r20790, r20791, r20792, r20793, r20794, r20795, r20796, r20797, r20798, r20799, r20800, r20801, r20802, r20803, r20804; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20789, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20790, "6", 10, MPFR_RNDN); mpfr_init_set_str(r20791, "-2", 10, MPFR_RNDN); mpfr_init(r20792); mpfr_init(r20793); mpfr_init(r20794); mpfr_init_set_str(r20795, "0.5", 10, MPFR_RNDN); mpfr_init(r20796); mpfr_init(r20797); mpfr_init(r20798); mpfr_init(r20799); mpfr_init(r20800); mpfr_init(r20801); mpfr_init(r20802); mpfr_init(r20803); mpfr_init(r20804); } double f_fm(double u1, double u2) { ; ; ; mpfr_set_d(r20792, u1, MPFR_RNDN); mpfr_log(r20793, r20792, MPFR_RNDN); mpfr_mul(r20794, r20791, r20793, MPFR_RNDN); ; mpfr_pow(r20796, r20794, r20795, MPFR_RNDN); mpfr_div(r20797, r20790, r20796, MPFR_RNDN); mpfr_div(r20798, r20789, r20797, MPFR_RNDN); mpfr_set_d(r20799, u2, MPFR_RNDN); mpfr_const_pi(r20800, MPFR_RNDN); mpfr_add(r20801, r20800, r20800, MPFR_RNDN); mpfr_mul(r20802, r20799, r20801, MPFR_RNDN); mpfr_cos(r20803, r20802, MPFR_RNDN); mpfr_fma(r20804, r20798, r20803, r20795, MPFR_RNDN); return mpfr_get_d(r20804, MPFR_RNDN); } static mpfr_t r20805, r20806, r20807, r20808, r20809, r20810, r20811, r20812, r20813, r20814, r20815, r20816, r20817, r20818, r20819, r20820; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20805, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20806, "6", 10, MPFR_RNDN); mpfr_init_set_str(r20807, "-2", 10, MPFR_RNDN); mpfr_init(r20808); mpfr_init(r20809); mpfr_init(r20810); mpfr_init_set_str(r20811, "0.5", 10, MPFR_RNDN); mpfr_init(r20812); mpfr_init(r20813); mpfr_init(r20814); mpfr_init(r20815); mpfr_init(r20816); mpfr_init(r20817); mpfr_init(r20818); mpfr_init(r20819); mpfr_init(r20820); } double f_dm(double u1, double u2) { ; ; ; mpfr_set_d(r20808, u1, MPFR_RNDN); mpfr_log(r20809, r20808, MPFR_RNDN); mpfr_mul(r20810, r20807, r20809, MPFR_RNDN); ; mpfr_pow(r20812, r20810, r20811, MPFR_RNDN); mpfr_div(r20813, r20806, r20812, MPFR_RNDN); mpfr_div(r20814, r20805, r20813, MPFR_RNDN); mpfr_set_d(r20815, u2, MPFR_RNDN); mpfr_const_pi(r20816, MPFR_RNDN); mpfr_add(r20817, r20816, r20816, MPFR_RNDN); mpfr_mul(r20818, r20815, r20817, MPFR_RNDN); mpfr_cos(r20819, r20818, MPFR_RNDN); mpfr_fma(r20820, r20814, r20819, r20811, MPFR_RNDN); return mpfr_get_d(r20820, MPFR_RNDN); }