#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 r19964 = 1;
        float r19965 = 6;
        float r19966 = r19964 / r19965;
        float r19967 = -2;
        float r19968 = u1;
        float r19969 = log(r19968);
        float r19970 = r19967 * r19969;
        float r19971 = 0.5;
        float r19972 = pow(r19970, r19971);
        float r19973 = r19966 * r19972;
        float r19974 = 2;
        float r19975 = atan2(1.0, 0.0);
        float r19976 = r19974 * r19975;
        float r19977 = u2;
        float r19978 = r19976 * r19977;
        float r19979 = cos(r19978);
        float r19980 = r19973 * r19979;
        float r19981 = r19980 + r19971;
        return r19981;
}

double f_id(double u1, double u2) {
        double r19982 = 1;
        double r19983 = 6;
        double r19984 = r19982 / r19983;
        double r19985 = -2;
        double r19986 = u1;
        double r19987 = log(r19986);
        double r19988 = r19985 * r19987;
        double r19989 = 0.5;
        double r19990 = pow(r19988, r19989);
        double r19991 = r19984 * r19990;
        double r19992 = 2;
        double r19993 = atan2(1.0, 0.0);
        double r19994 = r19992 * r19993;
        double r19995 = u2;
        double r19996 = r19994 * r19995;
        double r19997 = cos(r19996);
        double r19998 = r19991 * r19997;
        double r19999 = r19998 + r19989;
        return r19999;
}


double f_of(float u1, float u2) {
        float r20000 = 0.5;
        float r20001 = u2;
        float r20002 = atan2(1.0, 0.0);
        float r20003 = r20002 + r20002;
        float r20004 = r20001 * r20003;
        float r20005 = cos(r20004);
        float r20006 = 1;
        float r20007 = 6;
        float r20008 = -2;
        float r20009 = u1;
        float r20010 = log(r20009);
        float r20011 = r20008 * r20010;
        float r20012 = pow(r20011, r20000);
        float r20013 = r20007 / r20012;
        float r20014 = r20006 / r20013;
        float r20015 = r20005 * r20014;
        float r20016 = r20000 + r20015;
        return r20016;
}

double f_od(double u1, double u2) {
        double r20017 = 0.5;
        double r20018 = u2;
        double r20019 = atan2(1.0, 0.0);
        double r20020 = r20019 + r20019;
        double r20021 = r20018 * r20020;
        double r20022 = cos(r20021);
        double r20023 = 1;
        double r20024 = 6;
        double r20025 = -2;
        double r20026 = u1;
        double r20027 = log(r20026);
        double r20028 = r20025 * r20027;
        double r20029 = pow(r20028, r20017);
        double r20030 = r20024 / r20029;
        double r20031 = r20023 / r20030;
        double r20032 = r20022 * r20031;
        double r20033 = r20017 + r20032;
        return r20033;
}

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 r20034, r20035, r20036, r20037, r20038, r20039, r20040, r20041, r20042, r20043, r20044, r20045, r20046, r20047, r20048, r20049, r20050, r20051;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20034, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20035, "6", 10, MPFR_RNDN);
        mpfr_init(r20036);
        mpfr_init_set_str(r20037, "-2", 10, MPFR_RNDN);
        mpfr_init(r20038);
        mpfr_init(r20039);
        mpfr_init(r20040);
        mpfr_init_set_str(r20041, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20042);
        mpfr_init(r20043);
        mpfr_init_set_str(r20044, "2", 10, MPFR_RNDN);
        mpfr_init(r20045);
        mpfr_init(r20046);
        mpfr_init(r20047);
        mpfr_init(r20048);
        mpfr_init(r20049);
        mpfr_init(r20050);
        mpfr_init(r20051);
}

double f_im(double u1, double u2) {
        ;
        ;
        mpfr_div(r20036, r20034, r20035, MPFR_RNDN);
        ;
        mpfr_set_d(r20038, u1, MPFR_RNDN);
        mpfr_log(r20039, r20038, MPFR_RNDN);
        mpfr_mul(r20040, r20037, r20039, MPFR_RNDN);
        ;
        mpfr_pow(r20042, r20040, r20041, MPFR_RNDN);
        mpfr_mul(r20043, r20036, r20042, MPFR_RNDN);
        ;
        mpfr_const_pi(r20045, MPFR_RNDN);
        mpfr_mul(r20046, r20044, r20045, MPFR_RNDN);
        mpfr_set_d(r20047, u2, MPFR_RNDN);
        mpfr_mul(r20048, r20046, r20047, MPFR_RNDN);
        mpfr_cos(r20049, r20048, MPFR_RNDN);
        mpfr_mul(r20050, r20043, r20049, MPFR_RNDN);
        mpfr_add(r20051, r20050, r20041, MPFR_RNDN);
        return mpfr_get_d(r20051, MPFR_RNDN);
}

static mpfr_t r20052, r20053, r20054, r20055, r20056, r20057, r20058, r20059, r20060, r20061, r20062, r20063, r20064, r20065, r20066, r20067, r20068;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20052, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20053);
        mpfr_init(r20054);
        mpfr_init(r20055);
        mpfr_init(r20056);
        mpfr_init(r20057);
        mpfr_init_set_str(r20058, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20059, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20060, "-2", 10, MPFR_RNDN);
        mpfr_init(r20061);
        mpfr_init(r20062);
        mpfr_init(r20063);
        mpfr_init(r20064);
        mpfr_init(r20065);
        mpfr_init(r20066);
        mpfr_init(r20067);
        mpfr_init(r20068);
}

double f_fm(double u1, double u2) {
        ;
        mpfr_set_d(r20053, u2, MPFR_RNDN);
        mpfr_const_pi(r20054, MPFR_RNDN);
        mpfr_add(r20055, r20054, r20054, MPFR_RNDN);
        mpfr_mul(r20056, r20053, r20055, MPFR_RNDN);
        mpfr_cos(r20057, r20056, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20061, u1, MPFR_RNDN);
        mpfr_log(r20062, r20061, MPFR_RNDN);
        mpfr_mul(r20063, r20060, r20062, MPFR_RNDN);
        mpfr_pow(r20064, r20063, r20052, MPFR_RNDN);
        mpfr_div(r20065, r20059, r20064, MPFR_RNDN);
        mpfr_div(r20066, r20058, r20065, MPFR_RNDN);
        mpfr_mul(r20067, r20057, r20066, MPFR_RNDN);
        mpfr_add(r20068, r20052, r20067, MPFR_RNDN);
        return mpfr_get_d(r20068, MPFR_RNDN);
}

static mpfr_t r20069, r20070, r20071, r20072, r20073, r20074, r20075, r20076, r20077, r20078, r20079, r20080, r20081, r20082, r20083, r20084, r20085;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20069, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20070);
        mpfr_init(r20071);
        mpfr_init(r20072);
        mpfr_init(r20073);
        mpfr_init(r20074);
        mpfr_init_set_str(r20075, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20076, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r20077, "-2", 10, MPFR_RNDN);
        mpfr_init(r20078);
        mpfr_init(r20079);
        mpfr_init(r20080);
        mpfr_init(r20081);
        mpfr_init(r20082);
        mpfr_init(r20083);
        mpfr_init(r20084);
        mpfr_init(r20085);
}

double f_dm(double u1, double u2) {
        ;
        mpfr_set_d(r20070, u2, MPFR_RNDN);
        mpfr_const_pi(r20071, MPFR_RNDN);
        mpfr_add(r20072, r20071, r20071, MPFR_RNDN);
        mpfr_mul(r20073, r20070, r20072, MPFR_RNDN);
        mpfr_cos(r20074, r20073, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r20078, u1, MPFR_RNDN);
        mpfr_log(r20079, r20078, MPFR_RNDN);
        mpfr_mul(r20080, r20077, r20079, MPFR_RNDN);
        mpfr_pow(r20081, r20080, r20069, MPFR_RNDN);
        mpfr_div(r20082, r20076, r20081, MPFR_RNDN);
        mpfr_div(r20083, r20075, r20082, MPFR_RNDN);
        mpfr_mul(r20084, r20074, r20083, MPFR_RNDN);
        mpfr_add(r20085, r20069, r20084, MPFR_RNDN);
        return mpfr_get_d(r20085, MPFR_RNDN);
}