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