#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 r22837 = 1.7e+308;
        float r22838 = t;
        float r22839 = r22837 * r22838;
        float r22840 = r22839 - r22837;
        return r22840;
}

double f_id(double t) {
        double r22841 = 1.7e+308;
        double r22842 = t;
        double r22843 = r22841 * r22842;
        double r22844 = r22843 - r22841;
        return r22844;
}


double f_of(float t) {
        float r22845 = 1.7e+308;
        float r22846 = t;
        float r22847 = r22845 * r22846;
        float r22848 = r22847 - r22845;
        return r22848;
}

double f_od(double t) {
        double r22849 = 1.7e+308;
        double r22850 = t;
        double r22851 = r22849 * r22850;
        double r22852 = r22851 - r22849;
        return r22852;
}

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 r22853, r22854, r22855, r22856;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22853, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r22854);
        mpfr_init(r22855);
        mpfr_init(r22856);
}

double f_im(double t) {
        ;
        mpfr_set_d(r22854, t, MPFR_RNDN);
        mpfr_mul(r22855, r22853, r22854, MPFR_RNDN);
        mpfr_sub(r22856, r22855, r22853, MPFR_RNDN);
        return mpfr_get_d(r22856, MPFR_RNDN);
}

static mpfr_t r22857, r22858, r22859, r22860;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22857, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r22858);
        mpfr_init(r22859);
        mpfr_init(r22860);
}

double f_fm(double t) {
        ;
        mpfr_set_d(r22858, t, MPFR_RNDN);
        mpfr_mul(r22859, r22857, r22858, MPFR_RNDN);
        mpfr_sub(r22860, r22859, r22857, MPFR_RNDN);
        return mpfr_get_d(r22860, MPFR_RNDN);
}

static mpfr_t r22861, r22862, r22863, r22864;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r22861, "1.7e+308", 10, MPFR_RNDN);
        mpfr_init(r22862);
        mpfr_init(r22863);
        mpfr_init(r22864);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r22862, t, MPFR_RNDN);
        mpfr_mul(r22863, r22861, r22862, MPFR_RNDN);
        mpfr_sub(r22864, r22863, r22861, MPFR_RNDN);
        return mpfr_get_d(r22864, MPFR_RNDN);
}