#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 r20786 = 1;
        float r20787 = t;
        float r20788 = 2e-16;
        float r20789 = r20787 * r20788;
        float r20790 = r20786 + r20789;
        float r20791 = r20790 * r20790;
        float r20792 = -1;
        float r20793 = 2;
        float r20794 = r20793 * r20789;
        float r20795 = r20792 - r20794;
        float r20796 = r20791 + r20795;
        return r20796;
}

double f_id(double t) {
        double r20797 = 1;
        double r20798 = t;
        double r20799 = 2e-16;
        double r20800 = r20798 * r20799;
        double r20801 = r20797 + r20800;
        double r20802 = r20801 * r20801;
        double r20803 = -1;
        double r20804 = 2;
        double r20805 = r20804 * r20800;
        double r20806 = r20803 - r20805;
        double r20807 = r20802 + r20806;
        return r20807;
}


double f_of(float t) {
        float r20808 = 3.9999999999999997e-32;
        float r20809 = t;
        float r20810 = r20808 * r20809;
        float r20811 = r20810 * r20809;
        return r20811;
}

double f_od(double t) {
        double r20812 = 3.9999999999999997e-32;
        double r20813 = t;
        double r20814 = r20812 * r20813;
        double r20815 = r20814 * r20813;
        return r20815;
}

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 r20816, r20817, r20818, r20819, r20820, r20821, r20822, r20823, r20824, r20825, r20826;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r20816, "1", 10, MPFR_RNDN);
        mpfr_init(r20817);
        mpfr_init_set_str(r20818, "2e-16", 10, MPFR_RNDN);
        mpfr_init(r20819);
        mpfr_init(r20820);
        mpfr_init(r20821);
        mpfr_init_set_str(r20822, "-1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20823, "2", 10, MPFR_RNDN);
        mpfr_init(r20824);
        mpfr_init(r20825);
        mpfr_init(r20826);
}

double f_im(double t) {
        ;
        mpfr_set_d(r20817, t, MPFR_RNDN);
        ;
        mpfr_mul(r20819, r20817, r20818, MPFR_RNDN);
        mpfr_add(r20820, r20816, r20819, MPFR_RNDN);
        mpfr_mul(r20821, r20820, r20820, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r20824, r20823, r20819, MPFR_RNDN);
        mpfr_sub(r20825, r20822, r20824, MPFR_RNDN);
        mpfr_add(r20826, r20821, r20825, MPFR_RNDN);
        return mpfr_get_d(r20826, MPFR_RNDN);
}

static mpfr_t r20827, r20828, r20829, r20830;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r20827, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r20828);
        mpfr_init(r20829);
        mpfr_init(r20830);
}

double f_fm(double t) {
        ;
        mpfr_set_d(r20828, t, MPFR_RNDN);
        mpfr_mul(r20829, r20827, r20828, MPFR_RNDN);
        mpfr_mul(r20830, r20829, r20828, MPFR_RNDN);
        return mpfr_get_d(r20830, MPFR_RNDN);
}

static mpfr_t r20831, r20832, r20833, r20834;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(656);
        mpfr_init_set_str(r20831, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r20832);
        mpfr_init(r20833);
        mpfr_init(r20834);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r20832, t, MPFR_RNDN);
        mpfr_mul(r20833, r20831, r20832, MPFR_RNDN);
        mpfr_mul(r20834, r20833, r20832, MPFR_RNDN);
        return mpfr_get_d(r20834, MPFR_RNDN);
}