#include #include #include #include #include char *name = "normal distribution"; double f_if(float u1, float u2) { float r20788 = 1; float r20789 = 6; float r20790 = r20788 / r20789; float r20791 = -2; float r20792 = u1; float r20793 = log(r20792); float r20794 = r20791 * r20793; float r20795 = 0.5; float r20796 = pow(r20794, r20795); float r20797 = r20790 * r20796; float r20798 = 2; float r20799 = atan2(1.0, 0.0); float r20800 = r20798 * r20799; float r20801 = u2; float r20802 = r20800 * r20801; float r20803 = cos(r20802); float r20804 = r20797 * r20803; float r20805 = r20804 + r20795; return r20805; } double f_id(double u1, double u2) { double r20806 = 1; double r20807 = 6; double r20808 = r20806 / r20807; double r20809 = -2; double r20810 = u1; double r20811 = log(r20810); double r20812 = r20809 * r20811; double r20813 = 0.5; double r20814 = pow(r20812, r20813); double r20815 = r20808 * r20814; double r20816 = 2; double r20817 = atan2(1.0, 0.0); double r20818 = r20816 * r20817; double r20819 = u2; double r20820 = r20818 * r20819; double r20821 = cos(r20820); double r20822 = r20815 * r20821; double r20823 = r20822 + r20813; return r20823; } double f_of(float u1, float u2) { float r20824 = 1; float r20825 = 6; float r20826 = r20824 / r20825; float r20827 = sqrt(r20826); float r20828 = -2; float r20829 = u1; float r20830 = log(r20829); float r20831 = r20828 * r20830; float r20832 = 0.5; float r20833 = pow(r20831, r20832); float r20834 = r20827 * r20833; float r20835 = r20827 * r20834; float r20836 = 2; float r20837 = atan2(1.0, 0.0); float r20838 = r20836 * r20837; float r20839 = u2; float r20840 = r20838 * r20839; float r20841 = cos(r20840); float r20842 = r20835 * r20841; float r20843 = r20842 + r20832; return r20843; } double f_od(double u1, double u2) { double r20844 = 1; double r20845 = 6; double r20846 = r20844 / r20845; double r20847 = sqrt(r20846); double r20848 = -2; double r20849 = u1; double r20850 = log(r20849); double r20851 = r20848 * r20850; double r20852 = 0.5; double r20853 = pow(r20851, r20852); double r20854 = r20847 * r20853; double r20855 = r20847 * r20854; double r20856 = 2; double r20857 = atan2(1.0, 0.0); double r20858 = r20856 * r20857; double r20859 = u2; double r20860 = r20858 * r20859; double r20861 = cos(r20860); double r20862 = r20855 * r20861; double r20863 = r20862 + r20852; return r20863; } 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 r20864, r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872, r20873, r20874, r20875, r20876, r20877, r20878, r20879, r20880, r20881; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20864, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20865, "6", 10, MPFR_RNDN); mpfr_init(r20866); mpfr_init_set_str(r20867, "-2", 10, MPFR_RNDN); mpfr_init(r20868); mpfr_init(r20869); mpfr_init(r20870); mpfr_init_set_str(r20871, "0.5", 10, MPFR_RNDN); mpfr_init(r20872); mpfr_init(r20873); mpfr_init_set_str(r20874, "2", 10, MPFR_RNDN); mpfr_init(r20875); mpfr_init(r20876); mpfr_init(r20877); mpfr_init(r20878); mpfr_init(r20879); mpfr_init(r20880); mpfr_init(r20881); } double f_im(double u1, double u2) { ; ; mpfr_div(r20866, r20864, r20865, MPFR_RNDN); ; mpfr_set_d(r20868, u1, MPFR_RNDN); mpfr_log(r20869, r20868, MPFR_RNDN); mpfr_mul(r20870, r20867, r20869, MPFR_RNDN); ; mpfr_pow(r20872, r20870, r20871, MPFR_RNDN); mpfr_mul(r20873, r20866, r20872, MPFR_RNDN); ; mpfr_const_pi(r20875, MPFR_RNDN); mpfr_mul(r20876, r20874, r20875, MPFR_RNDN); mpfr_set_d(r20877, u2, MPFR_RNDN); mpfr_mul(r20878, r20876, r20877, MPFR_RNDN); mpfr_cos(r20879, r20878, MPFR_RNDN); mpfr_mul(r20880, r20873, r20879, MPFR_RNDN); mpfr_add(r20881, r20880, r20871, MPFR_RNDN); return mpfr_get_d(r20881, MPFR_RNDN); } static mpfr_t r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20882, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20883, "6", 10, MPFR_RNDN); mpfr_init(r20884); mpfr_init(r20885); mpfr_init_set_str(r20886, "-2", 10, MPFR_RNDN); mpfr_init(r20887); mpfr_init(r20888); mpfr_init(r20889); mpfr_init_set_str(r20890, "0.5", 10, MPFR_RNDN); mpfr_init(r20891); mpfr_init(r20892); mpfr_init(r20893); mpfr_init_set_str(r20894, "2", 10, MPFR_RNDN); mpfr_init(r20895); mpfr_init(r20896); mpfr_init(r20897); mpfr_init(r20898); mpfr_init(r20899); mpfr_init(r20900); mpfr_init(r20901); } double f_fm(double u1, double u2) { ; ; mpfr_div(r20884, r20882, r20883, MPFR_RNDN); mpfr_sqrt(r20885, r20884, MPFR_RNDN); ; mpfr_set_d(r20887, u1, MPFR_RNDN); mpfr_log(r20888, r20887, MPFR_RNDN); mpfr_mul(r20889, r20886, r20888, MPFR_RNDN); ; mpfr_pow(r20891, r20889, r20890, MPFR_RNDN); mpfr_mul(r20892, r20885, r20891, MPFR_RNDN); mpfr_mul(r20893, r20885, r20892, MPFR_RNDN); ; mpfr_const_pi(r20895, MPFR_RNDN); mpfr_mul(r20896, r20894, r20895, MPFR_RNDN); mpfr_set_d(r20897, u2, MPFR_RNDN); mpfr_mul(r20898, r20896, r20897, MPFR_RNDN); mpfr_cos(r20899, r20898, MPFR_RNDN); mpfr_mul(r20900, r20893, r20899, MPFR_RNDN); mpfr_add(r20901, r20900, r20890, MPFR_RNDN); return mpfr_get_d(r20901, MPFR_RNDN); } static mpfr_t r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20902, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20903, "6", 10, MPFR_RNDN); mpfr_init(r20904); mpfr_init(r20905); mpfr_init_set_str(r20906, "-2", 10, MPFR_RNDN); mpfr_init(r20907); mpfr_init(r20908); mpfr_init(r20909); mpfr_init_set_str(r20910, "0.5", 10, MPFR_RNDN); mpfr_init(r20911); mpfr_init(r20912); mpfr_init(r20913); mpfr_init_set_str(r20914, "2", 10, MPFR_RNDN); mpfr_init(r20915); mpfr_init(r20916); mpfr_init(r20917); mpfr_init(r20918); mpfr_init(r20919); mpfr_init(r20920); mpfr_init(r20921); } double f_dm(double u1, double u2) { ; ; mpfr_div(r20904, r20902, r20903, MPFR_RNDN); mpfr_sqrt(r20905, r20904, MPFR_RNDN); ; mpfr_set_d(r20907, u1, MPFR_RNDN); mpfr_log(r20908, r20907, MPFR_RNDN); mpfr_mul(r20909, r20906, r20908, MPFR_RNDN); ; mpfr_pow(r20911, r20909, r20910, MPFR_RNDN); mpfr_mul(r20912, r20905, r20911, MPFR_RNDN); mpfr_mul(r20913, r20905, r20912, MPFR_RNDN); ; mpfr_const_pi(r20915, MPFR_RNDN); mpfr_mul(r20916, r20914, r20915, MPFR_RNDN); mpfr_set_d(r20917, u2, MPFR_RNDN); mpfr_mul(r20918, r20916, r20917, MPFR_RNDN); mpfr_cos(r20919, r20918, MPFR_RNDN); mpfr_mul(r20920, r20913, r20919, MPFR_RNDN); mpfr_add(r20921, r20920, r20910, MPFR_RNDN); return mpfr_get_d(r20921, MPFR_RNDN); }