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

char *name = "FastMath test1";

double f_if(float d) {
        float r25679 = d;
        float r25680 = 10;
        float r25681 = r25679 * r25680;
        float r25682 = 20;
        float r25683 = r25679 * r25682;
        float r25684 = r25681 + r25683;
        return r25684;
}

double f_id(double d) {
        double r25685 = d;
        double r25686 = 10;
        double r25687 = r25685 * r25686;
        double r25688 = 20;
        double r25689 = r25685 * r25688;
        double r25690 = r25687 + r25689;
        return r25690;
}


double f_of(float d) {
        float r25691 = 20;
        float r25692 = 10;
        float r25693 = r25691 + r25692;
        float r25694 = d;
        float r25695 = r25693 * r25694;
        return r25695;
}

double f_od(double d) {
        double r25696 = 20;
        double r25697 = 10;
        double r25698 = r25696 + r25697;
        double r25699 = d;
        double r25700 = r25698 * r25699;
        return r25700;
}

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 r25701, r25702, r25703, r25704, r25705, r25706;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r25701);
        mpfr_init_set_str(r25702, "10", 10, MPFR_RNDN);
        mpfr_init(r25703);
        mpfr_init_set_str(r25704, "20", 10, MPFR_RNDN);
        mpfr_init(r25705);
        mpfr_init(r25706);
}

double f_im(double d) {
        mpfr_set_d(r25701, d, MPFR_RNDN);
        ;
        mpfr_mul(r25703, r25701, r25702, MPFR_RNDN);
        ;
        mpfr_mul(r25705, r25701, r25704, MPFR_RNDN);
        mpfr_add(r25706, r25703, r25705, MPFR_RNDN);
        return mpfr_get_d(r25706, MPFR_RNDN);
}

static mpfr_t r25707, r25708, r25709, r25710, r25711;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25707, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r25708, "10", 10, MPFR_RNDN);
        mpfr_init(r25709);
        mpfr_init(r25710);
        mpfr_init(r25711);
}

double f_fm(double d) {
        ;
        ;
        mpfr_add(r25709, r25707, r25708, MPFR_RNDN);
        mpfr_set_d(r25710, d, MPFR_RNDN);
        mpfr_mul(r25711, r25709, r25710, MPFR_RNDN);
        return mpfr_get_d(r25711, MPFR_RNDN);
}

static mpfr_t r25712, r25713, r25714, r25715, r25716;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r25712, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r25713, "10", 10, MPFR_RNDN);
        mpfr_init(r25714);
        mpfr_init(r25715);
        mpfr_init(r25716);
}

double f_dm(double d) {
        ;
        ;
        mpfr_add(r25714, r25712, r25713, MPFR_RNDN);
        mpfr_set_d(r25715, d, MPFR_RNDN);
        mpfr_mul(r25716, r25714, r25715, MPFR_RNDN);
        return mpfr_get_d(r25716, MPFR_RNDN);
}