#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 r25058 = d1;
        float r25059 = r25058 * r25058;
        float r25060 = r25059 * r25058;
        float r25061 = r25060 * r25058;
        return r25061;
}

double f_id(double d1) {
        double r25062 = d1;
        double r25063 = r25062 * r25062;
        double r25064 = r25063 * r25062;
        double r25065 = r25064 * r25062;
        return r25065;
}


double f_of(float d1) {
        float r25066 = d1;
        float r25067 = 3;
        float r25068 = 1;
        float r25069 = r25067 + r25068;
        float r25070 = pow(r25066, r25069);
        return r25070;
}

double f_od(double d1) {
        double r25071 = d1;
        double r25072 = 3;
        double r25073 = 1;
        double r25074 = r25072 + r25073;
        double r25075 = pow(r25071, r25074);
        return r25075;
}

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 r25076, r25077, r25078, r25079;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r25076);
        mpfr_init(r25077);
        mpfr_init(r25078);
        mpfr_init(r25079);
}

double f_im(double d1) {
        mpfr_set_d(r25076, d1, MPFR_RNDN);
        mpfr_mul(r25077, r25076, r25076, MPFR_RNDN);
        mpfr_mul(r25078, r25077, r25076, MPFR_RNDN);
        mpfr_mul(r25079, r25078, r25076, MPFR_RNDN);
        return mpfr_get_d(r25079, MPFR_RNDN);
}

static mpfr_t r25080, r25081, r25082, r25083, r25084;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25080);
        mpfr_init_set_str(r25081, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r25082, "1", 10, MPFR_RNDN);
        mpfr_init(r25083);
        mpfr_init(r25084);
}

double f_fm(double d1) {
        mpfr_set_d(r25080, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r25083, r25081, r25082, MPFR_RNDN);
        mpfr_pow(r25084, r25080, r25083, MPFR_RNDN);
        return mpfr_get_d(r25084, MPFR_RNDN);
}

static mpfr_t r25085, r25086, r25087, r25088, r25089;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25085);
        mpfr_init_set_str(r25086, "3", 10, MPFR_RNDN);
        mpfr_init_set_str(r25087, "1", 10, MPFR_RNDN);
        mpfr_init(r25088);
        mpfr_init(r25089);
}

double f_dm(double d1) {
        mpfr_set_d(r25085, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r25088, r25086, r25087, MPFR_RNDN);
        mpfr_pow(r25089, r25085, r25088, MPFR_RNDN);
        return mpfr_get_d(r25089, MPFR_RNDN);
}