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

char *name = "FastMath repmul";

double f_if(float d1) {
        float r26186 = d1;
        float r26187 = r26186 * r26186;
        float r26188 = r26187 * r26186;
        float r26189 = r26188 * r26186;
        return r26189;
}

double f_id(double d1) {
        double r26190 = d1;
        double r26191 = r26190 * r26190;
        double r26192 = r26191 * r26190;
        double r26193 = r26192 * r26190;
        return r26193;
}


double f_of(float d1) {
        float r26194 = d1;
        float r26195 = 3;
        float r26196 = 1;
        float r26197 = r26195 + r26196;
        float r26198 = pow(r26194, r26197);
        return r26198;
}

double f_od(double d1) {
        double r26199 = d1;
        double r26200 = 3;
        double r26201 = 1;
        double r26202 = r26200 + r26201;
        double r26203 = pow(r26199, r26202);
        return r26203;
}

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 r26204, r26205, r26206, r26207;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26204);
        mpfr_init(r26205);
        mpfr_init(r26206);
        mpfr_init(r26207);
}

double f_im(double d1) {
        mpfr_set_d(r26204, d1, MPFR_RNDN);
        mpfr_mul(r26205, r26204, r26204, MPFR_RNDN);
        mpfr_mul(r26206, r26205, r26204, MPFR_RNDN);
        mpfr_mul(r26207, r26206, r26204, MPFR_RNDN);
        return mpfr_get_d(r26207, MPFR_RNDN);
}

static mpfr_t r26208, r26209, r26210, r26211, r26212;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26208);
        mpfr_init_set_str(r26209, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26210, "1", 10, MPFR_RNDN);
        mpfr_init(r26211);
        mpfr_init(r26212);
}

double f_fm(double d1) {
        mpfr_set_d(r26208, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26211, r26209, r26210, MPFR_RNDN);
        mpfr_pow(r26212, r26208, r26211, MPFR_RNDN);
        return mpfr_get_d(r26212, MPFR_RNDN);
}

static mpfr_t r26213, r26214, r26215, r26216, r26217;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26213);
        mpfr_init_set_str(r26214, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26215, "1", 10, MPFR_RNDN);
        mpfr_init(r26216);
        mpfr_init(r26217);
}

double f_dm(double d1) {
        mpfr_set_d(r26213, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26216, r26214, r26215, MPFR_RNDN);
        mpfr_pow(r26217, r26213, r26216, MPFR_RNDN);
        return mpfr_get_d(r26217, MPFR_RNDN);
}