#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 r9681 = a;
        float r9682 = c;
        float r9683 = r9681 * r9682;
        float r9684 = b;
        float r9685 = d;
        float r9686 = r9684 * r9685;
        float r9687 = r9683 + r9686;
        float r9688 = r9682 * r9682;
        float r9689 = r9685 * r9685;
        float r9690 = r9688 + r9689;
        float r9691 = r9687 / r9690;
        return r9691;
}

double f_id(double a, double b, double c, double d) {
        double r9692 = a;
        double r9693 = c;
        double r9694 = r9692 * r9693;
        double r9695 = b;
        double r9696 = d;
        double r9697 = r9695 * r9696;
        double r9698 = r9694 + r9697;
        double r9699 = r9693 * r9693;
        double r9700 = r9696 * r9696;
        double r9701 = r9699 + r9700;
        double r9702 = r9698 / r9701;
        return r9702;
}


double f_of(float a, float b, float c, float d) {
        float r9703 = d;
        float r9704 = -1.4362146765870036e+89f;
        bool r9705 = r9703 <= r9704;
        float r9706 = 1.0f;
        float r9707 = -r9706;
        float r9708 = c;
        float r9709 = hypot(r9703, r9708);
        float r9710 = r9707 / r9709;
        float r9711 = b;
        float r9712 = r9710 * r9711;
        float r9713 = 3.682673386505032e+207f;
        bool r9714 = r9703 <= r9713;
        float r9715 = a;
        float r9716 = r9711 * r9703;
        float r9717 = fma(r9708, r9715, r9716);
        float r9718 = r9717 / r9709;
        float r9719 = r9718 / r9709;
        float r9720 = r9711 / r9709;
        float r9721 = r9714 ? r9719 : r9720;
        float r9722 = r9705 ? r9712 : r9721;
        return r9722;
}

double f_od(double a, double b, double c, double d) {
        double r9723 = d;
        double r9724 = -1.4362146765870036e+89;
        bool r9725 = r9723 <= r9724;
        double r9726 = 1.0;
        double r9727 = -r9726;
        double r9728 = c;
        double r9729 = hypot(r9723, r9728);
        double r9730 = r9727 / r9729;
        double r9731 = b;
        double r9732 = r9730 * r9731;
        double r9733 = 3.682673386505032e+207;
        bool r9734 = r9723 <= r9733;
        double r9735 = a;
        double r9736 = r9731 * r9723;
        double r9737 = fma(r9728, r9735, r9736);
        double r9738 = r9737 / r9729;
        double r9739 = r9738 / r9729;
        double r9740 = r9731 / r9729;
        double r9741 = r9734 ? r9739 : r9740;
        double r9742 = r9725 ? r9732 : r9741;
        return r9742;
}

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 r9743, r9744, r9745, r9746, r9747, r9748, r9749, r9750, r9751, r9752, r9753;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9743);
        mpfr_init(r9744);
        mpfr_init(r9745);
        mpfr_init(r9746);
        mpfr_init(r9747);
        mpfr_init(r9748);
        mpfr_init(r9749);
        mpfr_init(r9750);
        mpfr_init(r9751);
        mpfr_init(r9752);
        mpfr_init(r9753);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r9743, a, MPFR_RNDN);
        mpfr_set_d(r9744, c, MPFR_RNDN);
        mpfr_mul(r9745, r9743, r9744, MPFR_RNDN);
        mpfr_set_d(r9746, b, MPFR_RNDN);
        mpfr_set_d(r9747, d, MPFR_RNDN);
        mpfr_mul(r9748, r9746, r9747, MPFR_RNDN);
        mpfr_add(r9749, r9745, r9748, MPFR_RNDN);
        mpfr_mul(r9750, r9744, r9744, MPFR_RNDN);
        mpfr_mul(r9751, r9747, r9747, MPFR_RNDN);
        mpfr_add(r9752, r9750, r9751, MPFR_RNDN);
        mpfr_div(r9753, r9749, r9752, MPFR_RNDN);
        return mpfr_get_d(r9753, MPFR_RNDN);
}

static mpfr_t r9754, r9755, r9756, r9757, r9758, r9759, r9760, r9761, r9762, r9763, r9764, r9765, r9766, r9767, r9768, r9769, r9770, r9771, r9772, r9773;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9754);
        mpfr_init_set_str(r9755, "-1.4362146765870036e+89", 10, MPFR_RNDN);
        mpfr_init(r9756);
        mpfr_init_set_str(r9757, "1", 10, MPFR_RNDN);
        mpfr_init(r9758);
        mpfr_init(r9759);
        mpfr_init(r9760);
        mpfr_init(r9761);
        mpfr_init(r9762);
        mpfr_init(r9763);
        mpfr_init_set_str(r9764, "3.682673386505032e+207", 10, MPFR_RNDN);
        mpfr_init(r9765);
        mpfr_init(r9766);
        mpfr_init(r9767);
        mpfr_init(r9768);
        mpfr_init(r9769);
        mpfr_init(r9770);
        mpfr_init(r9771);
        mpfr_init(r9772);
        mpfr_init(r9773);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r9754, d, MPFR_RNDN);
        ;
        mpfr_set_si(r9756, mpfr_cmp(r9754, r9755) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r9758, r9757, MPFR_RNDN);
        mpfr_set_d(r9759, c, MPFR_RNDN);
        mpfr_hypot(r9760, r9754, r9759, MPFR_RNDN);
        mpfr_div(r9761, r9758, r9760, MPFR_RNDN);
        mpfr_set_d(r9762, b, MPFR_RNDN);
        mpfr_mul(r9763, r9761, r9762, MPFR_RNDN);
        ;
        mpfr_set_si(r9765, mpfr_cmp(r9754, r9764) <= 0, MPFR_RNDN);
        mpfr_set_d(r9766, a, MPFR_RNDN);
        mpfr_mul(r9767, r9762, r9754, MPFR_RNDN);
        mpfr_fma(r9768, r9759, r9766, r9767, MPFR_RNDN);
        mpfr_div(r9769, r9768, r9760, MPFR_RNDN);
        mpfr_div(r9770, r9769, r9760, MPFR_RNDN);
        mpfr_div(r9771, r9762, r9760, MPFR_RNDN);
        if (mpfr_get_si(r9765, MPFR_RNDN)) { mpfr_set(r9772, r9770, MPFR_RNDN); } else { mpfr_set(r9772, r9771, MPFR_RNDN); };
        if (mpfr_get_si(r9756, MPFR_RNDN)) { mpfr_set(r9773, r9763, MPFR_RNDN); } else { mpfr_set(r9773, r9772, MPFR_RNDN); };
        return mpfr_get_d(r9773, MPFR_RNDN);
}

static mpfr_t r9774, r9775, r9776, r9777, r9778, r9779, r9780, r9781, r9782, r9783, r9784, r9785, r9786, r9787, r9788, r9789, r9790, r9791, r9792, r9793;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9774);
        mpfr_init_set_str(r9775, "-1.4362146765870036e+89", 10, MPFR_RNDN);
        mpfr_init(r9776);
        mpfr_init_set_str(r9777, "1", 10, MPFR_RNDN);
        mpfr_init(r9778);
        mpfr_init(r9779);
        mpfr_init(r9780);
        mpfr_init(r9781);
        mpfr_init(r9782);
        mpfr_init(r9783);
        mpfr_init_set_str(r9784, "3.682673386505032e+207", 10, MPFR_RNDN);
        mpfr_init(r9785);
        mpfr_init(r9786);
        mpfr_init(r9787);
        mpfr_init(r9788);
        mpfr_init(r9789);
        mpfr_init(r9790);
        mpfr_init(r9791);
        mpfr_init(r9792);
        mpfr_init(r9793);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r9774, d, MPFR_RNDN);
        ;
        mpfr_set_si(r9776, mpfr_cmp(r9774, r9775) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r9778, r9777, MPFR_RNDN);
        mpfr_set_d(r9779, c, MPFR_RNDN);
        mpfr_hypot(r9780, r9774, r9779, MPFR_RNDN);
        mpfr_div(r9781, r9778, r9780, MPFR_RNDN);
        mpfr_set_d(r9782, b, MPFR_RNDN);
        mpfr_mul(r9783, r9781, r9782, MPFR_RNDN);
        ;
        mpfr_set_si(r9785, mpfr_cmp(r9774, r9784) <= 0, MPFR_RNDN);
        mpfr_set_d(r9786, a, MPFR_RNDN);
        mpfr_mul(r9787, r9782, r9774, MPFR_RNDN);
        mpfr_fma(r9788, r9779, r9786, r9787, MPFR_RNDN);
        mpfr_div(r9789, r9788, r9780, MPFR_RNDN);
        mpfr_div(r9790, r9789, r9780, MPFR_RNDN);
        mpfr_div(r9791, r9782, r9780, MPFR_RNDN);
        if (mpfr_get_si(r9785, MPFR_RNDN)) { mpfr_set(r9792, r9790, MPFR_RNDN); } else { mpfr_set(r9792, r9791, MPFR_RNDN); };
        if (mpfr_get_si(r9776, MPFR_RNDN)) { mpfr_set(r9793, r9783, MPFR_RNDN); } else { mpfr_set(r9793, r9792, MPFR_RNDN); };
        return mpfr_get_d(r9793, MPFR_RNDN);
}