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

char *name = "FastMath test5";

double f_if(float d1) {
        float r20056 = d1;
        float r20057 = r20056 * r20056;
        float r20058 = r20056 * r20057;
        float r20059 = r20058 * r20056;
        float r20060 = r20059 * r20056;
        float r20061 = r20060 * r20057;
        float r20062 = r20061 * r20056;
        float r20063 = r20056 * r20062;
        float r20064 = r20063 * r20056;
        return r20064;
}

double f_id(double d1) {
        double r20065 = d1;
        double r20066 = r20065 * r20065;
        double r20067 = r20065 * r20066;
        double r20068 = r20067 * r20065;
        double r20069 = r20068 * r20065;
        double r20070 = r20069 * r20066;
        double r20071 = r20070 * r20065;
        double r20072 = r20065 * r20071;
        double r20073 = r20072 * r20065;
        return r20073;
}


double f_of(float d1) {
        float r20074 = d1;
        float r20075 = 2.0f;
        float r20076 = 3.0f;
        float r20077 = r20075 + r20076;
        float r20078 = r20077 + r20077;
        float r20079 = pow(r20074, r20078);
        return r20079;
}

double f_od(double d1) {
        double r20080 = d1;
        double r20081 = 2.0;
        double r20082 = 3.0;
        double r20083 = r20081 + r20082;
        double r20084 = r20083 + r20083;
        double r20085 = pow(r20080, r20084);
        return r20085;
}

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 r20086, r20087, r20088, r20089, r20090, r20091, r20092, r20093, r20094;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20086);
        mpfr_init(r20087);
        mpfr_init(r20088);
        mpfr_init(r20089);
        mpfr_init(r20090);
        mpfr_init(r20091);
        mpfr_init(r20092);
        mpfr_init(r20093);
        mpfr_init(r20094);
}

double f_im(double d1) {
        mpfr_set_d(r20086, d1, MPFR_RNDN);
        mpfr_mul(r20087, r20086, r20086, MPFR_RNDN);
        mpfr_mul(r20088, r20086, r20087, MPFR_RNDN);
        mpfr_mul(r20089, r20088, r20086, MPFR_RNDN);
        mpfr_mul(r20090, r20089, r20086, MPFR_RNDN);
        mpfr_mul(r20091, r20090, r20087, MPFR_RNDN);
        mpfr_mul(r20092, r20091, r20086, MPFR_RNDN);
        mpfr_mul(r20093, r20086, r20092, MPFR_RNDN);
        mpfr_mul(r20094, r20093, r20086, MPFR_RNDN);
        return mpfr_get_d(r20094, MPFR_RNDN);
}

static mpfr_t r20095, r20096, r20097, r20098, r20099, r20100;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20095);
        mpfr_init_set_str(r20096, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r20097, "3", 10, MPFR_RNDN);
        mpfr_init(r20098);
        mpfr_init(r20099);
        mpfr_init(r20100);
}

double f_fm(double d1) {
        mpfr_set_d(r20095, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r20098, r20096, r20097, MPFR_RNDN);
        mpfr_add(r20099, r20098, r20098, MPFR_RNDN);
        mpfr_pow(r20100, r20095, r20099, MPFR_RNDN);
        return mpfr_get_d(r20100, MPFR_RNDN);
}

static mpfr_t r20101, r20102, r20103, r20104, r20105, r20106;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20101);
        mpfr_init_set_str(r20102, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r20103, "3", 10, MPFR_RNDN);
        mpfr_init(r20104);
        mpfr_init(r20105);
        mpfr_init(r20106);
}

double f_dm(double d1) {
        mpfr_set_d(r20101, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r20104, r20102, r20103, MPFR_RNDN);
        mpfr_add(r20105, r20104, r20104, MPFR_RNDN);
        mpfr_pow(r20106, r20101, r20105, MPFR_RNDN);
        return mpfr_get_d(r20106, MPFR_RNDN);
}