#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 r8744 = 1.0f;
        float r8745 = t;
        float r8746 = 2e-16f;
        float r8747 = r8745 * r8746;
        float r8748 = r8744 + r8747;
        float r8749 = r8748 * r8748;
        float r8750 = -1.0f;
        float r8751 = 2.0f;
        float r8752 = r8751 * r8747;
        float r8753 = r8750 - r8752;
        float r8754 = r8749 + r8753;
        return r8754;
}

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


double f_of(float t) {
        float r8766 = 3.9999999999999997e-32f;
        float r8767 = t;
        float r8768 = r8766 * r8767;
        float r8769 = r8768 * r8767;
        return r8769;
}

double f_od(double t) {
        double r8770 = 3.9999999999999997e-32;
        double r8771 = t;
        double r8772 = r8770 * r8771;
        double r8773 = r8772 * r8771;
        return r8773;
}

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 r8774, r8775, r8776, r8777, r8778, r8779, r8780, r8781, r8782, r8783, r8784;

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

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

static mpfr_t r8785, r8786, r8787, r8788;

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

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

static mpfr_t r8789, r8790, r8791, r8792;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8789, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r8790);
        mpfr_init(r8791);
        mpfr_init(r8792);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r8790, t, MPFR_RNDN);
        mpfr_mul(r8791, r8789, r8790, MPFR_RNDN);
        mpfr_mul(r8792, r8791, r8790, MPFR_RNDN);
        return mpfr_get_d(r8792, MPFR_RNDN);
}