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

char *name = "fma_test2";

double f_if(float t) {
        float r25514 = 1.7e+308;
        float r25515 = t;
        float r25516 = r25514 * r25515;
        float r25517 = r25516 - r25514;
        return r25517;
}

double f_id(double t) {
        double r25518 = 1.7e+308;
        double r25519 = t;
        double r25520 = r25518 * r25519;
        double r25521 = r25520 - r25518;
        return r25521;
}


double f_of(float t) {
        float r25522 = 1.7e+308;
        float r25523 = t;
        float r25524 = r25522 * r25523;
        float r25525 = r25524 - r25522;
        return r25525;
}

double f_od(double t) {
        double r25526 = 1.7e+308;
        double r25527 = t;
        double r25528 = r25526 * r25527;
        double r25529 = r25528 - r25526;
        return r25529;
}

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 r25530, r25531, r25532, r25533;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r25530, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r25531);
        mpfr_init(r25532);
        mpfr_init(r25533);
}

double f_im(double t) {
        ;
        mpfr_set_d(r25531, t, MPFR_RNDN);
        mpfr_mul(r25532, r25530, r25531, MPFR_RNDN);
        mpfr_sub(r25533, r25532, r25530, MPFR_RNDN);
        return mpfr_get_d(r25533, MPFR_RNDN);
}

static mpfr_t r25534, r25535, r25536, r25537;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r25534, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r25535);
        mpfr_init(r25536);
        mpfr_init(r25537);
}

double f_fm(double t) {
        ;
        mpfr_set_d(r25535, t, MPFR_RNDN);
        mpfr_mul(r25536, r25534, r25535, MPFR_RNDN);
        mpfr_sub(r25537, r25536, r25534, MPFR_RNDN);
        return mpfr_get_d(r25537, MPFR_RNDN);
}

static mpfr_t r25538, r25539, r25540, r25541;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r25538, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r25539);
        mpfr_init(r25540);
        mpfr_init(r25541);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r25539, t, MPFR_RNDN);
        mpfr_mul(r25540, r25538, r25539, MPFR_RNDN);
        mpfr_sub(r25541, r25540, r25538, MPFR_RNDN);
        return mpfr_get_d(r25541, MPFR_RNDN);
}