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