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