#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 = 1.0f; float r15255 = 6.0f; float r15256 = r15254 / r15255; float r15257 = -2.0f; float r15258 = u1; float r15259 = log(r15258); float r15260 = r15257 * r15259; float r15261 = pow(r15260, r15253); float r15262 = atan2(1.0, 0.0); float r15263 = u2; float r15264 = 2.0f; float r15265 = r15263 * r15264; float r15266 = r15262 * r15265; float r15267 = cos(r15266); float r15268 = r15254 / r15267; float r15269 = r15261 / r15268; float r15270 = exp(r15269); float r15271 = log(r15270); float r15272 = r15256 * r15271; float r15273 = r15253 + r15272; return r15273; } double f_od(double u1, double u2) { double r15274 = 0.5; double r15275 = 1.0; double r15276 = 6.0; double r15277 = r15275 / r15276; double r15278 = -2.0; double r15279 = u1; double r15280 = log(r15279); double r15281 = r15278 * r15280; double r15282 = pow(r15281, r15274); double r15283 = atan2(1.0, 0.0); double r15284 = u2; double r15285 = 2.0; double r15286 = r15284 * r15285; double r15287 = r15283 * r15286; double r15288 = cos(r15287); double r15289 = r15275 / r15288; double r15290 = r15282 / r15289; double r15291 = exp(r15290); double r15292 = log(r15291); double r15293 = r15277 * r15292; double r15294 = r15274 + r15293; return r15294; } 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 r15295, r15296, r15297, r15298, r15299, r15300, r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310, r15311, r15312; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r15295, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15296, "6", 10, MPFR_RNDN); mpfr_init(r15297); mpfr_init_set_str(r15298, "-2", 10, MPFR_RNDN); mpfr_init(r15299); mpfr_init(r15300); mpfr_init(r15301); mpfr_init_set_str(r15302, "0.5", 10, MPFR_RNDN); mpfr_init(r15303); mpfr_init(r15304); mpfr_init_set_str(r15305, "2", 10, MPFR_RNDN); mpfr_init(r15306); mpfr_init(r15307); mpfr_init(r15308); mpfr_init(r15309); mpfr_init(r15310); mpfr_init(r15311); mpfr_init(r15312); } double f_im(double u1, double u2) { ; ; mpfr_div(r15297, r15295, r15296, MPFR_RNDN); ; mpfr_set_d(r15299, u1, MPFR_RNDN); mpfr_log(r15300, r15299, MPFR_RNDN); mpfr_mul(r15301, r15298, r15300, MPFR_RNDN); ; mpfr_pow(r15303, r15301, r15302, MPFR_RNDN); mpfr_mul(r15304, r15297, r15303, MPFR_RNDN); ; mpfr_const_pi(r15306, MPFR_RNDN); mpfr_mul(r15307, r15305, r15306, MPFR_RNDN); mpfr_set_d(r15308, u2, MPFR_RNDN); mpfr_mul(r15309, r15307, r15308, MPFR_RNDN); mpfr_cos(r15310, r15309, MPFR_RNDN); mpfr_mul(r15311, r15304, r15310, MPFR_RNDN); mpfr_add(r15312, r15311, r15302, MPFR_RNDN); return mpfr_get_d(r15312, MPFR_RNDN); } static mpfr_t r15313, r15314, r15315, r15316, r15317, r15318, r15319, r15320, r15321, r15322, r15323, r15324, r15325, r15326, r15327, r15328, r15329, r15330, r15331, r15332, r15333; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r15313, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r15314, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15315, "6", 10, MPFR_RNDN); 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_set_str(r15324, "2", 10, MPFR_RNDN); mpfr_init(r15325); mpfr_init(r15326); mpfr_init(r15327); mpfr_init(r15328); mpfr_init(r15329); mpfr_init(r15330); mpfr_init(r15331); mpfr_init(r15332); mpfr_init(r15333); } double f_fm(double u1, double u2) { ; ; ; mpfr_div(r15316, r15314, r15315, MPFR_RNDN); ; mpfr_set_d(r15318, u1, MPFR_RNDN); mpfr_log(r15319, r15318, MPFR_RNDN); mpfr_mul(r15320, r15317, r15319, MPFR_RNDN); mpfr_pow(r15321, r15320, r15313, MPFR_RNDN); mpfr_const_pi(r15322, MPFR_RNDN); mpfr_set_d(r15323, u2, MPFR_RNDN); ; mpfr_mul(r15325, r15323, r15324, MPFR_RNDN); mpfr_mul(r15326, r15322, r15325, MPFR_RNDN); mpfr_cos(r15327, r15326, MPFR_RNDN); mpfr_div(r15328, r15314, r15327, MPFR_RNDN); mpfr_div(r15329, r15321, r15328, MPFR_RNDN); mpfr_exp(r15330, r15329, MPFR_RNDN); mpfr_log(r15331, r15330, MPFR_RNDN); mpfr_mul(r15332, r15316, r15331, MPFR_RNDN); mpfr_add(r15333, r15313, r15332, MPFR_RNDN); return mpfr_get_d(r15333, MPFR_RNDN); } static mpfr_t r15334, r15335, 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(r15334, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r15335, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15336, "6", 10, MPFR_RNDN); 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_set_str(r15345, "2", 10, MPFR_RNDN); mpfr_init(r15346); mpfr_init(r15347); 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(r15337, r15335, r15336, MPFR_RNDN); ; mpfr_set_d(r15339, u1, MPFR_RNDN); mpfr_log(r15340, r15339, MPFR_RNDN); mpfr_mul(r15341, r15338, r15340, MPFR_RNDN); mpfr_pow(r15342, r15341, r15334, MPFR_RNDN); mpfr_const_pi(r15343, MPFR_RNDN); mpfr_set_d(r15344, u2, MPFR_RNDN); ; mpfr_mul(r15346, r15344, r15345, MPFR_RNDN); mpfr_mul(r15347, r15343, r15346, MPFR_RNDN); mpfr_cos(r15348, r15347, MPFR_RNDN); mpfr_div(r15349, r15335, r15348, MPFR_RNDN); mpfr_div(r15350, r15342, r15349, MPFR_RNDN); mpfr_exp(r15351, r15350, MPFR_RNDN); mpfr_log(r15352, r15351, MPFR_RNDN); mpfr_mul(r15353, r15337, r15352, MPFR_RNDN); mpfr_add(r15354, r15334, r15353, MPFR_RNDN); return mpfr_get_d(r15354, MPFR_RNDN); }