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

char *name = "Complex division, real part";

double f_if(float a, float b, float c, float d) {
        float r23545 = a;
        float r23546 = c;
        float r23547 = r23545 * r23546;
        float r23548 = b;
        float r23549 = d;
        float r23550 = r23548 * r23549;
        float r23551 = r23547 + r23550;
        float r23552 = r23546 * r23546;
        float r23553 = r23549 * r23549;
        float r23554 = r23552 + r23553;
        float r23555 = r23551 / r23554;
        return r23555;
}

double f_id(double a, double b, double c, double d) {
        double r23556 = a;
        double r23557 = c;
        double r23558 = r23556 * r23557;
        double r23559 = b;
        double r23560 = d;
        double r23561 = r23559 * r23560;
        double r23562 = r23558 + r23561;
        double r23563 = r23557 * r23557;
        double r23564 = r23560 * r23560;
        double r23565 = r23563 + r23564;
        double r23566 = r23562 / r23565;
        return r23566;
}


double f_of(float a, float b, float c, float d) {
        float r23567 = d;
        float r23568 = -1.8681789896549758e+96;
        bool r23569 = r23567 <= r23568;
        float r23570 = b;
        float r23571 = -r23570;
        float r23572 = c;
        float r23573 = hypot(r23572, r23567);
        float r23574 = r23571 / r23573;
        float r23575 = 1.2175863528347458e+124;
        bool r23576 = r23567 <= r23575;
        float r23577 = a;
        float r23578 = r23577 * r23572;
        float r23579 = fma(r23567, r23570, r23578);
        float r23580 = r23579 / r23573;
        float r23581 = r23580 / r23573;
        float r23582 = r23572 / r23567;
        float r23583 = fma(r23582, r23577, r23570);
        float r23584 = r23583 / r23573;
        float r23585 = r23576 ? r23581 : r23584;
        float r23586 = r23569 ? r23574 : r23585;
        return r23586;
}

double f_od(double a, double b, double c, double d) {
        double r23587 = d;
        double r23588 = -1.8681789896549758e+96;
        bool r23589 = r23587 <= r23588;
        double r23590 = b;
        double r23591 = -r23590;
        double r23592 = c;
        double r23593 = hypot(r23592, r23587);
        double r23594 = r23591 / r23593;
        double r23595 = 1.2175863528347458e+124;
        bool r23596 = r23587 <= r23595;
        double r23597 = a;
        double r23598 = r23597 * r23592;
        double r23599 = fma(r23587, r23590, r23598);
        double r23600 = r23599 / r23593;
        double r23601 = r23600 / r23593;
        double r23602 = r23592 / r23587;
        double r23603 = fma(r23602, r23597, r23590);
        double r23604 = r23603 / r23593;
        double r23605 = r23596 ? r23601 : r23604;
        double r23606 = r23589 ? r23594 : r23605;
        return r23606;
}

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 r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r23607);
        mpfr_init(r23608);
        mpfr_init(r23609);
        mpfr_init(r23610);
        mpfr_init(r23611);
        mpfr_init(r23612);
        mpfr_init(r23613);
        mpfr_init(r23614);
        mpfr_init(r23615);
        mpfr_init(r23616);
        mpfr_init(r23617);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r23607, a, MPFR_RNDN);
        mpfr_set_d(r23608, c, MPFR_RNDN);
        mpfr_mul(r23609, r23607, r23608, MPFR_RNDN);
        mpfr_set_d(r23610, b, MPFR_RNDN);
        mpfr_set_d(r23611, d, MPFR_RNDN);
        mpfr_mul(r23612, r23610, r23611, MPFR_RNDN);
        mpfr_add(r23613, r23609, r23612, MPFR_RNDN);
        mpfr_mul(r23614, r23608, r23608, MPFR_RNDN);
        mpfr_mul(r23615, r23611, r23611, MPFR_RNDN);
        mpfr_add(r23616, r23614, r23615, MPFR_RNDN);
        mpfr_div(r23617, r23613, r23616, MPFR_RNDN);
        return mpfr_get_d(r23617, MPFR_RNDN);
}

static mpfr_t r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r23618);
        mpfr_init_set_str(r23619, "-1.8681789896549758e+96", 10, MPFR_RNDN);
        mpfr_init(r23620);
        mpfr_init(r23621);
        mpfr_init(r23622);
        mpfr_init(r23623);
        mpfr_init(r23624);
        mpfr_init(r23625);
        mpfr_init_set_str(r23626, "1.2175863528347458e+124", 10, MPFR_RNDN);
        mpfr_init(r23627);
        mpfr_init(r23628);
        mpfr_init(r23629);
        mpfr_init(r23630);
        mpfr_init(r23631);
        mpfr_init(r23632);
        mpfr_init(r23633);
        mpfr_init(r23634);
        mpfr_init(r23635);
        mpfr_init(r23636);
        mpfr_init(r23637);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r23618, d, MPFR_RNDN);
        ;
        mpfr_set_si(r23620, mpfr_cmp(r23618, r23619) <= 0, MPFR_RNDN);
        mpfr_set_d(r23621, b, MPFR_RNDN);
        mpfr_neg(r23622, r23621, MPFR_RNDN);
        mpfr_set_d(r23623, c, MPFR_RNDN);
        mpfr_hypot(r23624, r23623, r23618, MPFR_RNDN);
        mpfr_div(r23625, r23622, r23624, MPFR_RNDN);
        ;
        mpfr_set_si(r23627, mpfr_cmp(r23618, r23626) <= 0, MPFR_RNDN);
        mpfr_set_d(r23628, a, MPFR_RNDN);
        mpfr_mul(r23629, r23628, r23623, MPFR_RNDN);
        mpfr_fma(r23630, r23618, r23621, r23629, MPFR_RNDN);
        mpfr_div(r23631, r23630, r23624, MPFR_RNDN);
        mpfr_div(r23632, r23631, r23624, MPFR_RNDN);
        mpfr_div(r23633, r23623, r23618, MPFR_RNDN);
        mpfr_fma(r23634, r23633, r23628, r23621, MPFR_RNDN);
        mpfr_div(r23635, r23634, r23624, MPFR_RNDN);
        if (mpfr_get_si(r23627, MPFR_RNDN)) { mpfr_set(r23636, r23632, MPFR_RNDN); } else { mpfr_set(r23636, r23635, MPFR_RNDN); };
        if (mpfr_get_si(r23620, MPFR_RNDN)) { mpfr_set(r23637, r23625, MPFR_RNDN); } else { mpfr_set(r23637, r23636, MPFR_RNDN); };
        return mpfr_get_d(r23637, MPFR_RNDN);
}

static mpfr_t r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648, r23649, r23650, r23651, r23652, r23653, r23654, r23655, r23656, r23657;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r23638);
        mpfr_init_set_str(r23639, "-1.8681789896549758e+96", 10, MPFR_RNDN);
        mpfr_init(r23640);
        mpfr_init(r23641);
        mpfr_init(r23642);
        mpfr_init(r23643);
        mpfr_init(r23644);
        mpfr_init(r23645);
        mpfr_init_set_str(r23646, "1.2175863528347458e+124", 10, MPFR_RNDN);
        mpfr_init(r23647);
        mpfr_init(r23648);
        mpfr_init(r23649);
        mpfr_init(r23650);
        mpfr_init(r23651);
        mpfr_init(r23652);
        mpfr_init(r23653);
        mpfr_init(r23654);
        mpfr_init(r23655);
        mpfr_init(r23656);
        mpfr_init(r23657);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r23638, d, MPFR_RNDN);
        ;
        mpfr_set_si(r23640, mpfr_cmp(r23638, r23639) <= 0, MPFR_RNDN);
        mpfr_set_d(r23641, b, MPFR_RNDN);
        mpfr_neg(r23642, r23641, MPFR_RNDN);
        mpfr_set_d(r23643, c, MPFR_RNDN);
        mpfr_hypot(r23644, r23643, r23638, MPFR_RNDN);
        mpfr_div(r23645, r23642, r23644, MPFR_RNDN);
        ;
        mpfr_set_si(r23647, mpfr_cmp(r23638, r23646) <= 0, MPFR_RNDN);
        mpfr_set_d(r23648, a, MPFR_RNDN);
        mpfr_mul(r23649, r23648, r23643, MPFR_RNDN);
        mpfr_fma(r23650, r23638, r23641, r23649, MPFR_RNDN);
        mpfr_div(r23651, r23650, r23644, MPFR_RNDN);
        mpfr_div(r23652, r23651, r23644, MPFR_RNDN);
        mpfr_div(r23653, r23643, r23638, MPFR_RNDN);
        mpfr_fma(r23654, r23653, r23648, r23641, MPFR_RNDN);
        mpfr_div(r23655, r23654, r23644, MPFR_RNDN);
        if (mpfr_get_si(r23647, MPFR_RNDN)) { mpfr_set(r23656, r23652, MPFR_RNDN); } else { mpfr_set(r23656, r23655, MPFR_RNDN); };
        if (mpfr_get_si(r23640, MPFR_RNDN)) { mpfr_set(r23657, r23645, MPFR_RNDN); } else { mpfr_set(r23657, r23656, MPFR_RNDN); };
        return mpfr_get_d(r23657, MPFR_RNDN);
}