#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r21606 = 1; float r21607 = 6; float r21608 = r21606 / r21607; float r21609 = -2; float r21610 = u1; float r21611 = log(r21610); float r21612 = r21609 * r21611; float r21613 = 0.5; float r21614 = pow(r21612, r21613); float r21615 = r21608 * r21614; float r21616 = 2; float r21617 = atan2(1.0, 0.0); float r21618 = r21616 * r21617; float r21619 = u2; float r21620 = r21618 * r21619; float r21621 = cos(r21620); float r21622 = r21615 * r21621; float r21623 = r21622 + r21613; return r21623; } double f_id(double u1, double u2) { double r21624 = 1; double r21625 = 6; double r21626 = r21624 / r21625; double r21627 = -2; double r21628 = u1; double r21629 = log(r21628); double r21630 = r21627 * r21629; double r21631 = 0.5; double r21632 = pow(r21630, r21631); double r21633 = r21626 * r21632; double r21634 = 2; double r21635 = atan2(1.0, 0.0); double r21636 = r21634 * r21635; double r21637 = u2; double r21638 = r21636 * r21637; double r21639 = cos(r21638); double r21640 = r21633 * r21639; double r21641 = r21640 + r21631; return r21641; } double f_of(float u1, float u2) { float r21642 = 1; float r21643 = 6; float r21644 = r21642 / r21643; float r21645 = sqrt(r21644); float r21646 = -2; float r21647 = u1; float r21648 = log(r21647); float r21649 = r21646 * r21648; float r21650 = 0.5; float r21651 = pow(r21649, r21650); float r21652 = r21645 * r21651; float r21653 = r21645 * r21652; float r21654 = 2; float r21655 = atan2(1.0, 0.0); float r21656 = r21654 * r21655; float r21657 = u2; float r21658 = r21656 * r21657; float r21659 = cos(r21658); float r21660 = r21653 * r21659; float r21661 = r21660 + r21650; return r21661; } double f_od(double u1, double u2) { double r21662 = 1; double r21663 = 6; double r21664 = r21662 / r21663; double r21665 = sqrt(r21664); double r21666 = -2; double r21667 = u1; double r21668 = log(r21667); double r21669 = r21666 * r21668; double r21670 = 0.5; double r21671 = pow(r21669, r21670); double r21672 = r21665 * r21671; double r21673 = r21665 * r21672; double r21674 = 2; double r21675 = atan2(1.0, 0.0); double r21676 = r21674 * r21675; double r21677 = u2; double r21678 = r21676 * r21677; double r21679 = cos(r21678); double r21680 = r21673 * r21679; double r21681 = r21680 + r21670; return r21681; } 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 r21682, r21683, r21684, r21685, r21686, r21687, r21688, r21689, r21690, r21691, r21692, r21693, r21694, r21695, r21696, r21697, r21698, r21699; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21682, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21683, "6", 10, MPFR_RNDN); mpfr_init(r21684); mpfr_init_set_str(r21685, "-2", 10, MPFR_RNDN); mpfr_init(r21686); mpfr_init(r21687); mpfr_init(r21688); mpfr_init_set_str(r21689, "0.5", 10, MPFR_RNDN); mpfr_init(r21690); mpfr_init(r21691); mpfr_init_set_str(r21692, "2", 10, MPFR_RNDN); mpfr_init(r21693); mpfr_init(r21694); mpfr_init(r21695); mpfr_init(r21696); mpfr_init(r21697); mpfr_init(r21698); mpfr_init(r21699); } double f_im(double u1, double u2) { ; ; mpfr_div(r21684, r21682, r21683, MPFR_RNDN); ; mpfr_set_d(r21686, u1, MPFR_RNDN); mpfr_log(r21687, r21686, MPFR_RNDN); mpfr_mul(r21688, r21685, r21687, MPFR_RNDN); ; mpfr_pow(r21690, r21688, r21689, MPFR_RNDN); mpfr_mul(r21691, r21684, r21690, MPFR_RNDN); ; mpfr_const_pi(r21693, MPFR_RNDN); mpfr_mul(r21694, r21692, r21693, MPFR_RNDN); mpfr_set_d(r21695, u2, MPFR_RNDN); mpfr_mul(r21696, r21694, r21695, MPFR_RNDN); mpfr_cos(r21697, r21696, MPFR_RNDN); mpfr_mul(r21698, r21691, r21697, MPFR_RNDN); mpfr_add(r21699, r21698, r21689, MPFR_RNDN); return mpfr_get_d(r21699, MPFR_RNDN); } static mpfr_t r21700, r21701, r21702, r21703, r21704, r21705, r21706, r21707, r21708, r21709, r21710, r21711, r21712, r21713, r21714, r21715, r21716, r21717, r21718, r21719; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21700, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21701, "6", 10, MPFR_RNDN); mpfr_init(r21702); mpfr_init(r21703); mpfr_init_set_str(r21704, "-2", 10, MPFR_RNDN); mpfr_init(r21705); mpfr_init(r21706); mpfr_init(r21707); mpfr_init_set_str(r21708, "0.5", 10, MPFR_RNDN); mpfr_init(r21709); mpfr_init(r21710); mpfr_init(r21711); mpfr_init_set_str(r21712, "2", 10, MPFR_RNDN); mpfr_init(r21713); mpfr_init(r21714); mpfr_init(r21715); mpfr_init(r21716); mpfr_init(r21717); mpfr_init(r21718); mpfr_init(r21719); } double f_fm(double u1, double u2) { ; ; mpfr_div(r21702, r21700, r21701, MPFR_RNDN); mpfr_sqrt(r21703, r21702, MPFR_RNDN); ; mpfr_set_d(r21705, u1, MPFR_RNDN); mpfr_log(r21706, r21705, MPFR_RNDN); mpfr_mul(r21707, r21704, r21706, MPFR_RNDN); ; mpfr_pow(r21709, r21707, r21708, MPFR_RNDN); mpfr_mul(r21710, r21703, r21709, MPFR_RNDN); mpfr_mul(r21711, r21703, r21710, MPFR_RNDN); ; mpfr_const_pi(r21713, MPFR_RNDN); mpfr_mul(r21714, r21712, r21713, MPFR_RNDN); mpfr_set_d(r21715, u2, MPFR_RNDN); mpfr_mul(r21716, r21714, r21715, MPFR_RNDN); mpfr_cos(r21717, r21716, MPFR_RNDN); mpfr_mul(r21718, r21711, r21717, MPFR_RNDN); mpfr_add(r21719, r21718, r21708, MPFR_RNDN); return mpfr_get_d(r21719, MPFR_RNDN); } static mpfr_t r21720, r21721, r21722, r21723, r21724, r21725, r21726, r21727, r21728, r21729, r21730, r21731, r21732, r21733, r21734, r21735, r21736, r21737, r21738, r21739; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21720, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21721, "6", 10, MPFR_RNDN); mpfr_init(r21722); mpfr_init(r21723); mpfr_init_set_str(r21724, "-2", 10, MPFR_RNDN); mpfr_init(r21725); mpfr_init(r21726); mpfr_init(r21727); mpfr_init_set_str(r21728, "0.5", 10, MPFR_RNDN); mpfr_init(r21729); mpfr_init(r21730); mpfr_init(r21731); mpfr_init_set_str(r21732, "2", 10, MPFR_RNDN); mpfr_init(r21733); mpfr_init(r21734); mpfr_init(r21735); mpfr_init(r21736); mpfr_init(r21737); mpfr_init(r21738); mpfr_init(r21739); } double f_dm(double u1, double u2) { ; ; mpfr_div(r21722, r21720, r21721, MPFR_RNDN); mpfr_sqrt(r21723, r21722, MPFR_RNDN); ; mpfr_set_d(r21725, u1, MPFR_RNDN); mpfr_log(r21726, r21725, MPFR_RNDN); mpfr_mul(r21727, r21724, r21726, MPFR_RNDN); ; mpfr_pow(r21729, r21727, r21728, MPFR_RNDN); mpfr_mul(r21730, r21723, r21729, MPFR_RNDN); mpfr_mul(r21731, r21723, r21730, MPFR_RNDN); ; mpfr_const_pi(r21733, MPFR_RNDN); mpfr_mul(r21734, r21732, r21733, MPFR_RNDN); mpfr_set_d(r21735, u2, MPFR_RNDN); mpfr_mul(r21736, r21734, r21735, MPFR_RNDN); mpfr_cos(r21737, r21736, MPFR_RNDN); mpfr_mul(r21738, r21731, r21737, MPFR_RNDN); mpfr_add(r21739, r21738, r21728, MPFR_RNDN); return mpfr_get_d(r21739, MPFR_RNDN); }