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

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

double f_if(float a, float b, float c, float d) {
        float r26710 = b;
        float r26711 = c;
        float r26712 = r26710 * r26711;
        float r26713 = a;
        float r26714 = d;
        float r26715 = r26713 * r26714;
        float r26716 = r26712 - r26715;
        float r26717 = r26711 * r26711;
        float r26718 = r26714 * r26714;
        float r26719 = r26717 + r26718;
        float r26720 = r26716 / r26719;
        return r26720;
}

double f_id(double a, double b, double c, double d) {
        double r26721 = b;
        double r26722 = c;
        double r26723 = r26721 * r26722;
        double r26724 = a;
        double r26725 = d;
        double r26726 = r26724 * r26725;
        double r26727 = r26723 - r26726;
        double r26728 = r26722 * r26722;
        double r26729 = r26725 * r26725;
        double r26730 = r26728 + r26729;
        double r26731 = r26727 / r26730;
        return r26731;
}


double f_of(float a, float b, float c, float d) {
        float r26732 = c;
        float r26733 = -5.03571375883534e+131;
        bool r26734 = r26732 <= r26733;
        float r26735 = 1.5009880531967683e+116;
        bool r26736 = r26732 <= r26735;
        float r26737 = !r26736;
        bool r26738 = r26734 || r26737;
        float r26739 = b;
        float r26740 = r26739 / r26732;
        float r26741 = a;
        float r26742 = r26741 / r26732;
        float r26743 = d;
        float r26744 = r26743 / r26732;
        float r26745 = r26742 * r26744;
        float r26746 = r26740 - r26745;
        float r26747 = r26739 * r26732;
        float r26748 = r26743 * r26741;
        float r26749 = r26747 - r26748;
        float r26750 = r26732 * r26732;
        float r26751 = r26743 * r26743;
        float r26752 = r26750 + r26751;
        float r26753 = sqrt(r26752);
        float r26754 = r26749 / r26753;
        float r26755 = r26754 / r26753;
        float r26756 = r26738 ? r26746 : r26755;
        return r26756;
}

double f_od(double a, double b, double c, double d) {
        double r26757 = c;
        double r26758 = -5.03571375883534e+131;
        bool r26759 = r26757 <= r26758;
        double r26760 = 1.5009880531967683e+116;
        bool r26761 = r26757 <= r26760;
        double r26762 = !r26761;
        bool r26763 = r26759 || r26762;
        double r26764 = b;
        double r26765 = r26764 / r26757;
        double r26766 = a;
        double r26767 = r26766 / r26757;
        double r26768 = d;
        double r26769 = r26768 / r26757;
        double r26770 = r26767 * r26769;
        double r26771 = r26765 - r26770;
        double r26772 = r26764 * r26757;
        double r26773 = r26768 * r26766;
        double r26774 = r26772 - r26773;
        double r26775 = r26757 * r26757;
        double r26776 = r26768 * r26768;
        double r26777 = r26775 + r26776;
        double r26778 = sqrt(r26777);
        double r26779 = r26774 / r26778;
        double r26780 = r26779 / r26778;
        double r26781 = r26763 ? r26771 : r26780;
        return r26781;
}

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 r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26782);
        mpfr_init(r26783);
        mpfr_init(r26784);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init(r26787);
        mpfr_init(r26788);
        mpfr_init(r26789);
        mpfr_init(r26790);
        mpfr_init(r26791);
        mpfr_init(r26792);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26782, b, MPFR_RNDN);
        mpfr_set_d(r26783, c, MPFR_RNDN);
        mpfr_mul(r26784, r26782, r26783, MPFR_RNDN);
        mpfr_set_d(r26785, a, MPFR_RNDN);
        mpfr_set_d(r26786, d, MPFR_RNDN);
        mpfr_mul(r26787, r26785, r26786, MPFR_RNDN);
        mpfr_sub(r26788, r26784, r26787, MPFR_RNDN);
        mpfr_mul(r26789, r26783, r26783, MPFR_RNDN);
        mpfr_mul(r26790, r26786, r26786, MPFR_RNDN);
        mpfr_add(r26791, r26789, r26790, MPFR_RNDN);
        mpfr_div(r26792, r26788, r26791, MPFR_RNDN);
        return mpfr_get_d(r26792, MPFR_RNDN);
}

static mpfr_t r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807, r26808, r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26793);
        mpfr_init_set_str(r26794, "-5.03571375883534e+131", 10, MPFR_RNDN);
        mpfr_init(r26795);
        mpfr_init_set_str(r26796, "1.5009880531967683e+116", 10, MPFR_RNDN);
        mpfr_init(r26797);
        mpfr_init(r26798);
        mpfr_init(r26799);
        mpfr_init(r26800);
        mpfr_init(r26801);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init(r26804);
        mpfr_init(r26805);
        mpfr_init(r26806);
        mpfr_init(r26807);
        mpfr_init(r26808);
        mpfr_init(r26809);
        mpfr_init(r26810);
        mpfr_init(r26811);
        mpfr_init(r26812);
        mpfr_init(r26813);
        mpfr_init(r26814);
        mpfr_init(r26815);
        mpfr_init(r26816);
        mpfr_init(r26817);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26793, c, MPFR_RNDN);
        ;
        mpfr_set_si(r26795, mpfr_cmp(r26793, r26794) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r26797, mpfr_cmp(r26793, r26796) <= 0, MPFR_RNDN);
        mpfr_set_si(r26798, !mpfr_get_si(r26797, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r26799, mpfr_get_si(r26795, MPFR_RNDN) || mpfr_get_si(r26798, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r26800, b, MPFR_RNDN);
        mpfr_div(r26801, r26800, r26793, MPFR_RNDN);
        mpfr_set_d(r26802, a, MPFR_RNDN);
        mpfr_div(r26803, r26802, r26793, MPFR_RNDN);
        mpfr_set_d(r26804, d, MPFR_RNDN);
        mpfr_div(r26805, r26804, r26793, MPFR_RNDN);
        mpfr_mul(r26806, r26803, r26805, MPFR_RNDN);
        mpfr_sub(r26807, r26801, r26806, MPFR_RNDN);
        mpfr_mul(r26808, r26800, r26793, MPFR_RNDN);
        mpfr_mul(r26809, r26804, r26802, MPFR_RNDN);
        mpfr_sub(r26810, r26808, r26809, MPFR_RNDN);
        mpfr_mul(r26811, r26793, r26793, MPFR_RNDN);
        mpfr_mul(r26812, r26804, r26804, MPFR_RNDN);
        mpfr_add(r26813, r26811, r26812, MPFR_RNDN);
        mpfr_sqrt(r26814, r26813, MPFR_RNDN);
        mpfr_div(r26815, r26810, r26814, MPFR_RNDN);
        mpfr_div(r26816, r26815, r26814, MPFR_RNDN);
        if (mpfr_get_si(r26799, MPFR_RNDN)) { mpfr_set(r26817, r26807, MPFR_RNDN); } else { mpfr_set(r26817, r26816, MPFR_RNDN); };
        return mpfr_get_d(r26817, MPFR_RNDN);
}

static mpfr_t r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834, r26835, r26836, r26837, r26838, r26839, r26840, r26841, r26842;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26818);
        mpfr_init_set_str(r26819, "-5.03571375883534e+131", 10, MPFR_RNDN);
        mpfr_init(r26820);
        mpfr_init_set_str(r26821, "1.5009880531967683e+116", 10, MPFR_RNDN);
        mpfr_init(r26822);
        mpfr_init(r26823);
        mpfr_init(r26824);
        mpfr_init(r26825);
        mpfr_init(r26826);
        mpfr_init(r26827);
        mpfr_init(r26828);
        mpfr_init(r26829);
        mpfr_init(r26830);
        mpfr_init(r26831);
        mpfr_init(r26832);
        mpfr_init(r26833);
        mpfr_init(r26834);
        mpfr_init(r26835);
        mpfr_init(r26836);
        mpfr_init(r26837);
        mpfr_init(r26838);
        mpfr_init(r26839);
        mpfr_init(r26840);
        mpfr_init(r26841);
        mpfr_init(r26842);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26818, c, MPFR_RNDN);
        ;
        mpfr_set_si(r26820, mpfr_cmp(r26818, r26819) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r26822, mpfr_cmp(r26818, r26821) <= 0, MPFR_RNDN);
        mpfr_set_si(r26823, !mpfr_get_si(r26822, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r26824, mpfr_get_si(r26820, MPFR_RNDN) || mpfr_get_si(r26823, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r26825, b, MPFR_RNDN);
        mpfr_div(r26826, r26825, r26818, MPFR_RNDN);
        mpfr_set_d(r26827, a, MPFR_RNDN);
        mpfr_div(r26828, r26827, r26818, MPFR_RNDN);
        mpfr_set_d(r26829, d, MPFR_RNDN);
        mpfr_div(r26830, r26829, r26818, MPFR_RNDN);
        mpfr_mul(r26831, r26828, r26830, MPFR_RNDN);
        mpfr_sub(r26832, r26826, r26831, MPFR_RNDN);
        mpfr_mul(r26833, r26825, r26818, MPFR_RNDN);
        mpfr_mul(r26834, r26829, r26827, MPFR_RNDN);
        mpfr_sub(r26835, r26833, r26834, MPFR_RNDN);
        mpfr_mul(r26836, r26818, r26818, MPFR_RNDN);
        mpfr_mul(r26837, r26829, r26829, MPFR_RNDN);
        mpfr_add(r26838, r26836, r26837, MPFR_RNDN);
        mpfr_sqrt(r26839, r26838, MPFR_RNDN);
        mpfr_div(r26840, r26835, r26839, MPFR_RNDN);
        mpfr_div(r26841, r26840, r26839, MPFR_RNDN);
        if (mpfr_get_si(r26824, MPFR_RNDN)) { mpfr_set(r26842, r26832, MPFR_RNDN); } else { mpfr_set(r26842, r26841, MPFR_RNDN); };
        return mpfr_get_d(r26842, MPFR_RNDN);
}