#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "normal distribution";

double f_if(float u1, float u2) {
        float r20814 = 1;
        float r20815 = 6;
        float r20816 = r20814 / r20815;
        float r20817 = -2;
        float r20818 = u1;
        float r20819 = log(r20818);
        float r20820 = r20817 * r20819;
        float r20821 = 0.5;
        float r20822 = pow(r20820, r20821);
        float r20823 = r20816 * r20822;
        float r20824 = 2;
        float r20825 = atan2(1.0, 0.0);
        float r20826 = r20824 * r20825;
        float r20827 = u2;
        float r20828 = r20826 * r20827;
        float r20829 = cos(r20828);
        float r20830 = r20823 * r20829;
        float r20831 = r20830 + r20821;
        return r20831;
}

double f_id(double u1, double u2) {
        double r20832 = 1;
        double r20833 = 6;
        double r20834 = r20832 / r20833;
        double r20835 = -2;
        double r20836 = u1;
        double r20837 = log(r20836);
        double r20838 = r20835 * r20837;
        double r20839 = 0.5;
        double r20840 = pow(r20838, r20839);
        double r20841 = r20834 * r20840;
        double r20842 = 2;
        double r20843 = atan2(1.0, 0.0);
        double r20844 = r20842 * r20843;
        double r20845 = u2;
        double r20846 = r20844 * r20845;
        double r20847 = cos(r20846);
        double r20848 = r20841 * r20847;
        double r20849 = r20848 + r20839;
        return r20849;
}


double f_of(float u1, float u2) {
        float r20850 = 1;
        float r20851 = 6;
        float r20852 = r20850 / r20851;
        float r20853 = -2;
        float r20854 = u1;
        float r20855 = log(r20854);
        float r20856 = r20853 * r20855;
        float r20857 = 0.5;
        float r20858 = pow(r20856, r20857);
        float r20859 = r20852 * r20858;
        float r20860 = 2;
        float r20861 = atan2(1.0, 0.0);
        float r20862 = r20860 * r20861;
        float r20863 = u2;
        float r20864 = r20862 * r20863;
        float r20865 = cos(r20864);
        float r20866 = r20859 * r20865;
        float r20867 = r20866 + r20857;
        return r20867;
}

double f_od(double u1, double u2) {
        double r20868 = 1;
        double r20869 = 6;
        double r20870 = r20868 / r20869;
        double r20871 = -2;
        double r20872 = u1;
        double r20873 = log(r20872);
        double r20874 = r20871 * r20873;
        double r20875 = 0.5;
        double r20876 = pow(r20874, r20875);
        double r20877 = r20870 * r20876;
        double r20878 = 2;
        double r20879 = atan2(1.0, 0.0);
        double r20880 = r20878 * r20879;
        double r20881 = u2;
        double r20882 = r20880 * r20881;
        double r20883 = cos(r20882);
        double r20884 = r20877 * r20883;
        double r20885 = r20884 + r20875;
        return r20885;
}

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 r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20886, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20887, "6", 10, MPFR_RNDN);
        mpfr_init(r20888);
        mpfr_init_set_str(r20889, "-2", 10, MPFR_RNDN);
        mpfr_init(r20890);
        mpfr_init(r20891);
        mpfr_init(r20892);
        mpfr_init_set_str(r20893, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20894);
        mpfr_init(r20895);
        mpfr_init_set_str(r20896, "2", 10, MPFR_RNDN);
        mpfr_init(r20897);
        mpfr_init(r20898);
        mpfr_init(r20899);
        mpfr_init(r20900);
        mpfr_init(r20901);
        mpfr_init(r20902);
        mpfr_init(r20903);
}

double f_im(double u1, double u2) {
        ;
        ;
        mpfr_div(r20888, r20886, r20887, MPFR_RNDN);
        ;
        mpfr_set_d(r20890, u1, MPFR_RNDN);
        mpfr_log(r20891, r20890, MPFR_RNDN);
        mpfr_mul(r20892, r20889, r20891, MPFR_RNDN);
        ;
        mpfr_pow(r20894, r20892, r20893, MPFR_RNDN);
        mpfr_mul(r20895, r20888, r20894, MPFR_RNDN);
        ;
        mpfr_const_pi(r20897, MPFR_RNDN);
        mpfr_mul(r20898, r20896, r20897, MPFR_RNDN);
        mpfr_set_d(r20899, u2, MPFR_RNDN);
        mpfr_mul(r20900, r20898, r20899, MPFR_RNDN);
        mpfr_cos(r20901, r20900, MPFR_RNDN);
        mpfr_mul(r20902, r20895, r20901, MPFR_RNDN);
        mpfr_add(r20903, r20902, r20893, MPFR_RNDN);
        return mpfr_get_d(r20903, MPFR_RNDN);
}

static mpfr_t r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20904, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20905, "6", 10, MPFR_RNDN);
        mpfr_init(r20906);
        mpfr_init_set_str(r20907, "-2", 10, MPFR_RNDN);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init(r20910);
        mpfr_init_set_str(r20911, "0.5", 10, MPFR_RNDN);
        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_fm(double u1, double u2) {
        ;
        ;
        mpfr_div(r20906, r20904, r20905, MPFR_RNDN);
        ;
        mpfr_set_d(r20908, u1, MPFR_RNDN);
        mpfr_log(r20909, r20908, MPFR_RNDN);
        mpfr_mul(r20910, r20907, r20909, MPFR_RNDN);
        ;
        mpfr_pow(r20912, r20910, r20911, MPFR_RNDN);
        mpfr_mul(r20913, r20906, 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, r20911, MPFR_RNDN);
        return mpfr_get_d(r20921, MPFR_RNDN);
}

static mpfr_t r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937, r20938, r20939;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20922, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20923, "6", 10, MPFR_RNDN);
        mpfr_init(r20924);
        mpfr_init_set_str(r20925, "-2", 10, MPFR_RNDN);
        mpfr_init(r20926);
        mpfr_init(r20927);
        mpfr_init(r20928);
        mpfr_init_set_str(r20929, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20930);
        mpfr_init(r20931);
        mpfr_init_set_str(r20932, "2", 10, MPFR_RNDN);
        mpfr_init(r20933);
        mpfr_init(r20934);
        mpfr_init(r20935);
        mpfr_init(r20936);
        mpfr_init(r20937);
        mpfr_init(r20938);
        mpfr_init(r20939);
}

double f_dm(double u1, double u2) {
        ;
        ;
        mpfr_div(r20924, r20922, r20923, MPFR_RNDN);
        ;
        mpfr_set_d(r20926, u1, MPFR_RNDN);
        mpfr_log(r20927, r20926, MPFR_RNDN);
        mpfr_mul(r20928, r20925, r20927, MPFR_RNDN);
        ;
        mpfr_pow(r20930, r20928, r20929, MPFR_RNDN);
        mpfr_mul(r20931, r20924, r20930, MPFR_RNDN);
        ;
        mpfr_const_pi(r20933, MPFR_RNDN);
        mpfr_mul(r20934, r20932, r20933, MPFR_RNDN);
        mpfr_set_d(r20935, u2, MPFR_RNDN);
        mpfr_mul(r20936, r20934, r20935, MPFR_RNDN);
        mpfr_cos(r20937, r20936, MPFR_RNDN);
        mpfr_mul(r20938, r20931, r20937, MPFR_RNDN);
        mpfr_add(r20939, r20938, r20929, MPFR_RNDN);
        return mpfr_get_d(r20939, MPFR_RNDN);
}