#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 r26170 = d1;
        float r26171 = r26170 * r26170;
        float r26172 = r26171 * r26170;
        float r26173 = r26172 * r26170;
        return r26173;
}

double f_id(double d1) {
        double r26174 = d1;
        double r26175 = r26174 * r26174;
        double r26176 = r26175 * r26174;
        double r26177 = r26176 * r26174;
        return r26177;
}


double f_of(float d1) {
        float r26178 = d1;
        float r26179 = 3;
        float r26180 = 1;
        float r26181 = r26179 + r26180;
        float r26182 = pow(r26178, r26181);
        return r26182;
}

double f_od(double d1) {
        double r26183 = d1;
        double r26184 = 3;
        double r26185 = 1;
        double r26186 = r26184 + r26185;
        double r26187 = pow(r26183, r26186);
        return r26187;
}

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 r26188, r26189, r26190, r26191;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26188);
        mpfr_init(r26189);
        mpfr_init(r26190);
        mpfr_init(r26191);
}

double f_im(double d1) {
        mpfr_set_d(r26188, d1, MPFR_RNDN);
        mpfr_mul(r26189, r26188, r26188, MPFR_RNDN);
        mpfr_mul(r26190, r26189, r26188, MPFR_RNDN);
        mpfr_mul(r26191, r26190, r26188, MPFR_RNDN);
        return mpfr_get_d(r26191, MPFR_RNDN);
}

static mpfr_t r26192, r26193, r26194, r26195, r26196;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26192);
        mpfr_init_set_str(r26193, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26194, "1", 10, MPFR_RNDN);
        mpfr_init(r26195);
        mpfr_init(r26196);
}

double f_fm(double d1) {
        mpfr_set_d(r26192, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26195, r26193, r26194, MPFR_RNDN);
        mpfr_pow(r26196, r26192, r26195, MPFR_RNDN);
        return mpfr_get_d(r26196, MPFR_RNDN);
}

static mpfr_t r26197, r26198, r26199, r26200, r26201;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26197);
        mpfr_init_set_str(r26198, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26199, "1", 10, MPFR_RNDN);
        mpfr_init(r26200);
        mpfr_init(r26201);
}

double f_dm(double d1) {
        mpfr_set_d(r26197, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26200, r26198, r26199, MPFR_RNDN);
        mpfr_pow(r26201, r26197, r26200, MPFR_RNDN);
        return mpfr_get_d(r26201, MPFR_RNDN);
}