#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 r20851 = 1;
        float r20852 = 6;
        float r20853 = r20851 / r20852;
        float r20854 = -2;
        float r20855 = u1;
        float r20856 = log(r20855);
        float r20857 = r20854 * r20856;
        float r20858 = 0.5;
        float r20859 = pow(r20857, r20858);
        float r20860 = r20853 * r20859;
        float r20861 = 2;
        float r20862 = atan2(1.0, 0.0);
        float r20863 = r20861 * r20862;
        float r20864 = u2;
        float r20865 = r20863 * r20864;
        float r20866 = cos(r20865);
        float r20867 = r20860 * r20866;
        float r20868 = r20867 + r20858;
        return r20868;
}

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


double f_of(float u1, float u2) {
        float r20887 = 0.5;
        float r20888 = u2;
        float r20889 = atan2(1.0, 0.0);
        float r20890 = r20889 + r20889;
        float r20891 = r20888 * r20890;
        float r20892 = cos(r20891);
        float r20893 = 1;
        float r20894 = 6;
        float r20895 = -2;
        float r20896 = u1;
        float r20897 = log(r20896);
        float r20898 = r20895 * r20897;
        float r20899 = pow(r20898, r20887);
        float r20900 = r20894 / r20899;
        float r20901 = r20893 / r20900;
        float r20902 = r20892 * r20901;
        float r20903 = r20887 + r20902;
        return r20903;
}

double f_od(double u1, double u2) {
        double r20904 = 0.5;
        double r20905 = u2;
        double r20906 = atan2(1.0, 0.0);
        double r20907 = r20906 + r20906;
        double r20908 = r20905 * r20907;
        double r20909 = cos(r20908);
        double r20910 = 1;
        double r20911 = 6;
        double r20912 = -2;
        double r20913 = u1;
        double r20914 = log(r20913);
        double r20915 = r20912 * r20914;
        double r20916 = pow(r20915, r20904);
        double r20917 = r20911 / r20916;
        double r20918 = r20910 / r20917;
        double r20919 = r20909 * r20918;
        double r20920 = r20904 + r20919;
        return r20920;
}

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 r20921, r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937, r20938;

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

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

static mpfr_t r20939, r20940, r20941, r20942, r20943, r20944, r20945, r20946, r20947, r20948, r20949, r20950, r20951, r20952, r20953, r20954, r20955;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20939, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20940);
        mpfr_init(r20941);
        mpfr_init(r20942);
        mpfr_init(r20943);
        mpfr_init(r20944);
        mpfr_init_set_str(r20945, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20946, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20947, "-2", 10, MPFR_RNDN);
        mpfr_init(r20948);
        mpfr_init(r20949);
        mpfr_init(r20950);
        mpfr_init(r20951);
        mpfr_init(r20952);
        mpfr_init(r20953);
        mpfr_init(r20954);
        mpfr_init(r20955);
}

double f_fm(double u1, double u2) {
        ;
        mpfr_set_d(r20940, u2, MPFR_RNDN);
        mpfr_const_pi(r20941, MPFR_RNDN);
        mpfr_add(r20942, r20941, r20941, MPFR_RNDN);
        mpfr_mul(r20943, r20940, r20942, MPFR_RNDN);
        mpfr_cos(r20944, r20943, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20948, u1, MPFR_RNDN);
        mpfr_log(r20949, r20948, MPFR_RNDN);
        mpfr_mul(r20950, r20947, r20949, MPFR_RNDN);
        mpfr_pow(r20951, r20950, r20939, MPFR_RNDN);
        mpfr_div(r20952, r20946, r20951, MPFR_RNDN);
        mpfr_div(r20953, r20945, r20952, MPFR_RNDN);
        mpfr_mul(r20954, r20944, r20953, MPFR_RNDN);
        mpfr_add(r20955, r20939, r20954, MPFR_RNDN);
        return mpfr_get_d(r20955, MPFR_RNDN);
}

static mpfr_t r20956, r20957, r20958, r20959, r20960, r20961, r20962, r20963, r20964, r20965, r20966, r20967, r20968, r20969, r20970, r20971, r20972;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20956, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20957);
        mpfr_init(r20958);
        mpfr_init(r20959);
        mpfr_init(r20960);
        mpfr_init(r20961);
        mpfr_init_set_str(r20962, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20963, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20964, "-2", 10, MPFR_RNDN);
        mpfr_init(r20965);
        mpfr_init(r20966);
        mpfr_init(r20967);
        mpfr_init(r20968);
        mpfr_init(r20969);
        mpfr_init(r20970);
        mpfr_init(r20971);
        mpfr_init(r20972);
}

double f_dm(double u1, double u2) {
        ;
        mpfr_set_d(r20957, u2, MPFR_RNDN);
        mpfr_const_pi(r20958, MPFR_RNDN);
        mpfr_add(r20959, r20958, r20958, MPFR_RNDN);
        mpfr_mul(r20960, r20957, r20959, MPFR_RNDN);
        mpfr_cos(r20961, r20960, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20965, u1, MPFR_RNDN);
        mpfr_log(r20966, r20965, MPFR_RNDN);
        mpfr_mul(r20967, r20964, r20966, MPFR_RNDN);
        mpfr_pow(r20968, r20967, r20956, MPFR_RNDN);
        mpfr_div(r20969, r20963, r20968, MPFR_RNDN);
        mpfr_div(r20970, r20962, r20969, MPFR_RNDN);
        mpfr_mul(r20971, r20961, r20970, MPFR_RNDN);
        mpfr_add(r20972, r20956, r20971, MPFR_RNDN);
        return mpfr_get_d(r20972, MPFR_RNDN);
}