#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 r26765 = d1;
        float r26766 = r26765 * r26765;
        float r26767 = r26766 * r26765;
        float r26768 = r26767 * r26765;
        return r26768;
}

double f_id(double d1) {
        double r26769 = d1;
        double r26770 = r26769 * r26769;
        double r26771 = r26770 * r26769;
        double r26772 = r26771 * r26769;
        return r26772;
}


double f_of(float d1) {
        float r26773 = d1;
        float r26774 = 3;
        float r26775 = 1;
        float r26776 = r26774 + r26775;
        float r26777 = pow(r26773, r26776);
        return r26777;
}

double f_od(double d1) {
        double r26778 = d1;
        double r26779 = 3;
        double r26780 = 1;
        double r26781 = r26779 + r26780;
        double r26782 = pow(r26778, r26781);
        return r26782;
}

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 r26783, r26784, r26785, r26786;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26783);
        mpfr_init(r26784);
        mpfr_init(r26785);
        mpfr_init(r26786);
}

double f_im(double d1) {
        mpfr_set_d(r26783, d1, MPFR_RNDN);
        mpfr_mul(r26784, r26783, r26783, MPFR_RNDN);
        mpfr_mul(r26785, r26784, r26783, MPFR_RNDN);
        mpfr_mul(r26786, r26785, r26783, MPFR_RNDN);
        return mpfr_get_d(r26786, MPFR_RNDN);
}

static mpfr_t r26787, r26788, r26789, r26790, r26791;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26787);
        mpfr_init_set_str(r26788, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26789, "1", 10, MPFR_RNDN);
        mpfr_init(r26790);
        mpfr_init(r26791);
}

double f_fm(double d1) {
        mpfr_set_d(r26787, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26790, r26788, r26789, MPFR_RNDN);
        mpfr_pow(r26791, r26787, r26790, MPFR_RNDN);
        return mpfr_get_d(r26791, MPFR_RNDN);
}

static mpfr_t r26792, r26793, r26794, r26795, r26796;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26792);
        mpfr_init_set_str(r26793, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26794, "1", 10, MPFR_RNDN);
        mpfr_init(r26795);
        mpfr_init(r26796);
}

double f_dm(double d1) {
        mpfr_set_d(r26792, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r26795, r26793, r26794, MPFR_RNDN);
        mpfr_pow(r26796, r26792, r26795, MPFR_RNDN);
        return mpfr_get_d(r26796, MPFR_RNDN);
}