#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 r18475 = d1;
        float r18476 = r18475 * r18475;
        float r18477 = r18475 * r18476;
        float r18478 = r18477 * r18475;
        float r18479 = r18478 * r18475;
        float r18480 = r18479 * r18476;
        float r18481 = r18480 * r18475;
        float r18482 = r18475 * r18481;
        float r18483 = r18482 * r18475;
        return r18483;
}

double f_id(double d1) {
        double r18484 = d1;
        double r18485 = r18484 * r18484;
        double r18486 = r18484 * r18485;
        double r18487 = r18486 * r18484;
        double r18488 = r18487 * r18484;
        double r18489 = r18488 * r18485;
        double r18490 = r18489 * r18484;
        double r18491 = r18484 * r18490;
        double r18492 = r18491 * r18484;
        return r18492;
}


double f_of(float d1) {
        float r18493 = d1;
        float r18494 = 2.0f;
        float r18495 = 3.0f;
        float r18496 = r18494 + r18495;
        float r18497 = r18496 + r18496;
        float r18498 = pow(r18493, r18497);
        return r18498;
}

double f_od(double d1) {
        double r18499 = d1;
        double r18500 = 2.0;
        double r18501 = 3.0;
        double r18502 = r18500 + r18501;
        double r18503 = r18502 + r18502;
        double r18504 = pow(r18499, r18503);
        return r18504;
}

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 r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18505);
        mpfr_init(r18506);
        mpfr_init(r18507);
        mpfr_init(r18508);
        mpfr_init(r18509);
        mpfr_init(r18510);
        mpfr_init(r18511);
        mpfr_init(r18512);
        mpfr_init(r18513);
}

double f_im(double d1) {
        mpfr_set_d(r18505, d1, MPFR_RNDN);
        mpfr_mul(r18506, r18505, r18505, MPFR_RNDN);
        mpfr_mul(r18507, r18505, r18506, MPFR_RNDN);
        mpfr_mul(r18508, r18507, r18505, MPFR_RNDN);
        mpfr_mul(r18509, r18508, r18505, MPFR_RNDN);
        mpfr_mul(r18510, r18509, r18506, MPFR_RNDN);
        mpfr_mul(r18511, r18510, r18505, MPFR_RNDN);
        mpfr_mul(r18512, r18505, r18511, MPFR_RNDN);
        mpfr_mul(r18513, r18512, r18505, MPFR_RNDN);
        return mpfr_get_d(r18513, MPFR_RNDN);
}

static mpfr_t r18514, r18515, r18516, r18517, r18518, r18519;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18514);
        mpfr_init_set_str(r18515, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r18516, "3", 10, MPFR_RNDN);
        mpfr_init(r18517);
        mpfr_init(r18518);
        mpfr_init(r18519);
}

double f_fm(double d1) {
        mpfr_set_d(r18514, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r18517, r18515, r18516, MPFR_RNDN);
        mpfr_add(r18518, r18517, r18517, MPFR_RNDN);
        mpfr_pow(r18519, r18514, r18518, MPFR_RNDN);
        return mpfr_get_d(r18519, MPFR_RNDN);
}

static mpfr_t r18520, r18521, r18522, r18523, r18524, r18525;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18520);
        mpfr_init_set_str(r18521, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r18522, "3", 10, MPFR_RNDN);
        mpfr_init(r18523);
        mpfr_init(r18524);
        mpfr_init(r18525);
}

double f_dm(double d1) {
        mpfr_set_d(r18520, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r18523, r18521, r18522, MPFR_RNDN);
        mpfr_add(r18524, r18523, r18523, MPFR_RNDN);
        mpfr_pow(r18525, r18520, r18524, MPFR_RNDN);
        return mpfr_get_d(r18525, MPFR_RNDN);
}

