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

char *name = "_divideComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r9566 = x_im;
        float r9567 = y_re;
        float r9568 = r9566 * r9567;
        float r9569 = x_re;
        float r9570 = y_im;
        float r9571 = r9569 * r9570;
        float r9572 = r9568 - r9571;
        float r9573 = r9567 * r9567;
        float r9574 = r9570 * r9570;
        float r9575 = r9573 + r9574;
        float r9576 = r9572 / r9575;
        return r9576;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r9577 = x_im;
        double r9578 = y_re;
        double r9579 = r9577 * r9578;
        double r9580 = x_re;
        double r9581 = y_im;
        double r9582 = r9580 * r9581;
        double r9583 = r9579 - r9582;
        double r9584 = r9578 * r9578;
        double r9585 = r9581 * r9581;
        double r9586 = r9584 + r9585;
        double r9587 = r9583 / r9586;
        return r9587;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r9588 = x_im;
        float r9589 = y_re;
        float r9590 = r9588 * r9589;
        float r9591 = x_re;
        float r9592 = y_im;
        float r9593 = r9591 * r9592;
        float r9594 = r9590 - r9593;
        float r9595 = r9589 * r9589;
        float r9596 = r9592 * r9592;
        float r9597 = r9595 + r9596;
        float r9598 = sqrt(r9597);
        float r9599 = r9594 / r9598;
        float r9600 = r9599 / r9598;
        return r9600;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r9601 = x_im;
        double r9602 = y_re;
        double r9603 = r9601 * r9602;
        double r9604 = x_re;
        double r9605 = y_im;
        double r9606 = r9604 * r9605;
        double r9607 = r9603 - r9606;
        double r9608 = r9602 * r9602;
        double r9609 = r9605 * r9605;
        double r9610 = r9608 + r9609;
        double r9611 = sqrt(r9610);
        double r9612 = r9607 / r9611;
        double r9613 = r9612 / r9611;
        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, r9620, r9621, r9622, r9623, r9624;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9614);
        mpfr_init(r9615);
        mpfr_init(r9616);
        mpfr_init(r9617);
        mpfr_init(r9618);
        mpfr_init(r9619);
        mpfr_init(r9620);
        mpfr_init(r9621);
        mpfr_init(r9622);
        mpfr_init(r9623);
        mpfr_init(r9624);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9614, x_im, MPFR_RNDN);
        mpfr_set_d(r9615, y_re, MPFR_RNDN);
        mpfr_mul(r9616, r9614, r9615, MPFR_RNDN);
        mpfr_set_d(r9617, x_re, MPFR_RNDN);
        mpfr_set_d(r9618, y_im, MPFR_RNDN);
        mpfr_mul(r9619, r9617, r9618, MPFR_RNDN);
        mpfr_sub(r9620, r9616, r9619, MPFR_RNDN);
        mpfr_mul(r9621, r9615, r9615, MPFR_RNDN);
        mpfr_mul(r9622, r9618, r9618, MPFR_RNDN);
        mpfr_add(r9623, r9621, r9622, MPFR_RNDN);
        mpfr_div(r9624, r9620, r9623, MPFR_RNDN);
        return mpfr_get_d(r9624, MPFR_RNDN);
}

static mpfr_t r9625, r9626, r9627, r9628, r9629, r9630, r9631, r9632, r9633, r9634, r9635, r9636, r9637;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        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);
        mpfr_init(r9636);
        mpfr_init(r9637);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9625, x_im, MPFR_RNDN);
        mpfr_set_d(r9626, y_re, MPFR_RNDN);
        mpfr_mul(r9627, r9625, r9626, MPFR_RNDN);
        mpfr_set_d(r9628, x_re, MPFR_RNDN);
        mpfr_set_d(r9629, y_im, MPFR_RNDN);
        mpfr_mul(r9630, r9628, r9629, MPFR_RNDN);
        mpfr_sub(r9631, r9627, r9630, MPFR_RNDN);
        mpfr_mul(r9632, r9626, r9626, MPFR_RNDN);
        mpfr_mul(r9633, r9629, r9629, MPFR_RNDN);
        mpfr_add(r9634, r9632, r9633, MPFR_RNDN);
        mpfr_sqrt(r9635, r9634, MPFR_RNDN);
        mpfr_div(r9636, r9631, r9635, MPFR_RNDN);
        mpfr_div(r9637, r9636, r9635, MPFR_RNDN);
        return mpfr_get_d(r9637, MPFR_RNDN);
}

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

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        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);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r9638, x_im, MPFR_RNDN);
        mpfr_set_d(r9639, y_re, MPFR_RNDN);
        mpfr_mul(r9640, r9638, r9639, MPFR_RNDN);
        mpfr_set_d(r9641, x_re, MPFR_RNDN);
        mpfr_set_d(r9642, y_im, MPFR_RNDN);
        mpfr_mul(r9643, r9641, r9642, MPFR_RNDN);
        mpfr_sub(r9644, r9640, r9643, MPFR_RNDN);
        mpfr_mul(r9645, r9639, r9639, MPFR_RNDN);
        mpfr_mul(r9646, r9642, r9642, MPFR_RNDN);
        mpfr_add(r9647, r9645, r9646, MPFR_RNDN);
        mpfr_sqrt(r9648, r9647, MPFR_RNDN);
        mpfr_div(r9649, r9644, r9648, MPFR_RNDN);
        mpfr_div(r9650, r9649, r9648, MPFR_RNDN);
        return mpfr_get_d(r9650, MPFR_RNDN);
}