#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r21689 = 1; float r21690 = 6; float r21691 = r21689 / r21690; float r21692 = -2; float r21693 = u1; float r21694 = log(r21693); float r21695 = r21692 * r21694; float r21696 = 0.5; float r21697 = pow(r21695, r21696); float r21698 = r21691 * r21697; float r21699 = 2; float r21700 = atan2(1.0, 0.0); float r21701 = r21699 * r21700; float r21702 = u2; float r21703 = r21701 * r21702; float r21704 = cos(r21703); float r21705 = r21698 * r21704; float r21706 = r21705 + r21696; return r21706; } double f_id(double u1, double u2) { double r21707 = 1; double r21708 = 6; double r21709 = r21707 / r21708; double r21710 = -2; double r21711 = u1; double r21712 = log(r21711); double r21713 = r21710 * r21712; double r21714 = 0.5; double r21715 = pow(r21713, r21714); double r21716 = r21709 * r21715; double r21717 = 2; double r21718 = atan2(1.0, 0.0); double r21719 = r21717 * r21718; double r21720 = u2; double r21721 = r21719 * r21720; double r21722 = cos(r21721); double r21723 = r21716 * r21722; double r21724 = r21723 + r21714; return r21724; } double f_of(float u1, float u2) { float r21725 = u2; float r21726 = 2; float r21727 = r21725 * r21726; float r21728 = atan2(1.0, 0.0); float r21729 = r21727 * r21728; float r21730 = cos(r21729); float r21731 = 6; float r21732 = u1; float r21733 = log(r21732); float r21734 = -2; float r21735 = r21733 * r21734; float r21736 = 0.5; float r21737 = pow(r21735, r21736); float r21738 = r21731 / r21737; float r21739 = r21730 / r21738; float r21740 = r21739 + r21736; return r21740; } double f_od(double u1, double u2) { double r21741 = u2; double r21742 = 2; double r21743 = r21741 * r21742; double r21744 = atan2(1.0, 0.0); double r21745 = r21743 * r21744; double r21746 = cos(r21745); double r21747 = 6; double r21748 = u1; double r21749 = log(r21748); double r21750 = -2; double r21751 = r21749 * r21750; double r21752 = 0.5; double r21753 = pow(r21751, r21752); double r21754 = r21747 / r21753; double r21755 = r21746 / r21754; double r21756 = r21755 + r21752; return r21756; } 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 r21757, r21758, r21759, r21760, r21761, r21762, r21763, r21764, r21765, r21766, r21767, r21768, r21769, r21770, r21771, r21772, r21773, r21774; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r21757, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21758, "6", 10, MPFR_RNDN); mpfr_init(r21759); mpfr_init_set_str(r21760, "-2", 10, MPFR_RNDN); mpfr_init(r21761); mpfr_init(r21762); mpfr_init(r21763); mpfr_init_set_str(r21764, "0.5", 10, MPFR_RNDN); mpfr_init(r21765); mpfr_init(r21766); mpfr_init_set_str(r21767, "2", 10, MPFR_RNDN); mpfr_init(r21768); mpfr_init(r21769); mpfr_init(r21770); mpfr_init(r21771); mpfr_init(r21772); mpfr_init(r21773); mpfr_init(r21774); } double f_im(double u1, double u2) { ; ; mpfr_div(r21759, r21757, r21758, MPFR_RNDN); ; mpfr_set_d(r21761, u1, MPFR_RNDN); mpfr_log(r21762, r21761, MPFR_RNDN); mpfr_mul(r21763, r21760, r21762, MPFR_RNDN); ; mpfr_pow(r21765, r21763, r21764, MPFR_RNDN); mpfr_mul(r21766, r21759, r21765, MPFR_RNDN); ; mpfr_const_pi(r21768, MPFR_RNDN); mpfr_mul(r21769, r21767, r21768, MPFR_RNDN); mpfr_set_d(r21770, u2, MPFR_RNDN); mpfr_mul(r21771, r21769, r21770, MPFR_RNDN); mpfr_cos(r21772, r21771, MPFR_RNDN); mpfr_mul(r21773, r21766, r21772, MPFR_RNDN); mpfr_add(r21774, r21773, r21764, MPFR_RNDN); return mpfr_get_d(r21774, MPFR_RNDN); } static mpfr_t r21775, r21776, r21777, r21778, r21779, r21780, r21781, r21782, r21783, r21784, r21785, r21786, r21787, r21788, r21789, r21790; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r21775); mpfr_init_set_str(r21776, "2", 10, MPFR_RNDN); mpfr_init(r21777); mpfr_init(r21778); mpfr_init(r21779); mpfr_init(r21780); mpfr_init_set_str(r21781, "6", 10, MPFR_RNDN); mpfr_init(r21782); mpfr_init(r21783); mpfr_init_set_str(r21784, "-2", 10, MPFR_RNDN); mpfr_init(r21785); mpfr_init_set_str(r21786, "0.5", 10, MPFR_RNDN); mpfr_init(r21787); mpfr_init(r21788); mpfr_init(r21789); mpfr_init(r21790); } double f_fm(double u1, double u2) { mpfr_set_d(r21775, u2, MPFR_RNDN); ; mpfr_mul(r21777, r21775, r21776, MPFR_RNDN); mpfr_const_pi(r21778, MPFR_RNDN); mpfr_mul(r21779, r21777, r21778, MPFR_RNDN); mpfr_cos(r21780, r21779, MPFR_RNDN); ; mpfr_set_d(r21782, u1, MPFR_RNDN); mpfr_log(r21783, r21782, MPFR_RNDN); ; mpfr_mul(r21785, r21783, r21784, MPFR_RNDN); ; mpfr_pow(r21787, r21785, r21786, MPFR_RNDN); mpfr_div(r21788, r21781, r21787, MPFR_RNDN); mpfr_div(r21789, r21780, r21788, MPFR_RNDN); mpfr_add(r21790, r21789, r21786, MPFR_RNDN); return mpfr_get_d(r21790, MPFR_RNDN); } static mpfr_t r21791, r21792, r21793, r21794, r21795, r21796, r21797, r21798, r21799, r21800, r21801, r21802, r21803, r21804, r21805, r21806; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r21791); mpfr_init_set_str(r21792, "2", 10, MPFR_RNDN); mpfr_init(r21793); mpfr_init(r21794); mpfr_init(r21795); mpfr_init(r21796); mpfr_init_set_str(r21797, "6", 10, MPFR_RNDN); mpfr_init(r21798); mpfr_init(r21799); mpfr_init_set_str(r21800, "-2", 10, MPFR_RNDN); mpfr_init(r21801); mpfr_init_set_str(r21802, "0.5", 10, MPFR_RNDN); mpfr_init(r21803); mpfr_init(r21804); mpfr_init(r21805); mpfr_init(r21806); } double f_dm(double u1, double u2) { mpfr_set_d(r21791, u2, MPFR_RNDN); ; mpfr_mul(r21793, r21791, r21792, MPFR_RNDN); mpfr_const_pi(r21794, MPFR_RNDN); mpfr_mul(r21795, r21793, r21794, MPFR_RNDN); mpfr_cos(r21796, r21795, MPFR_RNDN); ; mpfr_set_d(r21798, u1, MPFR_RNDN); mpfr_log(r21799, r21798, MPFR_RNDN); ; mpfr_mul(r21801, r21799, r21800, MPFR_RNDN); ; mpfr_pow(r21803, r21801, r21802, MPFR_RNDN); mpfr_div(r21804, r21797, r21803, MPFR_RNDN); mpfr_div(r21805, r21796, r21804, MPFR_RNDN); mpfr_add(r21806, r21805, r21802, MPFR_RNDN); return mpfr_get_d(r21806, MPFR_RNDN); }