#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r20585 = 1; float r20586 = 6; float r20587 = r20585 / r20586; float r20588 = -2; float r20589 = u1; float r20590 = log(r20589); float r20591 = r20588 * r20590; float r20592 = 0.5; float r20593 = pow(r20591, r20592); float r20594 = r20587 * r20593; float r20595 = 2; float r20596 = atan2(1.0, 0.0); float r20597 = r20595 * r20596; float r20598 = u2; float r20599 = r20597 * r20598; float r20600 = cos(r20599); float r20601 = r20594 * r20600; float r20602 = r20601 + r20592; return r20602; } double f_id(double u1, double u2) { double r20603 = 1; double r20604 = 6; double r20605 = r20603 / r20604; double r20606 = -2; double r20607 = u1; double r20608 = log(r20607); double r20609 = r20606 * r20608; double r20610 = 0.5; double r20611 = pow(r20609, r20610); double r20612 = r20605 * r20611; double r20613 = 2; double r20614 = atan2(1.0, 0.0); double r20615 = r20613 * r20614; double r20616 = u2; double r20617 = r20615 * r20616; double r20618 = cos(r20617); double r20619 = r20612 * r20618; double r20620 = r20619 + r20610; return r20620; } double f_of(float u1, float u2) { float r20621 = 0.5; float r20622 = u2; float r20623 = r20622 + r20622; float r20624 = atan2(1.0, 0.0); float r20625 = r20623 * r20624; float r20626 = cos(r20625); float r20627 = 6; float r20628 = r20626 / r20627; float r20629 = 1; float r20630 = -2; float r20631 = 1.0; float r20632 = pow(r20630, r20631); float r20633 = r20629 / r20632; float r20634 = u1; float r20635 = log(r20634); float r20636 = pow(r20635, r20631); float r20637 = r20633 / r20636; float r20638 = pow(r20637, r20621); float r20639 = r20628 / r20638; float r20640 = r20621 + r20639; return r20640; } double f_od(double u1, double u2) { double r20641 = 0.5; double r20642 = u2; double r20643 = r20642 + r20642; double r20644 = atan2(1.0, 0.0); double r20645 = r20643 * r20644; double r20646 = cos(r20645); double r20647 = 6; double r20648 = r20646 / r20647; double r20649 = 1; double r20650 = -2; double r20651 = 1.0; double r20652 = pow(r20650, r20651); double r20653 = r20649 / r20652; double r20654 = u1; double r20655 = log(r20654); double r20656 = pow(r20655, r20651); double r20657 = r20653 / r20656; double r20658 = pow(r20657, r20641); double r20659 = r20648 / r20658; double r20660 = r20641 + r20659; return r20660; } 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 r20661, r20662, r20663, r20664, r20665, r20666, r20667, r20668, r20669, r20670, r20671, r20672, r20673, r20674, r20675, r20676, r20677, r20678; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20661, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20662, "6", 10, MPFR_RNDN); mpfr_init(r20663); mpfr_init_set_str(r20664, "-2", 10, MPFR_RNDN); mpfr_init(r20665); mpfr_init(r20666); mpfr_init(r20667); mpfr_init_set_str(r20668, "0.5", 10, MPFR_RNDN); mpfr_init(r20669); mpfr_init(r20670); mpfr_init_set_str(r20671, "2", 10, MPFR_RNDN); mpfr_init(r20672); mpfr_init(r20673); mpfr_init(r20674); mpfr_init(r20675); mpfr_init(r20676); mpfr_init(r20677); mpfr_init(r20678); } double f_im(double u1, double u2) { ; ; mpfr_div(r20663, r20661, r20662, MPFR_RNDN); ; mpfr_set_d(r20665, u1, MPFR_RNDN); mpfr_log(r20666, r20665, MPFR_RNDN); mpfr_mul(r20667, r20664, r20666, MPFR_RNDN); ; mpfr_pow(r20669, r20667, r20668, MPFR_RNDN); mpfr_mul(r20670, r20663, r20669, MPFR_RNDN); ; mpfr_const_pi(r20672, MPFR_RNDN); mpfr_mul(r20673, r20671, r20672, MPFR_RNDN); mpfr_set_d(r20674, u2, MPFR_RNDN); mpfr_mul(r20675, r20673, r20674, MPFR_RNDN); mpfr_cos(r20676, r20675, MPFR_RNDN); mpfr_mul(r20677, r20670, r20676, MPFR_RNDN); mpfr_add(r20678, r20677, r20668, MPFR_RNDN); return mpfr_get_d(r20678, MPFR_RNDN); } static mpfr_t r20679, r20680, r20681, r20682, r20683, r20684, r20685, r20686, r20687, r20688, r20689, r20690, r20691, r20692, r20693, r20694, r20695, r20696, r20697, r20698; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20679, "0.5", 10, MPFR_RNDN); mpfr_init(r20680); mpfr_init(r20681); mpfr_init(r20682); mpfr_init(r20683); mpfr_init(r20684); mpfr_init_set_str(r20685, "6", 10, MPFR_RNDN); mpfr_init(r20686); mpfr_init_set_str(r20687, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20688, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r20689, "1.0", 10, MPFR_RNDN); mpfr_init(r20690); mpfr_init(r20691); mpfr_init(r20692); mpfr_init(r20693); mpfr_init(r20694); mpfr_init(r20695); mpfr_init(r20696); mpfr_init(r20697); mpfr_init(r20698); } double f_fm(double u1, double u2) { ; mpfr_set_d(r20680, u2, MPFR_RNDN); mpfr_add(r20681, r20680, r20680, MPFR_RNDN); mpfr_const_pi(r20682, MPFR_RNDN); mpfr_mul(r20683, r20681, r20682, MPFR_RNDN); mpfr_cos(r20684, r20683, MPFR_RNDN); ; mpfr_div(r20686, r20684, r20685, MPFR_RNDN); ; ; ; mpfr_pow(r20690, r20688, r20689, MPFR_RNDN); mpfr_div(r20691, r20687, r20690, MPFR_RNDN); mpfr_set_d(r20692, u1, MPFR_RNDN); mpfr_log(r20693, r20692, MPFR_RNDN); mpfr_pow(r20694, r20693, r20689, MPFR_RNDN); mpfr_div(r20695, r20691, r20694, MPFR_RNDN); mpfr_pow(r20696, r20695, r20679, MPFR_RNDN); mpfr_div(r20697, r20686, r20696, MPFR_RNDN); mpfr_add(r20698, r20679, r20697, MPFR_RNDN); return mpfr_get_d(r20698, MPFR_RNDN); } static mpfr_t r20699, r20700, r20701, r20702, r20703, r20704, r20705, r20706, r20707, r20708, r20709, r20710, r20711, r20712, r20713, r20714, r20715, r20716, r20717, r20718; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20699, "0.5", 10, MPFR_RNDN); mpfr_init(r20700); mpfr_init(r20701); mpfr_init(r20702); mpfr_init(r20703); mpfr_init(r20704); mpfr_init_set_str(r20705, "6", 10, MPFR_RNDN); mpfr_init(r20706); mpfr_init_set_str(r20707, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20708, "-2", 10, MPFR_RNDN); mpfr_init_set_str(r20709, "1.0", 10, MPFR_RNDN); mpfr_init(r20710); mpfr_init(r20711); mpfr_init(r20712); mpfr_init(r20713); mpfr_init(r20714); mpfr_init(r20715); mpfr_init(r20716); mpfr_init(r20717); mpfr_init(r20718); } double f_dm(double u1, double u2) { ; mpfr_set_d(r20700, u2, MPFR_RNDN); mpfr_add(r20701, r20700, r20700, MPFR_RNDN); mpfr_const_pi(r20702, MPFR_RNDN); mpfr_mul(r20703, r20701, r20702, MPFR_RNDN); mpfr_cos(r20704, r20703, MPFR_RNDN); ; mpfr_div(r20706, r20704, r20705, MPFR_RNDN); ; ; ; mpfr_pow(r20710, r20708, r20709, MPFR_RNDN); mpfr_div(r20711, r20707, r20710, MPFR_RNDN); mpfr_set_d(r20712, u1, MPFR_RNDN); mpfr_log(r20713, r20712, MPFR_RNDN); mpfr_pow(r20714, r20713, r20709, MPFR_RNDN); mpfr_div(r20715, r20711, r20714, MPFR_RNDN); mpfr_pow(r20716, r20715, r20699, MPFR_RNDN); mpfr_div(r20717, r20706, r20716, MPFR_RNDN); mpfr_add(r20718, r20699, r20717, MPFR_RNDN); return mpfr_get_d(r20718, MPFR_RNDN); }