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

char *name = "2cbrt (problem 3.3.4)";

double f_if(float x) {
        float r9570 = x;
        float r9571 = 1;
        float r9572 = r9570 + r9571;
        float r9573 = cbrt(r9572);
        float r9574 = cbrt(r9570);
        float r9575 = r9573 - r9574;
        return r9575;
}

double f_id(double x) {
        double r9576 = x;
        double r9577 = 1;
        double r9578 = r9576 + r9577;
        double r9579 = cbrt(r9578);
        double r9580 = cbrt(r9576);
        double r9581 = r9579 - r9580;
        return r9581;
}


double f_of(float x) {
        float r9582 = 1;
        float r9583 = x;
        float r9584 = r9582 + r9583;
        float r9585 = cbrt(r9584);
        float r9586 = cbrt(r9583);
        float r9587 = r9585 + r9586;
        float r9588 = r9587 * r9586;
        float r9589 = r9585 * r9585;
        float r9590 = cbrt(r9589);
        float r9591 = cbrt(r9585);
        float r9592 = r9590 * r9591;
        float r9593 = cbrt(r9592);
        float r9594 = r9590 * r9593;
        float r9595 = r9585 * r9594;
        float r9596 = r9588 + r9595;
        float r9597 = r9582 / r9596;
        return r9597;
}

double f_od(double x) {
        double r9598 = 1;
        double r9599 = x;
        double r9600 = r9598 + r9599;
        double r9601 = cbrt(r9600);
        double r9602 = cbrt(r9599);
        double r9603 = r9601 + r9602;
        double r9604 = r9603 * r9602;
        double r9605 = r9601 * r9601;
        double r9606 = cbrt(r9605);
        double r9607 = cbrt(r9601);
        double r9608 = r9606 * r9607;
        double r9609 = cbrt(r9608);
        double r9610 = r9606 * r9609;
        double r9611 = r9601 * r9610;
        double r9612 = r9604 + r9611;
        double r9613 = r9598 / r9612;
        return r9613;
}

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 r9614, r9615, r9616, r9617, r9618, r9619;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9614);
        mpfr_init_set_str(r9615, "1", 10, MPFR_RNDN);
        mpfr_init(r9616);
        mpfr_init(r9617);
        mpfr_init(r9618);
        mpfr_init(r9619);
}

double f_im(double x) {
        mpfr_set_d(r9614, x, MPFR_RNDN);
        ;
        mpfr_add(r9616, r9614, r9615, MPFR_RNDN);
        mpfr_cbrt(r9617, r9616, MPFR_RNDN);
        mpfr_cbrt(r9618, r9614, MPFR_RNDN);
        mpfr_sub(r9619, r9617, r9618, MPFR_RNDN);
        return mpfr_get_d(r9619, MPFR_RNDN);
}

static mpfr_t r9620, r9621, r9622, r9623, r9624, r9625, r9626, r9627, r9628, r9629, r9630, r9631, r9632, r9633, r9634, r9635;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r9620, "1", 10, MPFR_RNDN);
        mpfr_init(r9621);
        mpfr_init(r9622);
        mpfr_init(r9623);
        mpfr_init(r9624);
        mpfr_init(r9625);
        mpfr_init(r9626);
        mpfr_init(r9627);
        mpfr_init(r9628);
        mpfr_init(r9629);
        mpfr_init(r9630);
        mpfr_init(r9631);
        mpfr_init(r9632);
        mpfr_init(r9633);
        mpfr_init(r9634);
        mpfr_init(r9635);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r9621, x, MPFR_RNDN);
        mpfr_add(r9622, r9620, r9621, MPFR_RNDN);
        mpfr_cbrt(r9623, r9622, MPFR_RNDN);
        mpfr_cbrt(r9624, r9621, MPFR_RNDN);
        mpfr_add(r9625, r9623, r9624, MPFR_RNDN);
        mpfr_mul(r9626, r9625, r9624, MPFR_RNDN);
        mpfr_mul(r9627, r9623, r9623, MPFR_RNDN);
        mpfr_cbrt(r9628, r9627, MPFR_RNDN);
        mpfr_cbrt(r9629, r9623, MPFR_RNDN);
        mpfr_mul(r9630, r9628, r9629, MPFR_RNDN);
        mpfr_cbrt(r9631, r9630, MPFR_RNDN);
        mpfr_mul(r9632, r9628, r9631, MPFR_RNDN);
        mpfr_mul(r9633, r9623, r9632, MPFR_RNDN);
        mpfr_add(r9634, r9626, r9633, MPFR_RNDN);
        mpfr_div(r9635, r9620, r9634, MPFR_RNDN);
        return mpfr_get_d(r9635, MPFR_RNDN);
}

static mpfr_t r9636, r9637, r9638, r9639, r9640, r9641, r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r9636, "1", 10, MPFR_RNDN);
        mpfr_init(r9637);
        mpfr_init(r9638);
        mpfr_init(r9639);
        mpfr_init(r9640);
        mpfr_init(r9641);
        mpfr_init(r9642);
        mpfr_init(r9643);
        mpfr_init(r9644);
        mpfr_init(r9645);
        mpfr_init(r9646);
        mpfr_init(r9647);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init(r9651);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r9637, x, MPFR_RNDN);
        mpfr_add(r9638, r9636, r9637, MPFR_RNDN);
        mpfr_cbrt(r9639, r9638, MPFR_RNDN);
        mpfr_cbrt(r9640, r9637, MPFR_RNDN);
        mpfr_add(r9641, r9639, r9640, MPFR_RNDN);
        mpfr_mul(r9642, r9641, r9640, MPFR_RNDN);
        mpfr_mul(r9643, r9639, r9639, MPFR_RNDN);
        mpfr_cbrt(r9644, r9643, MPFR_RNDN);
        mpfr_cbrt(r9645, r9639, MPFR_RNDN);
        mpfr_mul(r9646, r9644, r9645, MPFR_RNDN);
        mpfr_cbrt(r9647, r9646, MPFR_RNDN);
        mpfr_mul(r9648, r9644, r9647, MPFR_RNDN);
        mpfr_mul(r9649, r9639, r9648, MPFR_RNDN);
        mpfr_add(r9650, r9642, r9649, MPFR_RNDN);
        mpfr_div(r9651, r9636, r9650, MPFR_RNDN);
        return mpfr_get_d(r9651, MPFR_RNDN);
}