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

char *name = "_divideComplex, real part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r9489 = x_re;
        float r9490 = y_re;
        float r9491 = r9489 * r9490;
        float r9492 = x_im;
        float r9493 = y_im;
        float r9494 = r9492 * r9493;
        float r9495 = r9491 + r9494;
        float r9496 = r9490 * r9490;
        float r9497 = r9493 * r9493;
        float r9498 = r9496 + r9497;
        float r9499 = r9495 / r9498;
        return r9499;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r9500 = x_re;
        double r9501 = y_re;
        double r9502 = r9500 * r9501;
        double r9503 = x_im;
        double r9504 = y_im;
        double r9505 = r9503 * r9504;
        double r9506 = r9502 + r9505;
        double r9507 = r9501 * r9501;
        double r9508 = r9504 * r9504;
        double r9509 = r9507 + r9508;
        double r9510 = r9506 / r9509;
        return r9510;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r9511 = x_im;
        float r9512 = y_im;
        float r9513 = r9511 * r9512;
        float r9514 = x_re;
        float r9515 = y_re;
        float r9516 = r9514 * r9515;
        float r9517 = r9513 + r9516;
        float r9518 = r9515 * r9515;
        float r9519 = r9512 * r9512;
        float r9520 = r9518 + r9519;
        float r9521 = r9517 / r9520;
        return r9521;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r9522 = x_im;
        double r9523 = y_im;
        double r9524 = r9522 * r9523;
        double r9525 = x_re;
        double r9526 = y_re;
        double r9527 = r9525 * r9526;
        double r9528 = r9524 + r9527;
        double r9529 = r9526 * r9526;
        double r9530 = r9523 * r9523;
        double r9531 = r9529 + r9530;
        double r9532 = r9528 / r9531;
        return r9532;
}

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 r9533, r9534, r9535, r9536, r9537, r9538, r9539, r9540, r9541, r9542, r9543;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9533);
        mpfr_init(r9534);
        mpfr_init(r9535);
        mpfr_init(r9536);
        mpfr_init(r9537);
        mpfr_init(r9538);
        mpfr_init(r9539);
        mpfr_init(r9540);
        mpfr_init(r9541);
        mpfr_init(r9542);
        mpfr_init(r9543);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9533, x_re, MPFR_RNDN);
        mpfr_set_d(r9534, y_re, MPFR_RNDN);
        mpfr_mul(r9535, r9533, r9534, MPFR_RNDN);
        mpfr_set_d(r9536, x_im, MPFR_RNDN);
        mpfr_set_d(r9537, y_im, MPFR_RNDN);
        mpfr_mul(r9538, r9536, r9537, MPFR_RNDN);
        mpfr_add(r9539, r9535, r9538, MPFR_RNDN);
        mpfr_mul(r9540, r9534, r9534, MPFR_RNDN);
        mpfr_mul(r9541, r9537, r9537, MPFR_RNDN);
        mpfr_add(r9542, r9540, r9541, MPFR_RNDN);
        mpfr_div(r9543, r9539, r9542, MPFR_RNDN);
        return mpfr_get_d(r9543, MPFR_RNDN);
}

static mpfr_t r9544, r9545, r9546, r9547, r9548, r9549, r9550, r9551, r9552, r9553, r9554;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9544);
        mpfr_init(r9545);
        mpfr_init(r9546);
        mpfr_init(r9547);
        mpfr_init(r9548);
        mpfr_init(r9549);
        mpfr_init(r9550);
        mpfr_init(r9551);
        mpfr_init(r9552);
        mpfr_init(r9553);
        mpfr_init(r9554);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9544, x_im, MPFR_RNDN);
        mpfr_set_d(r9545, y_im, MPFR_RNDN);
        mpfr_mul(r9546, r9544, r9545, MPFR_RNDN);
        mpfr_set_d(r9547, x_re, MPFR_RNDN);
        mpfr_set_d(r9548, y_re, MPFR_RNDN);
        mpfr_mul(r9549, r9547, r9548, MPFR_RNDN);
        mpfr_add(r9550, r9546, r9549, MPFR_RNDN);
        mpfr_mul(r9551, r9548, r9548, MPFR_RNDN);
        mpfr_mul(r9552, r9545, r9545, MPFR_RNDN);
        mpfr_add(r9553, r9551, r9552, MPFR_RNDN);
        mpfr_div(r9554, r9550, r9553, MPFR_RNDN);
        return mpfr_get_d(r9554, MPFR_RNDN);
}

static mpfr_t r9555, r9556, r9557, r9558, r9559, r9560, r9561, r9562, r9563, r9564, r9565;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9555);
        mpfr_init(r9556);
        mpfr_init(r9557);
        mpfr_init(r9558);
        mpfr_init(r9559);
        mpfr_init(r9560);
        mpfr_init(r9561);
        mpfr_init(r9562);
        mpfr_init(r9563);
        mpfr_init(r9564);
        mpfr_init(r9565);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9555, x_im, MPFR_RNDN);
        mpfr_set_d(r9556, y_im, MPFR_RNDN);
        mpfr_mul(r9557, r9555, r9556, MPFR_RNDN);
        mpfr_set_d(r9558, x_re, MPFR_RNDN);
        mpfr_set_d(r9559, y_re, MPFR_RNDN);
        mpfr_mul(r9560, r9558, r9559, MPFR_RNDN);
        mpfr_add(r9561, r9557, r9560, MPFR_RNDN);
        mpfr_mul(r9562, r9559, r9559, MPFR_RNDN);
        mpfr_mul(r9563, r9556, r9556, MPFR_RNDN);
        mpfr_add(r9564, r9562, r9563, MPFR_RNDN);
        mpfr_div(r9565, r9561, r9564, MPFR_RNDN);
        return mpfr_get_d(r9565, MPFR_RNDN);
}