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

char *name = "simple fma test";

double f_if(float x, float y, float z) {
        float r15692 = x;
        float r15693 = y;
        float r15694 = z;
        float r15695 = fma(r15692, r15693, r15694);
        float r15696 = 1.0f;
        float r15697 = r15692 * r15693;
        float r15698 = r15697 + r15694;
        float r15699 = r15696 + r15698;
        float r15700 = r15695 - r15699;
        return r15700;
}

double f_id(double x, double y, double z) {
        double r15701 = x;
        double r15702 = y;
        double r15703 = z;
        double r15704 = fma(r15701, r15702, r15703);
        double r15705 = 1.0;
        double r15706 = r15701 * r15702;
        double r15707 = r15706 + r15703;
        double r15708 = r15705 + r15707;
        double r15709 = r15704 - r15708;
        return r15709;
}


double f_of(float x, float y, float z) {
        float r15710 = z;
        float r15711 = -1.6846166027071804e+162f;
        bool r15712 = r15710 <= r15711;
        float r15713 = 1.0f;
        float r15714 = x;
        float r15715 = r15713 / r15714;
        float r15716 = y;
        float r15717 = r15713 / r15716;
        float r15718 = r15713 / r15710;
        float r15719 = fma(r15715, r15717, r15718);
        float r15720 = r15719 - r15713;
        float r15721 = 2.8499054335704732e+197f;
        bool r15722 = r15710 <= r15721;
        float r15723 = fma(r15714, r15716, r15710);
        float r15724 = r15723 - r15710;
        float r15725 = r15716 * r15714;
        float r15726 = r15724 - r15725;
        float r15727 = r15726 - r15713;
        float r15728 = exp(r15727);
        float r15729 = log(r15728);
        float r15730 = r15722 ? r15729 : r15720;
        float r15731 = r15712 ? r15720 : r15730;
        return r15731;
}

double f_od(double x, double y, double z) {
        double r15732 = z;
        double r15733 = -1.6846166027071804e+162;
        bool r15734 = r15732 <= r15733;
        double r15735 = 1.0;
        double r15736 = x;
        double r15737 = r15735 / r15736;
        double r15738 = y;
        double r15739 = r15735 / r15738;
        double r15740 = r15735 / r15732;
        double r15741 = fma(r15737, r15739, r15740);
        double r15742 = r15741 - r15735;
        double r15743 = 2.8499054335704732e+197;
        bool r15744 = r15732 <= r15743;
        double r15745 = fma(r15736, r15738, r15732);
        double r15746 = r15745 - r15732;
        double r15747 = r15738 * r15736;
        double r15748 = r15746 - r15747;
        double r15749 = r15748 - r15735;
        double r15750 = exp(r15749);
        double r15751 = log(r15750);
        double r15752 = r15744 ? r15751 : r15742;
        double r15753 = r15734 ? r15742 : r15752;
        return r15753;
}

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 r15754, r15755, r15756, r15757, r15758, r15759, r15760, r15761, r15762;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r15754);
        mpfr_init(r15755);
        mpfr_init(r15756);
        mpfr_init(r15757);
        mpfr_init_set_str(r15758, "1", 10, MPFR_RNDN);
        mpfr_init(r15759);
        mpfr_init(r15760);
        mpfr_init(r15761);
        mpfr_init(r15762);
}

double f_im(double x, double y, double z) {
        mpfr_set_d(r15754, x, MPFR_RNDN);
        mpfr_set_d(r15755, y, MPFR_RNDN);
        mpfr_set_d(r15756, z, MPFR_RNDN);
        mpfr_fma(r15757, r15754, r15755, r15756, MPFR_RNDN);
        ;
        mpfr_mul(r15759, r15754, r15755, MPFR_RNDN);
        mpfr_add(r15760, r15759, r15756, MPFR_RNDN);
        mpfr_add(r15761, r15758, r15760, MPFR_RNDN);
        mpfr_sub(r15762, r15757, r15761, MPFR_RNDN);
        return mpfr_get_d(r15762, MPFR_RNDN);
}

static mpfr_t r15763, r15764, r15765, r15766, r15767, r15768, r15769, r15770, r15771, r15772, r15773, r15774, r15775, r15776, r15777, r15778, r15779, r15780, r15781, r15782, r15783, r15784;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r15763);
        mpfr_init_set_str(r15764, "-1.6846166027071804e+162", 10, MPFR_RNDN);
        mpfr_init(r15765);
        mpfr_init_set_str(r15766, "1", 10, MPFR_RNDN);
        mpfr_init(r15767);
        mpfr_init(r15768);
        mpfr_init(r15769);
        mpfr_init(r15770);
        mpfr_init(r15771);
        mpfr_init(r15772);
        mpfr_init(r15773);
        mpfr_init_set_str(r15774, "2.8499054335704732e+197", 10, MPFR_RNDN);
        mpfr_init(r15775);
        mpfr_init(r15776);
        mpfr_init(r15777);
        mpfr_init(r15778);
        mpfr_init(r15779);
        mpfr_init(r15780);
        mpfr_init(r15781);
        mpfr_init(r15782);
        mpfr_init(r15783);
        mpfr_init(r15784);
}

double f_fm(double x, double y, double z) {
        mpfr_set_d(r15763, z, MPFR_RNDN);
        ;
        mpfr_set_si(r15765, mpfr_cmp(r15763, r15764) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r15767, x, MPFR_RNDN);
        mpfr_div(r15768, r15766, r15767, MPFR_RNDN);
        mpfr_set_d(r15769, y, MPFR_RNDN);
        mpfr_div(r15770, r15766, r15769, MPFR_RNDN);
        mpfr_div(r15771, r15766, r15763, MPFR_RNDN);
        mpfr_fma(r15772, r15768, r15770, r15771, MPFR_RNDN);
        mpfr_sub(r15773, r15772, r15766, MPFR_RNDN);
        ;
        mpfr_set_si(r15775, mpfr_cmp(r15763, r15774) <= 0, MPFR_RNDN);
        mpfr_fma(r15776, r15767, r15769, r15763, MPFR_RNDN);
        mpfr_sub(r15777, r15776, r15763, MPFR_RNDN);
        mpfr_mul(r15778, r15769, r15767, MPFR_RNDN);
        mpfr_sub(r15779, r15777, r15778, MPFR_RNDN);
        mpfr_sub(r15780, r15779, r15766, MPFR_RNDN);
        mpfr_exp(r15781, r15780, MPFR_RNDN);
        mpfr_log(r15782, r15781, MPFR_RNDN);
        if (mpfr_get_si(r15775, MPFR_RNDN)) { mpfr_set(r15783, r15782, MPFR_RNDN); } else { mpfr_set(r15783, r15773, MPFR_RNDN); };
        if (mpfr_get_si(r15765, MPFR_RNDN)) { mpfr_set(r15784, r15773, MPFR_RNDN); } else { mpfr_set(r15784, r15783, MPFR_RNDN); };
        return mpfr_get_d(r15784, MPFR_RNDN);
}

static mpfr_t r15785, r15786, r15787, r15788, r15789, r15790, r15791, r15792, r15793, r15794, r15795, r15796, r15797, r15798, r15799, r15800, r15801, r15802, r15803, r15804, r15805, r15806;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r15785);
        mpfr_init_set_str(r15786, "-1.6846166027071804e+162", 10, MPFR_RNDN);
        mpfr_init(r15787);
        mpfr_init_set_str(r15788, "1", 10, MPFR_RNDN);
        mpfr_init(r15789);
        mpfr_init(r15790);
        mpfr_init(r15791);
        mpfr_init(r15792);
        mpfr_init(r15793);
        mpfr_init(r15794);
        mpfr_init(r15795);
        mpfr_init_set_str(r15796, "2.8499054335704732e+197", 10, MPFR_RNDN);
        mpfr_init(r15797);
        mpfr_init(r15798);
        mpfr_init(r15799);
        mpfr_init(r15800);
        mpfr_init(r15801);
        mpfr_init(r15802);
        mpfr_init(r15803);
        mpfr_init(r15804);
        mpfr_init(r15805);
        mpfr_init(r15806);
}

double f_dm(double x, double y, double z) {
        mpfr_set_d(r15785, z, MPFR_RNDN);
        ;
        mpfr_set_si(r15787, mpfr_cmp(r15785, r15786) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r15789, x, MPFR_RNDN);
        mpfr_div(r15790, r15788, r15789, MPFR_RNDN);
        mpfr_set_d(r15791, y, MPFR_RNDN);
        mpfr_div(r15792, r15788, r15791, MPFR_RNDN);
        mpfr_div(r15793, r15788, r15785, MPFR_RNDN);
        mpfr_fma(r15794, r15790, r15792, r15793, MPFR_RNDN);
        mpfr_sub(r15795, r15794, r15788, MPFR_RNDN);
        ;
        mpfr_set_si(r15797, mpfr_cmp(r15785, r15796) <= 0, MPFR_RNDN);
        mpfr_fma(r15798, r15789, r15791, r15785, MPFR_RNDN);
        mpfr_sub(r15799, r15798, r15785, MPFR_RNDN);
        mpfr_mul(r15800, r15791, r15789, MPFR_RNDN);
        mpfr_sub(r15801, r15799, r15800, MPFR_RNDN);
        mpfr_sub(r15802, r15801, r15788, MPFR_RNDN);
        mpfr_exp(r15803, r15802, MPFR_RNDN);
        mpfr_log(r15804, r15803, MPFR_RNDN);
        if (mpfr_get_si(r15797, MPFR_RNDN)) { mpfr_set(r15805, r15804, MPFR_RNDN); } else { mpfr_set(r15805, r15795, MPFR_RNDN); };
        if (mpfr_get_si(r15787, MPFR_RNDN)) { mpfr_set(r15806, r15795, MPFR_RNDN); } else { mpfr_set(r15806, r15805, MPFR_RNDN); };
        return mpfr_get_d(r15806, MPFR_RNDN);
}