#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 = 1; float r21726 = 6; float r21727 = r21725 / r21726; float r21728 = sqrt(r21727); float r21729 = -2; float r21730 = u1; float r21731 = log(r21730); float r21732 = r21729 * r21731; float r21733 = 0.5; float r21734 = pow(r21732, r21733); float r21735 = r21728 * r21734; float r21736 = r21728 * r21735; float r21737 = 2; float r21738 = atan2(1.0, 0.0); float r21739 = r21737 * r21738; float r21740 = u2; float r21741 = r21739 * r21740; float r21742 = cos(r21741); float r21743 = r21736 * r21742; float r21744 = r21743 + r21733; return r21744; } double f_od(double u1, double u2) { double r21745 = 1; double r21746 = 6; double r21747 = r21745 / r21746; double r21748 = sqrt(r21747); double r21749 = -2; double r21750 = u1; double r21751 = log(r21750); double r21752 = r21749 * r21751; double r21753 = 0.5; double r21754 = pow(r21752, r21753); double r21755 = r21748 * r21754; double r21756 = r21748 * r21755; double r21757 = 2; double r21758 = atan2(1.0, 0.0); double r21759 = r21757 * r21758; double r21760 = u2; double r21761 = r21759 * r21760; double r21762 = cos(r21761); double r21763 = r21756 * r21762; double r21764 = r21763 + r21753; return r21764; } 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 r21765, r21766, r21767, r21768, r21769, r21770, r21771, r21772, r21773, r21774, r21775, r21776, r21777, r21778, r21779, r21780, r21781, r21782; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r21765, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21766, "6", 10, MPFR_RNDN); mpfr_init(r21767); mpfr_init_set_str(r21768, "-2", 10, MPFR_RNDN); mpfr_init(r21769); mpfr_init(r21770); mpfr_init(r21771); mpfr_init_set_str(r21772, "0.5", 10, MPFR_RNDN); mpfr_init(r21773); mpfr_init(r21774); mpfr_init_set_str(r21775, "2", 10, MPFR_RNDN); mpfr_init(r21776); mpfr_init(r21777); mpfr_init(r21778); mpfr_init(r21779); mpfr_init(r21780); mpfr_init(r21781); mpfr_init(r21782); } double f_im(double u1, double u2) { ; ; mpfr_div(r21767, r21765, r21766, MPFR_RNDN); ; mpfr_set_d(r21769, u1, MPFR_RNDN); mpfr_log(r21770, r21769, MPFR_RNDN); mpfr_mul(r21771, r21768, r21770, MPFR_RNDN); ; mpfr_pow(r21773, r21771, r21772, MPFR_RNDN); mpfr_mul(r21774, r21767, r21773, MPFR_RNDN); ; mpfr_const_pi(r21776, MPFR_RNDN); mpfr_mul(r21777, r21775, r21776, MPFR_RNDN); mpfr_set_d(r21778, u2, MPFR_RNDN); mpfr_mul(r21779, r21777, r21778, MPFR_RNDN); mpfr_cos(r21780, r21779, MPFR_RNDN); mpfr_mul(r21781, r21774, r21780, MPFR_RNDN); mpfr_add(r21782, r21781, r21772, MPFR_RNDN); return mpfr_get_d(r21782, MPFR_RNDN); } static mpfr_t r21783, r21784, r21785, r21786, r21787, r21788, r21789, r21790, r21791, r21792, r21793, r21794, r21795, r21796, r21797, r21798, r21799, r21800, r21801, r21802; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r21783, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21784, "6", 10, MPFR_RNDN); mpfr_init(r21785); mpfr_init(r21786); mpfr_init_set_str(r21787, "-2", 10, MPFR_RNDN); mpfr_init(r21788); mpfr_init(r21789); mpfr_init(r21790); mpfr_init_set_str(r21791, "0.5", 10, MPFR_RNDN); mpfr_init(r21792); mpfr_init(r21793); mpfr_init(r21794); mpfr_init_set_str(r21795, "2", 10, MPFR_RNDN); mpfr_init(r21796); mpfr_init(r21797); mpfr_init(r21798); mpfr_init(r21799); mpfr_init(r21800); mpfr_init(r21801); mpfr_init(r21802); } double f_fm(double u1, double u2) { ; ; mpfr_div(r21785, r21783, r21784, MPFR_RNDN); mpfr_sqrt(r21786, r21785, MPFR_RNDN); ; mpfr_set_d(r21788, u1, MPFR_RNDN); mpfr_log(r21789, r21788, MPFR_RNDN); mpfr_mul(r21790, r21787, r21789, MPFR_RNDN); ; mpfr_pow(r21792, r21790, r21791, MPFR_RNDN); mpfr_mul(r21793, r21786, r21792, MPFR_RNDN); mpfr_mul(r21794, r21786, r21793, MPFR_RNDN); ; mpfr_const_pi(r21796, MPFR_RNDN); mpfr_mul(r21797, r21795, r21796, MPFR_RNDN); mpfr_set_d(r21798, u2, MPFR_RNDN); mpfr_mul(r21799, r21797, r21798, MPFR_RNDN); mpfr_cos(r21800, r21799, MPFR_RNDN); mpfr_mul(r21801, r21794, r21800, MPFR_RNDN); mpfr_add(r21802, r21801, r21791, MPFR_RNDN); return mpfr_get_d(r21802, MPFR_RNDN); } static mpfr_t r21803, r21804, r21805, r21806, r21807, r21808, r21809, r21810, r21811, r21812, r21813, r21814, r21815, r21816, r21817, r21818, r21819, r21820, r21821, r21822; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r21803, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21804, "6", 10, MPFR_RNDN); mpfr_init(r21805); mpfr_init(r21806); mpfr_init_set_str(r21807, "-2", 10, MPFR_RNDN); mpfr_init(r21808); mpfr_init(r21809); mpfr_init(r21810); mpfr_init_set_str(r21811, "0.5", 10, MPFR_RNDN); mpfr_init(r21812); mpfr_init(r21813); mpfr_init(r21814); mpfr_init_set_str(r21815, "2", 10, MPFR_RNDN); mpfr_init(r21816); mpfr_init(r21817); mpfr_init(r21818); mpfr_init(r21819); mpfr_init(r21820); mpfr_init(r21821); mpfr_init(r21822); } double f_dm(double u1, double u2) { ; ; mpfr_div(r21805, r21803, r21804, MPFR_RNDN); mpfr_sqrt(r21806, r21805, MPFR_RNDN); ; mpfr_set_d(r21808, u1, MPFR_RNDN); mpfr_log(r21809, r21808, MPFR_RNDN); mpfr_mul(r21810, r21807, r21809, MPFR_RNDN); ; mpfr_pow(r21812, r21810, r21811, MPFR_RNDN); mpfr_mul(r21813, r21806, r21812, MPFR_RNDN); mpfr_mul(r21814, r21806, r21813, MPFR_RNDN); ; mpfr_const_pi(r21816, MPFR_RNDN); mpfr_mul(r21817, r21815, r21816, MPFR_RNDN); mpfr_set_d(r21818, u2, MPFR_RNDN); mpfr_mul(r21819, r21817, r21818, MPFR_RNDN); mpfr_cos(r21820, r21819, MPFR_RNDN); mpfr_mul(r21821, r21814, r21820, MPFR_RNDN); mpfr_add(r21822, r21821, r21811, MPFR_RNDN); return mpfr_get_d(r21822, MPFR_RNDN); }