#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r15201 = 1.0f; float r15202 = 6.0f; float r15203 = r15201 / r15202; float r15204 = -2.0f; float r15205 = u1; float r15206 = log(r15205); float r15207 = r15204 * r15206; float r15208 = 0.5f; float r15209 = pow(r15207, r15208); float r15210 = r15203 * r15209; float r15211 = 2.0f; float r15212 = atan2(1.0, 0.0); float r15213 = r15211 * r15212; float r15214 = u2; float r15215 = r15213 * r15214; float r15216 = cos(r15215); float r15217 = r15210 * r15216; float r15218 = r15217 + r15208; return r15218; } double f_id(double u1, double u2) { double r15219 = 1.0; double r15220 = 6.0; double r15221 = r15219 / r15220; double r15222 = -2.0; double r15223 = u1; double r15224 = log(r15223); double r15225 = r15222 * r15224; double r15226 = 0.5; double r15227 = pow(r15225, r15226); double r15228 = r15221 * r15227; double r15229 = 2.0; double r15230 = atan2(1.0, 0.0); double r15231 = r15229 * r15230; double r15232 = u2; double r15233 = r15231 * r15232; double r15234 = cos(r15233); double r15235 = r15228 * r15234; double r15236 = r15235 + r15226; return r15236; } double f_of(float u1, float u2) { float r15237 = 0.5f; float r15238 = -2.0f; float r15239 = u1; float r15240 = log(r15239); float r15241 = r15238 * r15240; float r15242 = pow(r15241, r15237); float r15243 = 6.0f; float r15244 = atan2(1.0, 0.0); float r15245 = u2; float r15246 = 2.0f; float r15247 = r15245 * r15246; float r15248 = r15244 * r15247; float r15249 = cos(r15248); float r15250 = r15249 * (r15249 * r15249); float r15251 = cbrt(r15250); float r15252 = r15243 / r15251; float r15253 = r15242 / r15252; float r15254 = r15237 + r15253; return r15254; } double f_od(double u1, double u2) { double r15255 = 0.5; double r15256 = -2.0; double r15257 = u1; double r15258 = log(r15257); double r15259 = r15256 * r15258; double r15260 = pow(r15259, r15255); double r15261 = 6.0; double r15262 = atan2(1.0, 0.0); double r15263 = u2; double r15264 = 2.0; double r15265 = r15263 * r15264; double r15266 = r15262 * r15265; double r15267 = cos(r15266); double r15268 = r15267 * (r15267 * r15267); double r15269 = cbrt(r15268); double r15270 = r15261 / r15269; double r15271 = r15260 / r15270; double r15272 = r15255 + r15271; return r15272; } 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 r15273, r15274, r15275, r15276, r15277, r15278, r15279, r15280, r15281, r15282, r15283, r15284, r15285, r15286, r15287, r15288, r15289, r15290; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r15273, "1", 10, MPFR_RNDN); mpfr_init_set_str(r15274, "6", 10, MPFR_RNDN); mpfr_init(r15275); mpfr_init_set_str(r15276, "-2", 10, MPFR_RNDN); mpfr_init(r15277); mpfr_init(r15278); mpfr_init(r15279); mpfr_init_set_str(r15280, "0.5", 10, MPFR_RNDN); mpfr_init(r15281); mpfr_init(r15282); mpfr_init_set_str(r15283, "2", 10, MPFR_RNDN); mpfr_init(r15284); mpfr_init(r15285); mpfr_init(r15286); mpfr_init(r15287); mpfr_init(r15288); mpfr_init(r15289); mpfr_init(r15290); } double f_im(double u1, double u2) { ; ; mpfr_div(r15275, r15273, r15274, MPFR_RNDN); ; mpfr_set_d(r15277, u1, MPFR_RNDN); mpfr_log(r15278, r15277, MPFR_RNDN); mpfr_mul(r15279, r15276, r15278, MPFR_RNDN); ; mpfr_pow(r15281, r15279, r15280, MPFR_RNDN); mpfr_mul(r15282, r15275, r15281, MPFR_RNDN); ; mpfr_const_pi(r15284, MPFR_RNDN); mpfr_mul(r15285, r15283, r15284, MPFR_RNDN); mpfr_set_d(r15286, u2, MPFR_RNDN); mpfr_mul(r15287, r15285, r15286, MPFR_RNDN); mpfr_cos(r15288, r15287, MPFR_RNDN); mpfr_mul(r15289, r15282, r15288, MPFR_RNDN); mpfr_add(r15290, r15289, r15280, MPFR_RNDN); return mpfr_get_d(r15290, MPFR_RNDN); } 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_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r15291, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r15292, "-2", 10, MPFR_RNDN); mpfr_init(r15293); mpfr_init(r15294); mpfr_init(r15295); mpfr_init(r15296); mpfr_init_set_str(r15297, "6", 10, MPFR_RNDN); mpfr_init(r15298); mpfr_init(r15299); mpfr_init_set_str(r15300, "2", 10, MPFR_RNDN); mpfr_init(r15301); mpfr_init(r15302); mpfr_init(r15303); mpfr_init(r15304); mpfr_init(r15305); mpfr_init(r15306); mpfr_init(r15307); mpfr_init(r15308); } double f_fm(double u1, double u2) { ; ; mpfr_set_d(r15293, u1, MPFR_RNDN); mpfr_log(r15294, r15293, MPFR_RNDN); mpfr_mul(r15295, r15292, r15294, MPFR_RNDN); mpfr_pow(r15296, r15295, r15291, MPFR_RNDN); ; mpfr_const_pi(r15298, MPFR_RNDN); mpfr_set_d(r15299, u2, MPFR_RNDN); ; mpfr_mul(r15301, r15299, r15300, MPFR_RNDN); mpfr_mul(r15302, r15298, r15301, MPFR_RNDN); mpfr_cos(r15303, r15302, MPFR_RNDN); mpfr_mul(r15304, r15303, r15303, MPFR_RNDN); mpfr_mul(r15304, r15304, r15303, MPFR_RNDN); mpfr_cbrt(r15305, r15304, MPFR_RNDN); mpfr_div(r15306, r15297, r15305, MPFR_RNDN); mpfr_div(r15307, r15296, r15306, MPFR_RNDN); mpfr_add(r15308, r15291, r15307, 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; void setup_mpfr_f_dm() { 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_set_str(r15318, "2", 10, MPFR_RNDN); mpfr_init(r15319); mpfr_init(r15320); mpfr_init(r15321); mpfr_init(r15322); mpfr_init(r15323); mpfr_init(r15324); mpfr_init(r15325); mpfr_init(r15326); } double f_dm(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_set_d(r15317, u2, MPFR_RNDN); ; mpfr_mul(r15319, r15317, r15318, MPFR_RNDN); mpfr_mul(r15320, r15316, r15319, MPFR_RNDN); mpfr_cos(r15321, r15320, MPFR_RNDN); mpfr_mul(r15322, r15321, r15321, MPFR_RNDN); mpfr_mul(r15322, r15322, r15321, MPFR_RNDN); mpfr_cbrt(r15323, r15322, MPFR_RNDN); mpfr_div(r15324, r15315, r15323, MPFR_RNDN); mpfr_div(r15325, r15314, r15324, MPFR_RNDN); mpfr_add(r15326, r15309, r15325, MPFR_RNDN); return mpfr_get_d(r15326, MPFR_RNDN); }