#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r8268 = 1; float r8269 = 6; float r8270 = r8268 / r8269; float r8271 = -2; float r8272 = u1; float r8273 = log(r8272); float r8274 = r8271 * r8273; float r8275 = 0.5; float r8276 = pow(r8274, r8275); float r8277 = r8270 * r8276; float r8278 = 2; float r8279 = atan2(1.0, 0.0); float r8280 = r8278 * r8279; float r8281 = u2; float r8282 = r8280 * r8281; float r8283 = cos(r8282); float r8284 = r8277 * r8283; float r8285 = r8284 + r8275; return r8285; } double f_id(double u1, double u2) { double r8286 = 1; double r8287 = 6; double r8288 = r8286 / r8287; double r8289 = -2; double r8290 = u1; double r8291 = log(r8290); double r8292 = r8289 * r8291; double r8293 = 0.5; double r8294 = pow(r8292, r8293); double r8295 = r8288 * r8294; double r8296 = 2; double r8297 = atan2(1.0, 0.0); double r8298 = r8296 * r8297; double r8299 = u2; double r8300 = r8298 * r8299; double r8301 = cos(r8300); double r8302 = r8295 * r8301; double r8303 = r8302 + r8293; return r8303; } double f_of(float u1, float u2) { float r8304 = 1; float r8305 = 6; float r8306 = r8304 / r8305; float r8307 = -2; float r8308 = u1; float r8309 = log(r8308); float r8310 = r8307 * r8309; float r8311 = 0.5; float r8312 = pow(r8310, r8311); float r8313 = r8306 * r8312; float r8314 = 2; float r8315 = atan2(1.0, 0.0); float r8316 = r8314 * r8315; float r8317 = u2; float r8318 = r8316 * r8317; float r8319 = cos(r8318); float r8320 = r8313 * r8319; float r8321 = r8320 + r8311; return r8321; } double f_od(double u1, double u2) { double r8322 = 1; double r8323 = 6; double r8324 = r8322 / r8323; double r8325 = -2; double r8326 = u1; double r8327 = log(r8326); double r8328 = r8325 * r8327; double r8329 = 0.5; double r8330 = pow(r8328, r8329); double r8331 = r8324 * r8330; double r8332 = 2; double r8333 = atan2(1.0, 0.0); double r8334 = r8332 * r8333; double r8335 = u2; double r8336 = r8334 * r8335; double r8337 = cos(r8336); double r8338 = r8331 * r8337; double r8339 = r8338 + r8329; return r8339; } 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 r8340, r8341, r8342, r8343, r8344, r8345, r8346, r8347, r8348, r8349, r8350, r8351, r8352, r8353, r8354, r8355, r8356, r8357; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r8340, "1", 10, MPFR_RNDN); mpfr_init_set_str(r8341, "6", 10, MPFR_RNDN); mpfr_init(r8342); mpfr_init_set_str(r8343, "-2", 10, MPFR_RNDN); mpfr_init(r8344); mpfr_init(r8345); mpfr_init(r8346); mpfr_init_set_str(r8347, "0.5", 10, MPFR_RNDN); mpfr_init(r8348); mpfr_init(r8349); mpfr_init_set_str(r8350, "2", 10, MPFR_RNDN); mpfr_init(r8351); mpfr_init(r8352); mpfr_init(r8353); mpfr_init(r8354); mpfr_init(r8355); mpfr_init(r8356); mpfr_init(r8357); } double f_im(double u1, double u2) { ; ; mpfr_div(r8342, r8340, r8341, MPFR_RNDN); ; mpfr_set_d(r8344, u1, MPFR_RNDN); mpfr_log(r8345, r8344, MPFR_RNDN); mpfr_mul(r8346, r8343, r8345, MPFR_RNDN); ; mpfr_pow(r8348, r8346, r8347, MPFR_RNDN); mpfr_mul(r8349, r8342, r8348, MPFR_RNDN); ; mpfr_const_pi(r8351, MPFR_RNDN); mpfr_mul(r8352, r8350, r8351, MPFR_RNDN); mpfr_set_d(r8353, u2, MPFR_RNDN); mpfr_mul(r8354, r8352, r8353, MPFR_RNDN); mpfr_cos(r8355, r8354, MPFR_RNDN); mpfr_mul(r8356, r8349, r8355, MPFR_RNDN); mpfr_add(r8357, r8356, r8347, MPFR_RNDN); return mpfr_get_d(r8357, MPFR_RNDN); } static mpfr_t r8358, r8359, r8360, r8361, r8362, r8363, r8364, r8365, r8366, r8367, r8368, r8369, r8370, r8371, r8372, r8373, r8374, r8375; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r8358, "1", 10, MPFR_RNDN); mpfr_init_set_str(r8359, "6", 10, MPFR_RNDN); mpfr_init(r8360); mpfr_init_set_str(r8361, "-2", 10, MPFR_RNDN); mpfr_init(r8362); mpfr_init(r8363); mpfr_init(r8364); mpfr_init_set_str(r8365, "0.5", 10, MPFR_RNDN); mpfr_init(r8366); mpfr_init(r8367); mpfr_init_set_str(r8368, "2", 10, MPFR_RNDN); mpfr_init(r8369); mpfr_init(r8370); mpfr_init(r8371); mpfr_init(r8372); mpfr_init(r8373); mpfr_init(r8374); mpfr_init(r8375); } double f_fm(double u1, double u2) { ; ; mpfr_div(r8360, r8358, r8359, MPFR_RNDN); ; mpfr_set_d(r8362, u1, MPFR_RNDN); mpfr_log(r8363, r8362, MPFR_RNDN); mpfr_mul(r8364, r8361, r8363, MPFR_RNDN); ; mpfr_pow(r8366, r8364, r8365, MPFR_RNDN); mpfr_mul(r8367, r8360, r8366, MPFR_RNDN); ; mpfr_const_pi(r8369, MPFR_RNDN); mpfr_mul(r8370, r8368, r8369, MPFR_RNDN); mpfr_set_d(r8371, u2, MPFR_RNDN); mpfr_mul(r8372, r8370, r8371, MPFR_RNDN); mpfr_cos(r8373, r8372, MPFR_RNDN); mpfr_mul(r8374, r8367, r8373, MPFR_RNDN); mpfr_add(r8375, r8374, r8365, MPFR_RNDN); return mpfr_get_d(r8375, MPFR_RNDN); } static mpfr_t r8376, r8377, r8378, r8379, r8380, r8381, r8382, r8383, r8384, r8385, r8386, r8387, r8388, r8389, r8390, r8391, r8392, r8393; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r8376, "1", 10, MPFR_RNDN); mpfr_init_set_str(r8377, "6", 10, MPFR_RNDN); mpfr_init(r8378); mpfr_init_set_str(r8379, "-2", 10, MPFR_RNDN); mpfr_init(r8380); mpfr_init(r8381); mpfr_init(r8382); mpfr_init_set_str(r8383, "0.5", 10, MPFR_RNDN); mpfr_init(r8384); mpfr_init(r8385); mpfr_init_set_str(r8386, "2", 10, MPFR_RNDN); mpfr_init(r8387); mpfr_init(r8388); mpfr_init(r8389); mpfr_init(r8390); mpfr_init(r8391); mpfr_init(r8392); mpfr_init(r8393); } double f_dm(double u1, double u2) { ; ; mpfr_div(r8378, r8376, r8377, MPFR_RNDN); ; mpfr_set_d(r8380, u1, MPFR_RNDN); mpfr_log(r8381, r8380, MPFR_RNDN); mpfr_mul(r8382, r8379, r8381, MPFR_RNDN); ; mpfr_pow(r8384, r8382, r8383, MPFR_RNDN); mpfr_mul(r8385, r8378, r8384, MPFR_RNDN); ; mpfr_const_pi(r8387, MPFR_RNDN); mpfr_mul(r8388, r8386, r8387, MPFR_RNDN); mpfr_set_d(r8389, u2, MPFR_RNDN); mpfr_mul(r8390, r8388, r8389, MPFR_RNDN); mpfr_cos(r8391, r8390, MPFR_RNDN); mpfr_mul(r8392, r8385, r8391, MPFR_RNDN); mpfr_add(r8393, r8392, r8383, MPFR_RNDN); return mpfr_get_d(r8393, MPFR_RNDN); }