#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r20665 = 1; float r20666 = 6; float r20667 = r20665 / r20666; float r20668 = -2; float r20669 = u1; float r20670 = log(r20669); float r20671 = r20668 * r20670; float r20672 = 0.5; float r20673 = pow(r20671, r20672); float r20674 = r20667 * r20673; float r20675 = 2; float r20676 = atan2(1.0, 0.0); float r20677 = r20675 * r20676; float r20678 = u2; float r20679 = r20677 * r20678; float r20680 = cos(r20679); float r20681 = r20674 * r20680; float r20682 = r20681 + r20672; return r20682; } double f_id(double u1, double u2) { double r20683 = 1; double r20684 = 6; double r20685 = r20683 / r20684; double r20686 = -2; double r20687 = u1; double r20688 = log(r20687); double r20689 = r20686 * r20688; double r20690 = 0.5; double r20691 = pow(r20689, r20690); double r20692 = r20685 * r20691; double r20693 = 2; double r20694 = atan2(1.0, 0.0); double r20695 = r20693 * r20694; double r20696 = u2; double r20697 = r20695 * r20696; double r20698 = cos(r20697); double r20699 = r20692 * r20698; double r20700 = r20699 + r20690; return r20700; } double f_of(float u1, float u2) { float r20701 = 0.5; float r20702 = u2; float r20703 = r20702 + r20702; float r20704 = atan2(1.0, 0.0); float r20705 = r20703 * r20704; float r20706 = cos(r20705); float r20707 = 6; float r20708 = r20706 / r20707; float r20709 = 1; float r20710 = -2; float r20711 = 1.0; float r20712 = pow(r20710, r20711); float r20713 = r20709 / r20712; float r20714 = u1; float r20715 = log(r20714); float r20716 = pow(r20715, r20711); float r20717 = r20713 / r20716; float r20718 = pow(r20717, r20701); float r20719 = r20708 / r20718; float r20720 = r20701 + r20719; return r20720; } double f_od(double u1, double u2) { double r20721 = 0.5; double r20722 = u2; double r20723 = r20722 + r20722; double r20724 = atan2(1.0, 0.0); double r20725 = r20723 * r20724; double r20726 = cos(r20725); double r20727 = 6; double r20728 = r20726 / r20727; double r20729 = 1; double r20730 = -2; double r20731 = 1.0; double r20732 = pow(r20730, r20731); double r20733 = r20729 / r20732; double r20734 = u1; double r20735 = log(r20734); double r20736 = pow(r20735, r20731); double r20737 = r20733 / r20736; double r20738 = pow(r20737, r20721); double r20739 = r20728 / r20738; double r20740 = r20721 + r20739; return r20740; } 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 r20741, r20742, r20743, r20744, r20745, r20746, r20747, r20748, r20749, r20750, r20751, r20752, r20753, r20754, r20755, r20756, r20757, r20758; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20741, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20742, "6", 10, MPFR_RNDN); mpfr_init(r20743); mpfr_init_set_str(r20744, "-2", 10, MPFR_RNDN); mpfr_init(r20745); mpfr_init(r20746); mpfr_init(r20747); mpfr_init_set_str(r20748, "0.5", 10, MPFR_RNDN); mpfr_init(r20749); mpfr_init(r20750); mpfr_init_set_str(r20751, "2", 10, MPFR_RNDN); mpfr_init(r20752); mpfr_init(r20753); mpfr_init(r20754); mpfr_init(r20755); mpfr_init(r20756); mpfr_init(r20757); mpfr_init(r20758); } double f_im(double u1, double u2) { ; ; mpfr_div(r20743, r20741, r20742, MPFR_RNDN); ; mpfr_set_d(r20745, u1, MPFR_RNDN); mpfr_log(r20746, r20745, MPFR_RNDN); mpfr_mul(r20747, r20744, r20746, MPFR_RNDN); ; mpfr_pow(r20749, r20747, r20748, MPFR_RNDN); mpfr_mul(r20750, r20743, r20749, MPFR_RNDN); ; mpfr_const_pi(r20752, MPFR_RNDN); mpfr_mul(r20753, r20751, r20752, MPFR_RNDN); mpfr_set_d(r20754, u2, MPFR_RNDN); mpfr_mul(r20755, r20753, r20754, MPFR_RNDN); mpfr_cos(r20756, r20755, MPFR_RNDN); mpfr_mul(r20757, r20750, r20756, MPFR_RNDN); mpfr_add(r20758, r20757, r20748, MPFR_RNDN); return mpfr_get_d(r20758, MPFR_RNDN); } static mpfr_t r20759, r20760, r20761, r20762, r20763, r20764, r20765, r20766, r20767, r20768, r20769, r20770, r20771, r20772, r20773, r20774, r20775, r20776, r20777, r20778; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20759, "0.5", 10, MPFR_RNDN); mpfr_init(r20760); mpfr_init(r20761); mpfr_init(r20762); mpfr_init(r20763); mpfr_init(r20764); mpfr_init_set_str(r20765, "6", 10, MPFR_RNDN); mpfr_init(r20766); mpfr_init_set_str(r20767, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20768, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r20769, "1.0", 10, MPFR_RNDN); mpfr_init(r20770); mpfr_init(r20771); mpfr_init(r20772); mpfr_init(r20773); mpfr_init(r20774); mpfr_init(r20775); mpfr_init(r20776); mpfr_init(r20777); mpfr_init(r20778); } double f_fm(double u1, double u2) { ; mpfr_set_d(r20760, u2, MPFR_RNDN); mpfr_add(r20761, r20760, r20760, MPFR_RNDN); mpfr_const_pi(r20762, MPFR_RNDN); mpfr_mul(r20763, r20761, r20762, MPFR_RNDN); mpfr_cos(r20764, r20763, MPFR_RNDN); ; mpfr_div(r20766, r20764, r20765, MPFR_RNDN); ; ; ; mpfr_pow(r20770, r20768, r20769, MPFR_RNDN); mpfr_div(r20771, r20767, r20770, MPFR_RNDN); mpfr_set_d(r20772, u1, MPFR_RNDN); mpfr_log(r20773, r20772, MPFR_RNDN); mpfr_pow(r20774, r20773, r20769, MPFR_RNDN); mpfr_div(r20775, r20771, r20774, MPFR_RNDN); mpfr_pow(r20776, r20775, r20759, MPFR_RNDN); mpfr_div(r20777, r20766, r20776, MPFR_RNDN); mpfr_add(r20778, r20759, r20777, MPFR_RNDN); return mpfr_get_d(r20778, MPFR_RNDN); } static mpfr_t r20779, r20780, r20781, r20782, r20783, r20784, r20785, r20786, r20787, r20788, r20789, r20790, r20791, r20792, r20793, r20794, r20795, r20796, r20797, r20798; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20779, "0.5", 10, MPFR_RNDN); mpfr_init(r20780); mpfr_init(r20781); mpfr_init(r20782); mpfr_init(r20783); mpfr_init(r20784); mpfr_init_set_str(r20785, "6", 10, MPFR_RNDN); mpfr_init(r20786); mpfr_init_set_str(r20787, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20788, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r20789, "1.0", 10, MPFR_RNDN); mpfr_init(r20790); mpfr_init(r20791); mpfr_init(r20792); mpfr_init(r20793); mpfr_init(r20794); mpfr_init(r20795); mpfr_init(r20796); mpfr_init(r20797); mpfr_init(r20798); } double f_dm(double u1, double u2) { ; mpfr_set_d(r20780, u2, MPFR_RNDN); mpfr_add(r20781, r20780, r20780, MPFR_RNDN); mpfr_const_pi(r20782, MPFR_RNDN); mpfr_mul(r20783, r20781, r20782, MPFR_RNDN); mpfr_cos(r20784, r20783, MPFR_RNDN); ; mpfr_div(r20786, r20784, r20785, MPFR_RNDN); ; ; ; mpfr_pow(r20790, r20788, r20789, MPFR_RNDN); mpfr_div(r20791, r20787, r20790, MPFR_RNDN); mpfr_set_d(r20792, u1, MPFR_RNDN); mpfr_log(r20793, r20792, MPFR_RNDN); mpfr_pow(r20794, r20793, r20789, MPFR_RNDN); mpfr_div(r20795, r20791, r20794, MPFR_RNDN); mpfr_pow(r20796, r20795, r20779, MPFR_RNDN); mpfr_div(r20797, r20786, r20796, MPFR_RNDN); mpfr_add(r20798, r20779, r20797, MPFR_RNDN); return mpfr_get_d(r20798, MPFR_RNDN); }