#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 r21761 = 1;
        float r21762 = 6;
        float r21763 = r21761 / r21762;
        float r21764 = -2;
        float r21765 = u1;
        float r21766 = log(r21765);
        float r21767 = r21764 * r21766;
        float r21768 = 0.5;
        float r21769 = pow(r21767, r21768);
        float r21770 = r21763 * r21769;
        float r21771 = 2;
        float r21772 = atan2(1.0, 0.0);
        float r21773 = r21771 * r21772;
        float r21774 = u2;
        float r21775 = r21773 * r21774;
        float r21776 = cos(r21775);
        float r21777 = r21770 * r21776;
        float r21778 = r21777 + r21768;
        return r21778;
}

double f_id(double u1, double u2) {
        double r21779 = 1;
        double r21780 = 6;
        double r21781 = r21779 / r21780;
        double r21782 = -2;
        double r21783 = u1;
        double r21784 = log(r21783);
        double r21785 = r21782 * r21784;
        double r21786 = 0.5;
        double r21787 = pow(r21785, r21786);
        double r21788 = r21781 * r21787;
        double r21789 = 2;
        double r21790 = atan2(1.0, 0.0);
        double r21791 = r21789 * r21790;
        double r21792 = u2;
        double r21793 = r21791 * r21792;
        double r21794 = cos(r21793);
        double r21795 = r21788 * r21794;
        double r21796 = r21795 + r21786;
        return r21796;
}


double f_of(float u1, float u2) {
        float r21797 = atan2(1.0, 0.0);
        float r21798 = u2;
        float r21799 = r21797 * r21798;
        float r21800 = 2;
        float r21801 = r21799 * r21800;
        float r21802 = cos(r21801);
        float r21803 = 1;
        float r21804 = 6;
        float r21805 = -2;
        float r21806 = u1;
        float r21807 = log(r21806);
        float r21808 = r21805 * r21807;
        float r21809 = 0.5;
        float r21810 = pow(r21808, r21809);
        float r21811 = r21804 / r21810;
        float r21812 = r21803 / r21811;
        float r21813 = fma(r21802, r21812, r21809);
        return r21813;
}

double f_od(double u1, double u2) {
        double r21814 = atan2(1.0, 0.0);
        double r21815 = u2;
        double r21816 = r21814 * r21815;
        double r21817 = 2;
        double r21818 = r21816 * r21817;
        double r21819 = cos(r21818);
        double r21820 = 1;
        double r21821 = 6;
        double r21822 = -2;
        double r21823 = u1;
        double r21824 = log(r21823);
        double r21825 = r21822 * r21824;
        double r21826 = 0.5;
        double r21827 = pow(r21825, r21826);
        double r21828 = r21821 / r21827;
        double r21829 = r21820 / r21828;
        double r21830 = fma(r21819, r21829, r21826);
        return r21830;
}

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 r21831, r21832, r21833, r21834, r21835, r21836, r21837, r21838, r21839, r21840, r21841, r21842, r21843, r21844, r21845, r21846, r21847, r21848;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r21831, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r21832, "6", 10, MPFR_RNDN);
        mpfr_init(r21833);
        mpfr_init_set_str(r21834, "-2", 10, MPFR_RNDN);
        mpfr_init(r21835);
        mpfr_init(r21836);
        mpfr_init(r21837);
        mpfr_init_set_str(r21838, "0.5", 10, MPFR_RNDN);
        mpfr_init(r21839);
        mpfr_init(r21840);
        mpfr_init_set_str(r21841, "2", 10, MPFR_RNDN);
        mpfr_init(r21842);
        mpfr_init(r21843);
        mpfr_init(r21844);
        mpfr_init(r21845);
        mpfr_init(r21846);
        mpfr_init(r21847);
        mpfr_init(r21848);
}

double f_im(double u1, double u2) {
        ;
        ;
        mpfr_div(r21833, r21831, r21832, MPFR_RNDN);
        ;
        mpfr_set_d(r21835, u1, MPFR_RNDN);
        mpfr_log(r21836, r21835, MPFR_RNDN);
        mpfr_mul(r21837, r21834, r21836, MPFR_RNDN);
        ;
        mpfr_pow(r21839, r21837, r21838, MPFR_RNDN);
        mpfr_mul(r21840, r21833, r21839, MPFR_RNDN);
        ;
        mpfr_const_pi(r21842, MPFR_RNDN);
        mpfr_mul(r21843, r21841, r21842, MPFR_RNDN);
        mpfr_set_d(r21844, u2, MPFR_RNDN);
        mpfr_mul(r21845, r21843, r21844, MPFR_RNDN);
        mpfr_cos(r21846, r21845, MPFR_RNDN);
        mpfr_mul(r21847, r21840, r21846, MPFR_RNDN);
        mpfr_add(r21848, r21847, r21838, MPFR_RNDN);
        return mpfr_get_d(r21848, MPFR_RNDN);
}

static mpfr_t r21849, r21850, r21851, r21852, r21853, r21854, r21855, r21856, r21857, r21858, r21859, r21860, r21861, r21862, r21863, r21864, r21865;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r21849);
        mpfr_init(r21850);
        mpfr_init(r21851);
        mpfr_init_set_str(r21852, "2", 10, MPFR_RNDN);
        mpfr_init(r21853);
        mpfr_init(r21854);
        mpfr_init_set_str(r21855, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r21856, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r21857, "-2", 10, MPFR_RNDN);
        mpfr_init(r21858);
        mpfr_init(r21859);
        mpfr_init(r21860);
        mpfr_init_set_str(r21861, "0.5", 10, MPFR_RNDN);
        mpfr_init(r21862);
        mpfr_init(r21863);
        mpfr_init(r21864);
        mpfr_init(r21865);
}

double f_fm(double u1, double u2) {
        mpfr_const_pi(r21849, MPFR_RNDN);
        mpfr_set_d(r21850, u2, MPFR_RNDN);
        mpfr_mul(r21851, r21849, r21850, MPFR_RNDN);
        ;
        mpfr_mul(r21853, r21851, r21852, MPFR_RNDN);
        mpfr_cos(r21854, r21853, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r21858, u1, MPFR_RNDN);
        mpfr_log(r21859, r21858, MPFR_RNDN);
        mpfr_mul(r21860, r21857, r21859, MPFR_RNDN);
        ;
        mpfr_pow(r21862, r21860, r21861, MPFR_RNDN);
        mpfr_div(r21863, r21856, r21862, MPFR_RNDN);
        mpfr_div(r21864, r21855, r21863, MPFR_RNDN);
        mpfr_fma(r21865, r21854, r21864, r21861, MPFR_RNDN);
        return mpfr_get_d(r21865, MPFR_RNDN);
}

static mpfr_t r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878, r21879, r21880, r21881, r21882;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r21866);
        mpfr_init(r21867);
        mpfr_init(r21868);
        mpfr_init_set_str(r21869, "2", 10, MPFR_RNDN);
        mpfr_init(r21870);
        mpfr_init(r21871);
        mpfr_init_set_str(r21872, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r21873, "6", 10, MPFR_RNDN);
        mpfr_init_set_str(r21874, "-2", 10, MPFR_RNDN);
        mpfr_init(r21875);
        mpfr_init(r21876);
        mpfr_init(r21877);
        mpfr_init_set_str(r21878, "0.5", 10, MPFR_RNDN);
        mpfr_init(r21879);
        mpfr_init(r21880);
        mpfr_init(r21881);
        mpfr_init(r21882);
}

double f_dm(double u1, double u2) {
        mpfr_const_pi(r21866, MPFR_RNDN);
        mpfr_set_d(r21867, u2, MPFR_RNDN);
        mpfr_mul(r21868, r21866, r21867, MPFR_RNDN);
        ;
        mpfr_mul(r21870, r21868, r21869, MPFR_RNDN);
        mpfr_cos(r21871, r21870, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_set_d(r21875, u1, MPFR_RNDN);
        mpfr_log(r21876, r21875, MPFR_RNDN);
        mpfr_mul(r21877, r21874, r21876, MPFR_RNDN);
        ;
        mpfr_pow(r21879, r21877, r21878, MPFR_RNDN);
        mpfr_div(r21880, r21873, r21879, MPFR_RNDN);
        mpfr_div(r21881, r21872, r21880, MPFR_RNDN);
        mpfr_fma(r21882, r21871, r21881, r21878, MPFR_RNDN);
        return mpfr_get_d(r21882, MPFR_RNDN);
}