#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 r50499 = d;
        float r50500 = 10;
        float r50501 = r50499 * r50500;
        float r50502 = 20;
        float r50503 = r50499 * r50502;
        float r50504 = r50501 + r50503;
        return r50504;
}

double f_id(double d) {
        double r50505 = d;
        double r50506 = 10;
        double r50507 = r50505 * r50506;
        double r50508 = 20;
        double r50509 = r50505 * r50508;
        double r50510 = r50507 + r50509;
        return r50510;
}


double f_of(float d) {
        float r50511 = d;
        float r50512 = 20;
        float r50513 = 10;
        float r50514 = r50512 + r50513;
        float r50515 = r50511 * r50514;
        return r50515;
}

double f_od(double d) {
        double r50516 = d;
        double r50517 = 20;
        double r50518 = 10;
        double r50519 = r50517 + r50518;
        double r50520 = r50516 * r50519;
        return r50520;
}

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 r50521, r50522, r50523, r50524, r50525, r50526;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r50521);
        mpfr_init_set_str(r50522, "10", 10, MPFR_RNDN);
        mpfr_init(r50523);
        mpfr_init_set_str(r50524, "20", 10, MPFR_RNDN);
        mpfr_init(r50525);
        mpfr_init(r50526);
}

double f_im(double d) {
        mpfr_set_d(r50521, d, MPFR_RNDN);
        ;
        mpfr_mul(r50523, r50521, r50522, MPFR_RNDN);
        ;
        mpfr_mul(r50525, r50521, r50524, MPFR_RNDN);
        mpfr_add(r50526, r50523, r50525, MPFR_RNDN);
        return mpfr_get_d(r50526, MPFR_RNDN);
}

static mpfr_t r50527, r50528, r50529, r50530, r50531;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50527);
        mpfr_init_set_str(r50528, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r50529, "10", 10, MPFR_RNDN);
        mpfr_init(r50530);
        mpfr_init(r50531);
}

double f_fm(double d) {
        mpfr_set_d(r50527, d, MPFR_RNDN);
        ;
        ;
        mpfr_add(r50530, r50528, r50529, MPFR_RNDN);
        mpfr_mul(r50531, r50527, r50530, MPFR_RNDN);
        return mpfr_get_d(r50531, MPFR_RNDN);
}

static mpfr_t r50532, r50533, r50534, r50535, r50536;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r50532);
        mpfr_init_set_str(r50533, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r50534, "10", 10, MPFR_RNDN);
        mpfr_init(r50535);
        mpfr_init(r50536);
}

double f_dm(double d) {
        mpfr_set_d(r50532, d, MPFR_RNDN);
        ;
        ;
        mpfr_add(r50535, r50533, r50534, MPFR_RNDN);
        mpfr_mul(r50536, r50532, r50535, MPFR_RNDN);
        return mpfr_get_d(r50536, MPFR_RNDN);
}