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