#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r15225 = 1.0f; float r15226 = 6.0f; float r15227 = r15225 / r15226; float r15228 = -2.0f; float r15229 = u1; float r15230 = log(r15229); float r15231 = r15228 * r15230; float r15232 = 0.5f; float r15233 = pow(r15231, r15232); float r15234 = r15227 * r15233; float r15235 = 2.0f; float r15236 = atan2(1.0, 0.0); float r15237 = r15235 * r15236; float r15238 = u2; float r15239 = r15237 * r15238; float r15240 = cos(r15239); float r15241 = r15234 * r15240; float r15242 = r15241 + r15232; return r15242; } double f_id(double u1, double u2) { double r15243 = 1.0; double r15244 = 6.0; double r15245 = r15243 / r15244; double r15246 = -2.0; double r15247 = u1; double r15248 = log(r15247); double r15249 = r15246 * r15248; double r15250 = 0.5; double r15251 = pow(r15249, r15250); double r15252 = r15245 * r15251; double r15253 = 2.0; double r15254 = atan2(1.0, 0.0); double r15255 = r15253 * r15254; double r15256 = u2; double r15257 = r15255 * r15256; double r15258 = cos(r15257); double r15259 = r15252 * r15258; double r15260 = r15259 + r15250; return r15260; } double f_of(float u1, float u2) { float r15261 = 0.5f; float r15262 = 1.0f; float r15263 = 6.0f; float r15264 = r15262 / r15263; float r15265 = -2.0f; float r15266 = u1; float r15267 = log(r15266); float r15268 = r15265 * r15267; float r15269 = pow(r15268, r15261); float r15270 = atan2(1.0, 0.0); float r15271 = u2; float r15272 = 2.0f; float r15273 = r15271 * r15272; float r15274 = r15270 * r15273; float r15275 = cos(r15274); float r15276 = r15262 / r15275; float r15277 = r15269 / r15276; float r15278 = r15264 * r15277; float r15279 = r15261 + r15278; return r15279; } double f_od(double u1, double u2) { double r15280 = 0.5; double r15281 = 1.0; double r15282 = 6.0; double r15283 = r15281 / r15282; double r15284 = -2.0; double r15285 = u1; double r15286 = log(r15285); double r15287 = r15284 * r15286; double r15288 = pow(r15287, r15280); double r15289 = atan2(1.0, 0.0); double r15290 = u2; double r15291 = 2.0; double r15292 = r15290 * r15291; double r15293 = r15289 * r15292; double r15294 = cos(r15293); double r15295 = r15281 / r15294; double r15296 = r15288 / r15295; double r15297 = r15283 * r15296; double r15298 = r15280 + r15297; return r15298; } 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 r15299, r15300, r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310, r15311, r15312, r15313, r15314, r15315, r15316; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r15299, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15300, "6", 10, MPFR_RNDN); mpfr_init(r15301); mpfr_init_set_str(r15302, "-2", 10, MPFR_RNDN); mpfr_init(r15303); mpfr_init(r15304); mpfr_init(r15305); mpfr_init_set_str(r15306, "0.5", 10, MPFR_RNDN); mpfr_init(r15307); mpfr_init(r15308); mpfr_init_set_str(r15309, "2", 10, MPFR_RNDN); mpfr_init(r15310); mpfr_init(r15311); mpfr_init(r15312); mpfr_init(r15313); mpfr_init(r15314); mpfr_init(r15315); mpfr_init(r15316); } double f_im(double u1, double u2) { ; ; mpfr_div(r15301, r15299, r15300, MPFR_RNDN); ; mpfr_set_d(r15303, u1, MPFR_RNDN); mpfr_log(r15304, r15303, MPFR_RNDN); mpfr_mul(r15305, r15302, r15304, MPFR_RNDN); ; mpfr_pow(r15307, r15305, r15306, MPFR_RNDN); mpfr_mul(r15308, r15301, r15307, MPFR_RNDN); ; mpfr_const_pi(r15310, MPFR_RNDN); mpfr_mul(r15311, r15309, r15310, MPFR_RNDN); mpfr_set_d(r15312, u2, MPFR_RNDN); mpfr_mul(r15313, r15311, r15312, MPFR_RNDN); mpfr_cos(r15314, r15313, MPFR_RNDN); mpfr_mul(r15315, r15308, r15314, MPFR_RNDN); mpfr_add(r15316, r15315, r15306, MPFR_RNDN); return mpfr_get_d(r15316, MPFR_RNDN); } static mpfr_t r15317, r15318, r15319, r15320, r15321, r15322, r15323, r15324, r15325, r15326, r15327, r15328, r15329, r15330, r15331, r15332, r15333, r15334, r15335; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r15317, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r15318, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15319, "6", 10, MPFR_RNDN); mpfr_init(r15320); mpfr_init_set_str(r15321, "-2", 10, MPFR_RNDN); mpfr_init(r15322); mpfr_init(r15323); mpfr_init(r15324); mpfr_init(r15325); mpfr_init(r15326); mpfr_init(r15327); mpfr_init_set_str(r15328, "2", 10, MPFR_RNDN); mpfr_init(r15329); mpfr_init(r15330); mpfr_init(r15331); mpfr_init(r15332); mpfr_init(r15333); mpfr_init(r15334); mpfr_init(r15335); } double f_fm(double u1, double u2) { ; ; ; mpfr_div(r15320, r15318, r15319, MPFR_RNDN); ; mpfr_set_d(r15322, u1, MPFR_RNDN); mpfr_log(r15323, r15322, MPFR_RNDN); mpfr_mul(r15324, r15321, r15323, MPFR_RNDN); mpfr_pow(r15325, r15324, r15317, MPFR_RNDN); mpfr_const_pi(r15326, MPFR_RNDN); mpfr_set_d(r15327, u2, MPFR_RNDN); ; mpfr_mul(r15329, r15327, r15328, MPFR_RNDN); mpfr_mul(r15330, r15326, r15329, MPFR_RNDN); mpfr_cos(r15331, r15330, MPFR_RNDN); mpfr_div(r15332, r15318, r15331, MPFR_RNDN); mpfr_div(r15333, r15325, r15332, MPFR_RNDN); mpfr_mul(r15334, r15320, r15333, MPFR_RNDN); mpfr_add(r15335, r15317, r15334, MPFR_RNDN); return mpfr_get_d(r15335, MPFR_RNDN); } static mpfr_t r15336, r15337, r15338, r15339, r15340, r15341, r15342, r15343, r15344, r15345, r15346, r15347, r15348, r15349, r15350, r15351, r15352, r15353, r15354; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r15336, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r15337, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15338, "6", 10, MPFR_RNDN); mpfr_init(r15339); mpfr_init_set_str(r15340, "-2", 10, MPFR_RNDN); mpfr_init(r15341); mpfr_init(r15342); mpfr_init(r15343); mpfr_init(r15344); mpfr_init(r15345); mpfr_init(r15346); mpfr_init_set_str(r15347, "2", 10, MPFR_RNDN); mpfr_init(r15348); mpfr_init(r15349); mpfr_init(r15350); mpfr_init(r15351); mpfr_init(r15352); mpfr_init(r15353); mpfr_init(r15354); } double f_dm(double u1, double u2) { ; ; ; mpfr_div(r15339, r15337, r15338, MPFR_RNDN); ; mpfr_set_d(r15341, u1, MPFR_RNDN); mpfr_log(r15342, r15341, MPFR_RNDN); mpfr_mul(r15343, r15340, r15342, MPFR_RNDN); mpfr_pow(r15344, r15343, r15336, MPFR_RNDN); mpfr_const_pi(r15345, MPFR_RNDN); mpfr_set_d(r15346, u2, MPFR_RNDN); ; mpfr_mul(r15348, r15346, r15347, MPFR_RNDN); mpfr_mul(r15349, r15345, r15348, MPFR_RNDN); mpfr_cos(r15350, r15349, MPFR_RNDN); mpfr_div(r15351, r15337, r15350, MPFR_RNDN); mpfr_div(r15352, r15344, r15351, MPFR_RNDN); mpfr_mul(r15353, r15339, r15352, MPFR_RNDN); mpfr_add(r15354, r15336, r15353, MPFR_RNDN); return mpfr_get_d(r15354, MPFR_RNDN); }