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

char *name = "math.cube on complex, imaginary part";

double f_if(float x_re, float x_im) {
        float r11364 = x_re;
        float r11365 = r11364 * r11364;
        float r11366 = x_im;
        float r11367 = r11366 * r11366;
        float r11368 = r11365 - r11367;
        float r11369 = r11368 * r11366;
        float r11370 = r11364 * r11366;
        float r11371 = r11366 * r11364;
        float r11372 = r11370 + r11371;
        float r11373 = r11372 * r11364;
        float r11374 = r11369 + r11373;
        return r11374;
}

double f_id(double x_re, double x_im) {
        double r11375 = x_re;
        double r11376 = r11375 * r11375;
        double r11377 = x_im;
        double r11378 = r11377 * r11377;
        double r11379 = r11376 - r11378;
        double r11380 = r11379 * r11377;
        double r11381 = r11375 * r11377;
        double r11382 = r11377 * r11375;
        double r11383 = r11381 + r11382;
        double r11384 = r11383 * r11375;
        double r11385 = r11380 + r11384;
        return r11385;
}


double f_of(float x_re, float x_im) {
        float r11386 = x_re;
        float r11387 = 3.0f;
        float r11388 = x_im;
        float r11389 = r11387 * r11388;
        float r11390 = r11389 * r11386;
        float r11391 = r11386 * r11390;
        float r11392 = pow(r11388, r11387);
        float r11393 = r11391 - r11392;
        return r11393;
}

double f_od(double x_re, double x_im) {
        double r11394 = x_re;
        double r11395 = 3.0;
        double r11396 = x_im;
        double r11397 = r11395 * r11396;
        double r11398 = r11397 * r11394;
        double r11399 = r11394 * r11398;
        double r11400 = pow(r11396, r11395);
        double r11401 = r11399 - r11400;
        return r11401;
}

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 r11402, r11403, r11404, r11405, r11406, r11407, r11408, r11409, r11410, r11411, r11412;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r11402);
        mpfr_init(r11403);
        mpfr_init(r11404);
        mpfr_init(r11405);
        mpfr_init(r11406);
        mpfr_init(r11407);
        mpfr_init(r11408);
        mpfr_init(r11409);
        mpfr_init(r11410);
        mpfr_init(r11411);
        mpfr_init(r11412);
}

double f_im(double x_re, double x_im) {
        mpfr_set_d(r11402, x_re, MPFR_RNDN);
        mpfr_mul(r11403, r11402, r11402, MPFR_RNDN);
        mpfr_set_d(r11404, x_im, MPFR_RNDN);
        mpfr_mul(r11405, r11404, r11404, MPFR_RNDN);
        mpfr_sub(r11406, r11403, r11405, MPFR_RNDN);
        mpfr_mul(r11407, r11406, r11404, MPFR_RNDN);
        mpfr_mul(r11408, r11402, r11404, MPFR_RNDN);
        mpfr_mul(r11409, r11404, r11402, MPFR_RNDN);
        mpfr_add(r11410, r11408, r11409, MPFR_RNDN);
        mpfr_mul(r11411, r11410, r11402, MPFR_RNDN);
        mpfr_add(r11412, r11407, r11411, MPFR_RNDN);
        return mpfr_get_d(r11412, MPFR_RNDN);
}

static mpfr_t r11413, r11414, r11415, r11416, r11417, r11418, r11419, r11420;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r11413);
        mpfr_init_set_str(r11414, "3", 10, MPFR_RNDN);
        mpfr_init(r11415);
        mpfr_init(r11416);
        mpfr_init(r11417);
        mpfr_init(r11418);
        mpfr_init(r11419);
        mpfr_init(r11420);
}

double f_fm(double x_re, double x_im) {
        mpfr_set_d(r11413, x_re, MPFR_RNDN);
        ;
        mpfr_set_d(r11415, x_im, MPFR_RNDN);
        mpfr_mul(r11416, r11414, r11415, MPFR_RNDN);
        mpfr_mul(r11417, r11416, r11413, MPFR_RNDN);
        mpfr_mul(r11418, r11413, r11417, MPFR_RNDN);
        mpfr_pow(r11419, r11415, r11414, MPFR_RNDN);
        mpfr_sub(r11420, r11418, r11419, MPFR_RNDN);
        return mpfr_get_d(r11420, MPFR_RNDN);
}

static mpfr_t r11421, r11422, r11423, r11424, r11425, r11426, r11427, r11428;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r11421);
        mpfr_init_set_str(r11422, "3", 10, MPFR_RNDN);
        mpfr_init(r11423);
        mpfr_init(r11424);
        mpfr_init(r11425);
        mpfr_init(r11426);
        mpfr_init(r11427);
        mpfr_init(r11428);
}

double f_dm(double x_re, double x_im) {
        mpfr_set_d(r11421, x_re, MPFR_RNDN);
        ;
        mpfr_set_d(r11423, x_im, MPFR_RNDN);
        mpfr_mul(r11424, r11422, r11423, MPFR_RNDN);
        mpfr_mul(r11425, r11424, r11421, MPFR_RNDN);
        mpfr_mul(r11426, r11421, r11425, MPFR_RNDN);
        mpfr_pow(r11427, r11423, r11422, MPFR_RNDN);
        mpfr_sub(r11428, r11426, r11427, MPFR_RNDN);
        return mpfr_get_d(r11428, MPFR_RNDN);
}