#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 r9122 = 1;
        float r9123 = t;
        float r9124 = 2e-16;
        float r9125 = r9123 * r9124;
        float r9126 = r9122 + r9125;
        float r9127 = r9126 * r9126;
        float r9128 = -1;
        float r9129 = 2;
        float r9130 = r9129 * r9125;
        float r9131 = r9128 - r9130;
        float r9132 = r9127 + r9131;
        return r9132;
}

double f_id(double t) {
        double r9133 = 1;
        double r9134 = t;
        double r9135 = 2e-16;
        double r9136 = r9134 * r9135;
        double r9137 = r9133 + r9136;
        double r9138 = r9137 * r9137;
        double r9139 = -1;
        double r9140 = 2;
        double r9141 = r9140 * r9136;
        double r9142 = r9139 - r9141;
        double r9143 = r9138 + r9142;
        return r9143;
}


double f_of(float t) {
        float r9144 = 3.9999999999999997e-32;
        float r9145 = t;
        float r9146 = r9144 * r9145;
        float r9147 = r9146 * r9145;
        return r9147;
}

double f_od(double t) {
        double r9148 = 3.9999999999999997e-32;
        double r9149 = t;
        double r9150 = r9148 * r9149;
        double r9151 = r9150 * r9149;
        return r9151;
}

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 r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159, r9160, r9161, r9162;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9152, "1", 10, MPFR_RNDN);
        mpfr_init(r9153);
        mpfr_init_set_str(r9154, "2e-16", 10, MPFR_RNDN);
        mpfr_init(r9155);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init_set_str(r9158, "-1", 10, MPFR_RNDN);
        mpfr_init_set_str(r9159, "2", 10, MPFR_RNDN);
        mpfr_init(r9160);
        mpfr_init(r9161);
        mpfr_init(r9162);
}

double f_im(double t) {
        ;
        mpfr_set_d(r9153, t, MPFR_RNDN);
        ;
        mpfr_mul(r9155, r9153, r9154, MPFR_RNDN);
        mpfr_add(r9156, r9152, r9155, MPFR_RNDN);
        mpfr_mul(r9157, r9156, r9156, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r9160, r9159, r9155, MPFR_RNDN);
        mpfr_sub(r9161, r9158, r9160, MPFR_RNDN);
        mpfr_add(r9162, r9157, r9161, MPFR_RNDN);
        return mpfr_get_d(r9162, MPFR_RNDN);
}

static mpfr_t r9163, r9164, r9165, r9166;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9163, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
}

double f_fm(double t) {
        ;
        mpfr_set_d(r9164, t, MPFR_RNDN);
        mpfr_mul(r9165, r9163, r9164, MPFR_RNDN);
        mpfr_mul(r9166, r9165, r9164, MPFR_RNDN);
        return mpfr_get_d(r9166, MPFR_RNDN);
}

static mpfr_t r9167, r9168, r9169, r9170;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9167, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init(r9170);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r9168, t, MPFR_RNDN);
        mpfr_mul(r9169, r9167, r9168, MPFR_RNDN);
        mpfr_mul(r9170, r9169, r9168, MPFR_RNDN);
        return mpfr_get_d(r9170, MPFR_RNDN);
}