#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 r19992 = 1;
        float r19993 = 6;
        float r19994 = r19992 / r19993;
        float r19995 = -2;
        float r19996 = u1;
        float r19997 = log(r19996);
        float r19998 = r19995 * r19997;
        float r19999 = 0.5;
        float r20000 = pow(r19998, r19999);
        float r20001 = r19994 * r20000;
        float r20002 = 2;
        float r20003 = atan2(1.0, 0.0);
        float r20004 = r20002 * r20003;
        float r20005 = u2;
        float r20006 = r20004 * r20005;
        float r20007 = cos(r20006);
        float r20008 = r20001 * r20007;
        float r20009 = r20008 + r19999;
        return r20009;
}

double f_id(double u1, double u2) {
        double r20010 = 1;
        double r20011 = 6;
        double r20012 = r20010 / r20011;
        double r20013 = -2;
        double r20014 = u1;
        double r20015 = log(r20014);
        double r20016 = r20013 * r20015;
        double r20017 = 0.5;
        double r20018 = pow(r20016, r20017);
        double r20019 = r20012 * r20018;
        double r20020 = 2;
        double r20021 = atan2(1.0, 0.0);
        double r20022 = r20020 * r20021;
        double r20023 = u2;
        double r20024 = r20022 * r20023;
        double r20025 = cos(r20024);
        double r20026 = r20019 * r20025;
        double r20027 = r20026 + r20017;
        return r20027;
}


double f_of(float u1, float u2) {
        float r20028 = 0.5;
        float r20029 = u2;
        float r20030 = atan2(1.0, 0.0);
        float r20031 = r20030 + r20030;
        float r20032 = r20029 * r20031;
        float r20033 = cos(r20032);
        float r20034 = 1;
        float r20035 = 6;
        float r20036 = -2;
        float r20037 = u1;
        float r20038 = log(r20037);
        float r20039 = r20036 * r20038;
        float r20040 = pow(r20039, r20028);
        float r20041 = r20035 / r20040;
        float r20042 = r20034 / r20041;
        float r20043 = r20033 * r20042;
        float r20044 = r20028 + r20043;
        return r20044;
}

double f_od(double u1, double u2) {
        double r20045 = 0.5;
        double r20046 = u2;
        double r20047 = atan2(1.0, 0.0);
        double r20048 = r20047 + r20047;
        double r20049 = r20046 * r20048;
        double r20050 = cos(r20049);
        double r20051 = 1;
        double r20052 = 6;
        double r20053 = -2;
        double r20054 = u1;
        double r20055 = log(r20054);
        double r20056 = r20053 * r20055;
        double r20057 = pow(r20056, r20045);
        double r20058 = r20052 / r20057;
        double r20059 = r20051 / r20058;
        double r20060 = r20050 * r20059;
        double r20061 = r20045 + r20060;
        return r20061;
}

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 r20062, r20063, r20064, r20065, r20066, r20067, r20068, r20069, r20070, r20071, r20072, r20073, r20074, r20075, r20076, r20077, r20078, r20079;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20062, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20063, "6", 10, MPFR_RNDN);
        mpfr_init(r20064);
        mpfr_init_set_str(r20065, "-2", 10, MPFR_RNDN);
        mpfr_init(r20066);
        mpfr_init(r20067);
        mpfr_init(r20068);
        mpfr_init_set_str(r20069, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20070);
        mpfr_init(r20071);
        mpfr_init_set_str(r20072, "2", 10, MPFR_RNDN);
        mpfr_init(r20073);
        mpfr_init(r20074);
        mpfr_init(r20075);
        mpfr_init(r20076);
        mpfr_init(r20077);
        mpfr_init(r20078);
        mpfr_init(r20079);
}

double f_im(double u1, double u2) {
        ;
        ;
        mpfr_div(r20064, r20062, r20063, MPFR_RNDN);
        ;
        mpfr_set_d(r20066, u1, MPFR_RNDN);
        mpfr_log(r20067, r20066, MPFR_RNDN);
        mpfr_mul(r20068, r20065, r20067, MPFR_RNDN);
        ;
        mpfr_pow(r20070, r20068, r20069, MPFR_RNDN);
        mpfr_mul(r20071, r20064, r20070, MPFR_RNDN);
        ;
        mpfr_const_pi(r20073, MPFR_RNDN);
        mpfr_mul(r20074, r20072, r20073, MPFR_RNDN);
        mpfr_set_d(r20075, u2, MPFR_RNDN);
        mpfr_mul(r20076, r20074, r20075, MPFR_RNDN);
        mpfr_cos(r20077, r20076, MPFR_RNDN);
        mpfr_mul(r20078, r20071, r20077, MPFR_RNDN);
        mpfr_add(r20079, r20078, r20069, MPFR_RNDN);
        return mpfr_get_d(r20079, MPFR_RNDN);
}

static mpfr_t r20080, r20081, r20082, r20083, r20084, r20085, r20086, r20087, r20088, r20089, r20090, r20091, r20092, r20093, r20094, r20095, r20096;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20080, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20081);
        mpfr_init(r20082);
        mpfr_init(r20083);
        mpfr_init(r20084);
        mpfr_init(r20085);
        mpfr_init_set_str(r20086, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20087, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20088, "-2", 10, MPFR_RNDN);
        mpfr_init(r20089);
        mpfr_init(r20090);
        mpfr_init(r20091);
        mpfr_init(r20092);
        mpfr_init(r20093);
        mpfr_init(r20094);
        mpfr_init(r20095);
        mpfr_init(r20096);
}

double f_fm(double u1, double u2) {
        ;
        mpfr_set_d(r20081, u2, MPFR_RNDN);
        mpfr_const_pi(r20082, MPFR_RNDN);
        mpfr_add(r20083, r20082, r20082, MPFR_RNDN);
        mpfr_mul(r20084, r20081, r20083, MPFR_RNDN);
        mpfr_cos(r20085, r20084, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20089, u1, MPFR_RNDN);
        mpfr_log(r20090, r20089, MPFR_RNDN);
        mpfr_mul(r20091, r20088, r20090, MPFR_RNDN);
        mpfr_pow(r20092, r20091, r20080, MPFR_RNDN);
        mpfr_div(r20093, r20087, r20092, MPFR_RNDN);
        mpfr_div(r20094, r20086, r20093, MPFR_RNDN);
        mpfr_mul(r20095, r20085, r20094, MPFR_RNDN);
        mpfr_add(r20096, r20080, r20095, MPFR_RNDN);
        return mpfr_get_d(r20096, MPFR_RNDN);
}

static mpfr_t r20097, r20098, r20099, r20100, r20101, r20102, r20103, r20104, r20105, r20106, r20107, r20108, r20109, r20110, r20111, r20112, r20113;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20097, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20098);
        mpfr_init(r20099);
        mpfr_init(r20100);
        mpfr_init(r20101);
        mpfr_init(r20102);
        mpfr_init_set_str(r20103, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20104, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20105, "-2", 10, MPFR_RNDN);
        mpfr_init(r20106);
        mpfr_init(r20107);
        mpfr_init(r20108);
        mpfr_init(r20109);
        mpfr_init(r20110);
        mpfr_init(r20111);
        mpfr_init(r20112);
        mpfr_init(r20113);
}

double f_dm(double u1, double u2) {
        ;
        mpfr_set_d(r20098, u2, MPFR_RNDN);
        mpfr_const_pi(r20099, MPFR_RNDN);
        mpfr_add(r20100, r20099, r20099, MPFR_RNDN);
        mpfr_mul(r20101, r20098, r20100, MPFR_RNDN);
        mpfr_cos(r20102, r20101, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20106, u1, MPFR_RNDN);
        mpfr_log(r20107, r20106, MPFR_RNDN);
        mpfr_mul(r20108, r20105, r20107, MPFR_RNDN);
        mpfr_pow(r20109, r20108, r20097, MPFR_RNDN);
        mpfr_div(r20110, r20104, r20109, MPFR_RNDN);
        mpfr_div(r20111, r20103, r20110, MPFR_RNDN);
        mpfr_mul(r20112, r20102, r20111, MPFR_RNDN);
        mpfr_add(r20113, r20097, r20112, MPFR_RNDN);
        return mpfr_get_d(r20113, MPFR_RNDN);
}