#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r21793 = 1; float r21794 = 6; float r21795 = r21793 / r21794; float r21796 = -2; float r21797 = u1; float r21798 = log(r21797); float r21799 = r21796 * r21798; float r21800 = 0.5; float r21801 = pow(r21799, r21800); float r21802 = r21795 * r21801; float r21803 = 2; float r21804 = atan2(1.0, 0.0); float r21805 = r21803 * r21804; float r21806 = u2; float r21807 = r21805 * r21806; float r21808 = cos(r21807); float r21809 = r21802 * r21808; float r21810 = r21809 + r21800; return r21810; } double f_id(double u1, double u2) { double r21811 = 1; double r21812 = 6; double r21813 = r21811 / r21812; double r21814 = -2; double r21815 = u1; double r21816 = log(r21815); double r21817 = r21814 * r21816; double r21818 = 0.5; double r21819 = pow(r21817, r21818); double r21820 = r21813 * r21819; double r21821 = 2; double r21822 = atan2(1.0, 0.0); double r21823 = r21821 * r21822; double r21824 = u2; double r21825 = r21823 * r21824; double r21826 = cos(r21825); double r21827 = r21820 * r21826; double r21828 = r21827 + r21818; return r21828; } double f_of(float u1, float u2) { float r21829 = 0.5; float r21830 = atan2(1.0, 0.0); float r21831 = u2; float r21832 = r21830 * r21831; float r21833 = 2; float r21834 = r21832 * r21833; float r21835 = cos(r21834); float r21836 = -2; float r21837 = u1; float r21838 = log(r21837); float r21839 = r21836 * r21838; float r21840 = pow(r21839, r21829); float r21841 = 6; float r21842 = r21840 / r21841; float r21843 = r21835 * r21842; float r21844 = r21829 + r21843; return r21844; } double f_od(double u1, double u2) { double r21845 = 0.5; double r21846 = atan2(1.0, 0.0); double r21847 = u2; double r21848 = r21846 * r21847; double r21849 = 2; double r21850 = r21848 * r21849; double r21851 = cos(r21850); double r21852 = -2; double r21853 = u1; double r21854 = log(r21853); double r21855 = r21852 * r21854; double r21856 = pow(r21855, r21845); double r21857 = 6; double r21858 = r21856 / r21857; double r21859 = r21851 * r21858; double r21860 = r21845 + r21859; return r21860; } 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 r21861, r21862, r21863, r21864, r21865, r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r21861, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21862, "6", 10, MPFR_RNDN); mpfr_init(r21863); mpfr_init_set_str(r21864, "-2", 10, MPFR_RNDN); mpfr_init(r21865); mpfr_init(r21866); mpfr_init(r21867); mpfr_init_set_str(r21868, "0.5", 10, MPFR_RNDN); mpfr_init(r21869); mpfr_init(r21870); mpfr_init_set_str(r21871, "2", 10, MPFR_RNDN); mpfr_init(r21872); mpfr_init(r21873); mpfr_init(r21874); mpfr_init(r21875); mpfr_init(r21876); mpfr_init(r21877); mpfr_init(r21878); } double f_im(double u1, double u2) { ; ; mpfr_div(r21863, r21861, r21862, MPFR_RNDN); ; mpfr_set_d(r21865, u1, MPFR_RNDN); mpfr_log(r21866, r21865, MPFR_RNDN); mpfr_mul(r21867, r21864, r21866, MPFR_RNDN); ; mpfr_pow(r21869, r21867, r21868, MPFR_RNDN); mpfr_mul(r21870, r21863, r21869, MPFR_RNDN); ; mpfr_const_pi(r21872, MPFR_RNDN); mpfr_mul(r21873, r21871, r21872, MPFR_RNDN); mpfr_set_d(r21874, u2, MPFR_RNDN); mpfr_mul(r21875, r21873, r21874, MPFR_RNDN); mpfr_cos(r21876, r21875, MPFR_RNDN); mpfr_mul(r21877, r21870, r21876, MPFR_RNDN); mpfr_add(r21878, r21877, r21868, MPFR_RNDN); return mpfr_get_d(r21878, MPFR_RNDN); } static mpfr_t r21879, r21880, r21881, r21882, r21883, r21884, r21885, r21886, r21887, r21888, r21889, r21890, r21891, r21892, r21893, r21894; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r21879, "0.5", 10, MPFR_RNDN); mpfr_init(r21880); mpfr_init(r21881); mpfr_init(r21882); mpfr_init_set_str(r21883, "2", 10, MPFR_RNDN); mpfr_init(r21884); mpfr_init(r21885); mpfr_init_set_str(r21886, "-2", 10, MPFR_RNDN); mpfr_init(r21887); mpfr_init(r21888); mpfr_init(r21889); mpfr_init(r21890); mpfr_init_set_str(r21891, "6", 10, MPFR_RNDN); mpfr_init(r21892); mpfr_init(r21893); mpfr_init(r21894); } double f_fm(double u1, double u2) { ; mpfr_const_pi(r21880, MPFR_RNDN); mpfr_set_d(r21881, u2, MPFR_RNDN); mpfr_mul(r21882, r21880, r21881, MPFR_RNDN); ; mpfr_mul(r21884, r21882, r21883, MPFR_RNDN); mpfr_cos(r21885, r21884, MPFR_RNDN); ; mpfr_set_d(r21887, u1, MPFR_RNDN); mpfr_log(r21888, r21887, MPFR_RNDN); mpfr_mul(r21889, r21886, r21888, MPFR_RNDN); mpfr_pow(r21890, r21889, r21879, MPFR_RNDN); ; mpfr_div(r21892, r21890, r21891, MPFR_RNDN); mpfr_mul(r21893, r21885, r21892, MPFR_RNDN); mpfr_add(r21894, r21879, r21893, MPFR_RNDN); return mpfr_get_d(r21894, MPFR_RNDN); } static mpfr_t r21895, r21896, r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907, r21908, r21909, r21910; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r21895, "0.5", 10, MPFR_RNDN); mpfr_init(r21896); mpfr_init(r21897); mpfr_init(r21898); mpfr_init_set_str(r21899, "2", 10, MPFR_RNDN); mpfr_init(r21900); mpfr_init(r21901); mpfr_init_set_str(r21902, "-2", 10, MPFR_RNDN); mpfr_init(r21903); mpfr_init(r21904); mpfr_init(r21905); mpfr_init(r21906); mpfr_init_set_str(r21907, "6", 10, MPFR_RNDN); mpfr_init(r21908); mpfr_init(r21909); mpfr_init(r21910); } double f_dm(double u1, double u2) { ; mpfr_const_pi(r21896, MPFR_RNDN); mpfr_set_d(r21897, u2, MPFR_RNDN); mpfr_mul(r21898, r21896, r21897, MPFR_RNDN); ; mpfr_mul(r21900, r21898, r21899, MPFR_RNDN); mpfr_cos(r21901, r21900, MPFR_RNDN); ; mpfr_set_d(r21903, u1, MPFR_RNDN); mpfr_log(r21904, r21903, MPFR_RNDN); mpfr_mul(r21905, r21902, r21904, MPFR_RNDN); mpfr_pow(r21906, r21905, r21895, MPFR_RNDN); ; mpfr_div(r21908, r21906, r21907, MPFR_RNDN); mpfr_mul(r21909, r21901, r21908, MPFR_RNDN); mpfr_add(r21910, r21895, r21909, MPFR_RNDN); return mpfr_get_d(r21910, MPFR_RNDN); }