#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 r9304 = 1;
        float r9305 = t;
        float r9306 = 2e-16;
        float r9307 = r9305 * r9306;
        float r9308 = r9304 + r9307;
        float r9309 = r9308 * r9308;
        float r9310 = -1;
        float r9311 = 2;
        float r9312 = r9311 * r9307;
        float r9313 = r9310 - r9312;
        float r9314 = r9309 + r9313;
        return r9314;
}

double f_id(double t) {
        double r9315 = 1;
        double r9316 = t;
        double r9317 = 2e-16;
        double r9318 = r9316 * r9317;
        double r9319 = r9315 + r9318;
        double r9320 = r9319 * r9319;
        double r9321 = -1;
        double r9322 = 2;
        double r9323 = r9322 * r9318;
        double r9324 = r9321 - r9323;
        double r9325 = r9320 + r9324;
        return r9325;
}


double f_of(float t) {
        float r9326 = 3.9999999999999997e-32;
        float r9327 = t;
        float r9328 = r9326 * r9327;
        float r9329 = r9328 * r9327;
        return r9329;
}

double f_od(double t) {
        double r9330 = 3.9999999999999997e-32;
        double r9331 = t;
        double r9332 = r9330 * r9331;
        double r9333 = r9332 * r9331;
        return r9333;
}

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 r9334, r9335, r9336, r9337, r9338, r9339, r9340, r9341, r9342, r9343, r9344;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9334, "1", 10, MPFR_RNDN);
        mpfr_init(r9335);
        mpfr_init_set_str(r9336, "2e-16", 10, MPFR_RNDN);
        mpfr_init(r9337);
        mpfr_init(r9338);
        mpfr_init(r9339);
        mpfr_init_set_str(r9340, "-1", 10, MPFR_RNDN);
        mpfr_init_set_str(r9341, "2", 10, MPFR_RNDN);
        mpfr_init(r9342);
        mpfr_init(r9343);
        mpfr_init(r9344);
}

double f_im(double t) {
        ;
        mpfr_set_d(r9335, t, MPFR_RNDN);
        ;
        mpfr_mul(r9337, r9335, r9336, MPFR_RNDN);
        mpfr_add(r9338, r9334, r9337, MPFR_RNDN);
        mpfr_mul(r9339, r9338, r9338, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r9342, r9341, r9337, MPFR_RNDN);
        mpfr_sub(r9343, r9340, r9342, MPFR_RNDN);
        mpfr_add(r9344, r9339, r9343, MPFR_RNDN);
        return mpfr_get_d(r9344, MPFR_RNDN);
}

static mpfr_t r9345, r9346, r9347, r9348;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9345, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r9346);
        mpfr_init(r9347);
        mpfr_init(r9348);
}

double f_fm(double t) {
        ;
        mpfr_set_d(r9346, t, MPFR_RNDN);
        mpfr_mul(r9347, r9345, r9346, MPFR_RNDN);
        mpfr_mul(r9348, r9347, r9346, MPFR_RNDN);
        return mpfr_get_d(r9348, MPFR_RNDN);
}

static mpfr_t r9349, r9350, r9351, r9352;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9349, "3.9999999999999997e-32", 10, MPFR_RNDN);
        mpfr_init(r9350);
        mpfr_init(r9351);
        mpfr_init(r9352);
}

double f_dm(double t) {
        ;
        mpfr_set_d(r9350, t, MPFR_RNDN);
        mpfr_mul(r9351, r9349, r9350, MPFR_RNDN);
        mpfr_mul(r9352, r9351, r9350, MPFR_RNDN);
        return mpfr_get_d(r9352, MPFR_RNDN);
}