#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 r26001 = x;
        float r26002 = 1;
        float r26003 = r26001 + r26002;
        float r26004 = cbrt(r26003);
        float r26005 = cbrt(r26001);
        float r26006 = r26004 - r26005;
        return r26006;
}

double f_id(double x) {
        double r26007 = x;
        double r26008 = 1;
        double r26009 = r26007 + r26008;
        double r26010 = cbrt(r26009);
        double r26011 = cbrt(r26007);
        double r26012 = r26010 - r26011;
        return r26012;
}


double f_of(float x) {
        float r26013 = 1;
        float r26014 = x;
        float r26015 = r26013 + r26014;
        float r26016 = cbrt(r26015);
        float r26017 = r26016 * r26016;
        float r26018 = cbrt(r26014);
        float r26019 = r26018 + r26016;
        float r26020 = cbrt(r26018);
        float r26021 = r26020 * r26020;
        float r26022 = r26019 * r26021;
        float r26023 = r26022 * r26020;
        float r26024 = r26017 + r26023;
        float r26025 = r26013 / r26024;
        return r26025;
}

double f_od(double x) {
        double r26026 = 1;
        double r26027 = x;
        double r26028 = r26026 + r26027;
        double r26029 = cbrt(r26028);
        double r26030 = r26029 * r26029;
        double r26031 = cbrt(r26027);
        double r26032 = r26031 + r26029;
        double r26033 = cbrt(r26031);
        double r26034 = r26033 * r26033;
        double r26035 = r26032 * r26034;
        double r26036 = r26035 * r26033;
        double r26037 = r26030 + r26036;
        double r26038 = r26026 / r26037;
        return r26038;
}

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 r26039, r26040, r26041, r26042, r26043, r26044;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26039);
        mpfr_init_set_str(r26040, "1", 10, MPFR_RNDN);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init(r26044);
}

double f_im(double x) {
        mpfr_set_d(r26039, x, MPFR_RNDN);
        ;
        mpfr_add(r26041, r26039, r26040, MPFR_RNDN);
        mpfr_cbrt(r26042, r26041, MPFR_RNDN);
        mpfr_cbrt(r26043, r26039, MPFR_RNDN);
        mpfr_sub(r26044, r26042, r26043, MPFR_RNDN);
        return mpfr_get_d(r26044, MPFR_RNDN);
}

static mpfr_t r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26045, "1", 10, MPFR_RNDN);
        mpfr_init(r26046);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
        mpfr_init(r26056);
        mpfr_init(r26057);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r26046, x, MPFR_RNDN);
        mpfr_add(r26047, r26045, r26046, MPFR_RNDN);
        mpfr_cbrt(r26048, r26047, MPFR_RNDN);
        mpfr_mul(r26049, r26048, r26048, MPFR_RNDN);
        mpfr_cbrt(r26050, r26046, MPFR_RNDN);
        mpfr_add(r26051, r26050, r26048, MPFR_RNDN);
        mpfr_cbrt(r26052, r26050, MPFR_RNDN);
        mpfr_mul(r26053, r26052, r26052, MPFR_RNDN);
        mpfr_mul(r26054, r26051, r26053, MPFR_RNDN);
        mpfr_mul(r26055, r26054, r26052, MPFR_RNDN);
        mpfr_add(r26056, r26049, r26055, MPFR_RNDN);
        mpfr_div(r26057, r26045, r26056, MPFR_RNDN);
        return mpfr_get_d(r26057, MPFR_RNDN);
}

static mpfr_t r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26058, "1", 10, MPFR_RNDN);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init(r26061);
        mpfr_init(r26062);
        mpfr_init(r26063);
        mpfr_init(r26064);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init(r26068);
        mpfr_init(r26069);
        mpfr_init(r26070);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r26059, x, MPFR_RNDN);
        mpfr_add(r26060, r26058, r26059, MPFR_RNDN);
        mpfr_cbrt(r26061, r26060, MPFR_RNDN);
        mpfr_mul(r26062, r26061, r26061, MPFR_RNDN);
        mpfr_cbrt(r26063, r26059, MPFR_RNDN);
        mpfr_add(r26064, r26063, r26061, MPFR_RNDN);
        mpfr_cbrt(r26065, r26063, MPFR_RNDN);
        mpfr_mul(r26066, r26065, r26065, MPFR_RNDN);
        mpfr_mul(r26067, r26064, r26066, MPFR_RNDN);
        mpfr_mul(r26068, r26067, r26065, MPFR_RNDN);
        mpfr_add(r26069, r26062, r26068, MPFR_RNDN);
        mpfr_div(r26070, r26058, r26069, MPFR_RNDN);
        return mpfr_get_d(r26070, MPFR_RNDN);
}