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

char *name = "fma_test1";

double f_if(float t) {
        float r8741 = 1.0f;
        float r8742 = t;
        float r8743 = 2e-16f;
        float r8744 = r8742 * r8743;
        float r8745 = r8741 + r8744;
        float r8746 = r8745 * r8745;
        float r8747 = -1.0f;
        float r8748 = 2.0f;
        float r8749 = r8748 * r8744;
        float r8750 = r8747 - r8749;
        float r8751 = r8746 + r8750;
        return r8751;
}

double f_id(double t) {
        double r8752 = 1.0;
        double r8753 = t;
        double r8754 = 2e-16;
        double r8755 = r8753 * r8754;
        double r8756 = r8752 + r8755;
        double r8757 = r8756 * r8756;
        double r8758 = -1.0;
        double r8759 = 2.0;
        double r8760 = r8759 * r8755;
        double r8761 = r8758 - r8760;
        double r8762 = r8757 + r8761;
        return r8762;
}


double f_of(float t) {
        float r8763 = t;
        float r8764 = 3.9999999999999997e-32f;
        float r8765 = r8763 * r8764;
        float r8766 = r8763 * r8765;
        return r8766;
}

double f_od(double t) {
        double r8767 = t;
        double r8768 = 3.9999999999999997e-32;
        double r8769 = r8767 * r8768;
        double r8770 = r8767 * r8769;
        return r8770;
}

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 r8771, r8772, r8773, r8774, r8775, r8776, r8777, r8778, r8779, r8780, r8781;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8771, "1", 10, MPFR_RNDN);
        mpfr_init(r8772);
        mpfr_init_set_str(r8773, "2e-16", 10, MPFR_RNDN);
        mpfr_init(r8774);
        mpfr_init(r8775);
        mpfr_init(r8776);
        mpfr_init_set_str(r8777, "-1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8778, "2", 10, MPFR_RNDN);
        mpfr_init(r8779);
        mpfr_init(r8780);
        mpfr_init(r8781);
}

double f_im(double t) {
        ;
        mpfr_set_d(r8772, t, MPFR_RNDN);
        ;
        mpfr_mul(r8774, r8772, r8773, MPFR_RNDN);
        mpfr_add(r8775, r8771, r8774, MPFR_RNDN);
        mpfr_mul(r8776, r8775, r8775, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r8779, r8778, r8774, MPFR_RNDN);
        mpfr_sub(r8780, r8777, r8779, MPFR_RNDN);
        mpfr_add(r8781, r8776, r8780, MPFR_RNDN);
        return mpfr_get_d(r8781, MPFR_RNDN);
}

static mpfr_t r8782, r8783, r8784, r8785;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8782);
        mpfr_init_set_str(r8783, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r8784);
        mpfr_init(r8785);
}

double f_fm(double t) {
        mpfr_set_d(r8782, t, MPFR_RNDN);
        ;
        mpfr_mul(r8784, r8782, r8783, MPFR_RNDN);
        mpfr_mul(r8785, r8782, r8784, MPFR_RNDN);
        return mpfr_get_d(r8785, MPFR_RNDN);
}

static mpfr_t r8786, r8787, r8788, r8789;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8786);
        mpfr_init_set_str(r8787, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r8788);
        mpfr_init(r8789);
}

double f_dm(double t) {
        mpfr_set_d(r8786, t, MPFR_RNDN);
        ;
        mpfr_mul(r8788, r8786, r8787, MPFR_RNDN);
        mpfr_mul(r8789, r8786, r8788, MPFR_RNDN);
        return mpfr_get_d(r8789, MPFR_RNDN);
}