#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 r9434 = x_re;
        float r9435 = r9434 * r9434;
        float r9436 = x_im;
        float r9437 = r9436 * r9436;
        float r9438 = r9435 - r9437;
        float r9439 = r9438 * r9436;
        float r9440 = r9434 * r9436;
        float r9441 = r9436 * r9434;
        float r9442 = r9440 + r9441;
        float r9443 = r9442 * r9434;
        float r9444 = r9439 + r9443;
        return r9444;
}

double f_id(double x_re, double x_im) {
        double r9445 = x_re;
        double r9446 = r9445 * r9445;
        double r9447 = x_im;
        double r9448 = r9447 * r9447;
        double r9449 = r9446 - r9448;
        double r9450 = r9449 * r9447;
        double r9451 = r9445 * r9447;
        double r9452 = r9447 * r9445;
        double r9453 = r9451 + r9452;
        double r9454 = r9453 * r9445;
        double r9455 = r9450 + r9454;
        return r9455;
}


double f_of(float x_re, float x_im) {
        float r9456 = x_re;
        float r9457 = x_im;
        float r9458 = r9456 - r9457;
        float r9459 = r9457 + r9456;
        float r9460 = r9459 * r9457;
        float r9461 = r9457 + r9457;
        float r9462 = r9461 * r9456;
        float r9463 = r9462 * r9456;
        float r9464 = fma(r9458, r9460, r9463);
        return r9464;
}

double f_od(double x_re, double x_im) {
        double r9465 = x_re;
        double r9466 = x_im;
        double r9467 = r9465 - r9466;
        double r9468 = r9466 + r9465;
        double r9469 = r9468 * r9466;
        double r9470 = r9466 + r9466;
        double r9471 = r9470 * r9465;
        double r9472 = r9471 * r9465;
        double r9473 = fma(r9467, r9469, r9472);
        return r9473;
}

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 r9474, r9475, r9476, r9477, r9478, r9479, r9480, r9481, r9482, r9483, r9484;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9474);
        mpfr_init(r9475);
        mpfr_init(r9476);
        mpfr_init(r9477);
        mpfr_init(r9478);
        mpfr_init(r9479);
        mpfr_init(r9480);
        mpfr_init(r9481);
        mpfr_init(r9482);
        mpfr_init(r9483);
        mpfr_init(r9484);
}

double f_im(double x_re, double x_im) {
        mpfr_set_d(r9474, x_re, MPFR_RNDN);
        mpfr_mul(r9475, r9474, r9474, MPFR_RNDN);
        mpfr_set_d(r9476, x_im, MPFR_RNDN);
        mpfr_mul(r9477, r9476, r9476, MPFR_RNDN);
        mpfr_sub(r9478, r9475, r9477, MPFR_RNDN);
        mpfr_mul(r9479, r9478, r9476, MPFR_RNDN);
        mpfr_mul(r9480, r9474, r9476, MPFR_RNDN);
        mpfr_mul(r9481, r9476, r9474, MPFR_RNDN);
        mpfr_add(r9482, r9480, r9481, MPFR_RNDN);
        mpfr_mul(r9483, r9482, r9474, MPFR_RNDN);
        mpfr_add(r9484, r9479, r9483, MPFR_RNDN);
        return mpfr_get_d(r9484, MPFR_RNDN);
}

static mpfr_t r9485, r9486, r9487, r9488, r9489, r9490, r9491, r9492, r9493;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9485);
        mpfr_init(r9486);
        mpfr_init(r9487);
        mpfr_init(r9488);
        mpfr_init(r9489);
        mpfr_init(r9490);
        mpfr_init(r9491);
        mpfr_init(r9492);
        mpfr_init(r9493);
}

double f_fm(double x_re, double x_im) {
        mpfr_set_d(r9485, x_re, MPFR_RNDN);
        mpfr_set_d(r9486, x_im, MPFR_RNDN);
        mpfr_sub(r9487, r9485, r9486, MPFR_RNDN);
        mpfr_add(r9488, r9486, r9485, MPFR_RNDN);
        mpfr_mul(r9489, r9488, r9486, MPFR_RNDN);
        mpfr_add(r9490, r9486, r9486, MPFR_RNDN);
        mpfr_mul(r9491, r9490, r9485, MPFR_RNDN);
        mpfr_mul(r9492, r9491, r9485, MPFR_RNDN);
        mpfr_fma(r9493, r9487, r9489, r9492, MPFR_RNDN);
        return mpfr_get_d(r9493, MPFR_RNDN);
}

static mpfr_t r9494, r9495, r9496, r9497, r9498, r9499, r9500, r9501, r9502;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9494);
        mpfr_init(r9495);
        mpfr_init(r9496);
        mpfr_init(r9497);
        mpfr_init(r9498);
        mpfr_init(r9499);
        mpfr_init(r9500);
        mpfr_init(r9501);
        mpfr_init(r9502);
}

double f_dm(double x_re, double x_im) {
        mpfr_set_d(r9494, x_re, MPFR_RNDN);
        mpfr_set_d(r9495, x_im, MPFR_RNDN);
        mpfr_sub(r9496, r9494, r9495, MPFR_RNDN);
        mpfr_add(r9497, r9495, r9494, MPFR_RNDN);
        mpfr_mul(r9498, r9497, r9495, MPFR_RNDN);
        mpfr_add(r9499, r9495, r9495, MPFR_RNDN);
        mpfr_mul(r9500, r9499, r9494, MPFR_RNDN);
        mpfr_mul(r9501, r9500, r9494, MPFR_RNDN);
        mpfr_fma(r9502, r9496, r9498, r9501, MPFR_RNDN);
        return mpfr_get_d(r9502, MPFR_RNDN);
}