#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 r17127 = b;
        float r17128 = c;
        float r17129 = r17127 * r17128;
        float r17130 = a;
        float r17131 = d;
        float r17132 = r17130 * r17131;
        float r17133 = r17129 - r17132;
        float r17134 = r17128 * r17128;
        float r17135 = r17131 * r17131;
        float r17136 = r17134 + r17135;
        float r17137 = r17133 / r17136;
        return r17137;
}

double f_id(double a, double b, double c, double d) {
        double r17138 = b;
        double r17139 = c;
        double r17140 = r17138 * r17139;
        double r17141 = a;
        double r17142 = d;
        double r17143 = r17141 * r17142;
        double r17144 = r17140 - r17143;
        double r17145 = r17139 * r17139;
        double r17146 = r17142 * r17142;
        double r17147 = r17145 + r17146;
        double r17148 = r17144 / r17147;
        return r17148;
}


double f_of(float a, float b, float c, float d) {
        float r17149 = d;
        float r17150 = -2.0753310964069658e+55f;
        bool r17151 = r17149 <= r17150;
        float r17152 = b;
        float r17153 = c;
        float r17154 = r17152 * r17153;
        float r17155 = r17153 * r17153;
        float r17156 = r17149 * r17149;
        float r17157 = r17155 + r17156;
        float r17158 = r17154 / r17157;
        float r17159 = a;
        float r17160 = r17155 / r17149;
        float r17161 = r17149 + r17160;
        float r17162 = r17159 / r17161;
        float r17163 = r17158 - r17162;
        float r17164 = 2.9373351910676443e+99f;
        bool r17165 = r17149 <= r17164;
        float r17166 = r17156 / r17153;
        float r17167 = r17166 + r17153;
        float r17168 = r17152 / r17167;
        float r17169 = r17159 * r17149;
        float r17170 = r17153 * r17153;
        float r17171 = r17156 + r17170;
        float r17172 = r17169 / r17171;
        float r17173 = r17168 - r17172;
        float r17174 = r17165 ? r17173 : r17163;
        float r17175 = r17151 ? r17163 : r17174;
        return r17175;
}

double f_od(double a, double b, double c, double d) {
        double r17176 = d;
        double r17177 = -2.0753310964069658e+55;
        bool r17178 = r17176 <= r17177;
        double r17179 = b;
        double r17180 = c;
        double r17181 = r17179 * r17180;
        double r17182 = r17180 * r17180;
        double r17183 = r17176 * r17176;
        double r17184 = r17182 + r17183;
        double r17185 = r17181 / r17184;
        double r17186 = a;
        double r17187 = r17182 / r17176;
        double r17188 = r17176 + r17187;
        double r17189 = r17186 / r17188;
        double r17190 = r17185 - r17189;
        double r17191 = 2.9373351910676443e+99;
        bool r17192 = r17176 <= r17191;
        double r17193 = r17183 / r17180;
        double r17194 = r17193 + r17180;
        double r17195 = r17179 / r17194;
        double r17196 = r17186 * r17176;
        double r17197 = r17180 * r17180;
        double r17198 = r17183 + r17197;
        double r17199 = r17196 / r17198;
        double r17200 = r17195 - r17199;
        double r17201 = r17192 ? r17200 : r17190;
        double r17202 = r17178 ? r17190 : r17201;
        return r17202;
}

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 r17203, r17204, r17205, r17206, r17207, r17208, r17209, r17210, r17211, r17212, r17213;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17203);
        mpfr_init(r17204);
        mpfr_init(r17205);
        mpfr_init(r17206);
        mpfr_init(r17207);
        mpfr_init(r17208);
        mpfr_init(r17209);
        mpfr_init(r17210);
        mpfr_init(r17211);
        mpfr_init(r17212);
        mpfr_init(r17213);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r17203, b, MPFR_RNDN);
        mpfr_set_d(r17204, c, MPFR_RNDN);
        mpfr_mul(r17205, r17203, r17204, MPFR_RNDN);
        mpfr_set_d(r17206, a, MPFR_RNDN);
        mpfr_set_d(r17207, d, MPFR_RNDN);
        mpfr_mul(r17208, r17206, r17207, MPFR_RNDN);
        mpfr_sub(r17209, r17205, r17208, MPFR_RNDN);
        mpfr_sqr(r17210, r17204, MPFR_RNDN);
        mpfr_sqr(r17211, r17207, MPFR_RNDN);
        mpfr_add(r17212, r17210, r17211, MPFR_RNDN);
        mpfr_div(r17213, r17209, r17212, MPFR_RNDN);
        return mpfr_get_d(r17213, MPFR_RNDN);
}

static mpfr_t r17214, r17215, r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229, r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17214);
        mpfr_init_set_str(r17215, "-2.0753310964069658e+55", 10, MPFR_RNDN);
        mpfr_init(r17216);
        mpfr_init(r17217);
        mpfr_init(r17218);
        mpfr_init(r17219);
        mpfr_init(r17220);
        mpfr_init(r17221);
        mpfr_init(r17222);
        mpfr_init(r17223);
        mpfr_init(r17224);
        mpfr_init(r17225);
        mpfr_init(r17226);
        mpfr_init(r17227);
        mpfr_init(r17228);
        mpfr_init_set_str(r17229, "2.9373351910676443e+99", 10, MPFR_RNDN);
        mpfr_init(r17230);
        mpfr_init(r17231);
        mpfr_init(r17232);
        mpfr_init(r17233);
        mpfr_init(r17234);
        mpfr_init(r17235);
        mpfr_init(r17236);
        mpfr_init(r17237);
        mpfr_init(r17238);
        mpfr_init(r17239);
        mpfr_init(r17240);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r17214, d, MPFR_RNDN);
        ;
        mpfr_set_si(r17216, mpfr_cmp(r17214, r17215) <= 0, MPFR_RNDN);
        mpfr_set_d(r17217, b, MPFR_RNDN);
        mpfr_set_d(r17218, c, MPFR_RNDN);
        mpfr_mul(r17219, r17217, r17218, MPFR_RNDN);
        mpfr_sqr(r17220, r17218, MPFR_RNDN);
        mpfr_mul(r17221, r17214, r17214, MPFR_RNDN);
        mpfr_add(r17222, r17220, r17221, MPFR_RNDN);
        mpfr_div(r17223, r17219, r17222, MPFR_RNDN);
        mpfr_set_d(r17224, a, MPFR_RNDN);
        mpfr_div(r17225, r17220, r17214, MPFR_RNDN);
        mpfr_add(r17226, r17214, r17225, MPFR_RNDN);
        mpfr_div(r17227, r17224, r17226, MPFR_RNDN);
        mpfr_sub(r17228, r17223, r17227, MPFR_RNDN);
        ;
        mpfr_set_si(r17230, mpfr_cmp(r17214, r17229) <= 0, MPFR_RNDN);
        mpfr_div(r17231, r17221, r17218, MPFR_RNDN);
        mpfr_add(r17232, r17231, r17218, MPFR_RNDN);
        mpfr_div(r17233, r17217, r17232, MPFR_RNDN);
        mpfr_mul(r17234, r17224, r17214, MPFR_RNDN);
        mpfr_mul(r17235, r17218, r17218, MPFR_RNDN);
        mpfr_add(r17236, r17221, r17235, MPFR_RNDN);
        mpfr_div(r17237, r17234, r17236, MPFR_RNDN);
        mpfr_sub(r17238, r17233, r17237, MPFR_RNDN);
        if (mpfr_get_si(r17230, MPFR_RNDN)) { mpfr_set(r17239, r17238, MPFR_RNDN); } else { mpfr_set(r17239, r17228, MPFR_RNDN); };
        if (mpfr_get_si(r17216, MPFR_RNDN)) { mpfr_set(r17240, r17228, MPFR_RNDN); } else { mpfr_set(r17240, r17239, MPFR_RNDN); };
        return mpfr_get_d(r17240, MPFR_RNDN);
}

static mpfr_t r17241, r17242, r17243, r17244, r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257, r17258, r17259, r17260, r17261, r17262, r17263, r17264, r17265, r17266, r17267;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17241);
        mpfr_init_set_str(r17242, "-2.0753310964069658e+55", 10, MPFR_RNDN);
        mpfr_init(r17243);
        mpfr_init(r17244);
        mpfr_init(r17245);
        mpfr_init(r17246);
        mpfr_init(r17247);
        mpfr_init(r17248);
        mpfr_init(r17249);
        mpfr_init(r17250);
        mpfr_init(r17251);
        mpfr_init(r17252);
        mpfr_init(r17253);
        mpfr_init(r17254);
        mpfr_init(r17255);
        mpfr_init_set_str(r17256, "2.9373351910676443e+99", 10, MPFR_RNDN);
        mpfr_init(r17257);
        mpfr_init(r17258);
        mpfr_init(r17259);
        mpfr_init(r17260);
        mpfr_init(r17261);
        mpfr_init(r17262);
        mpfr_init(r17263);
        mpfr_init(r17264);
        mpfr_init(r17265);
        mpfr_init(r17266);
        mpfr_init(r17267);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r17241, d, MPFR_RNDN);
        ;
        mpfr_set_si(r17243, mpfr_cmp(r17241, r17242) <= 0, MPFR_RNDN);
        mpfr_set_d(r17244, b, MPFR_RNDN);
        mpfr_set_d(r17245, c, MPFR_RNDN);
        mpfr_mul(r17246, r17244, r17245, MPFR_RNDN);
        mpfr_sqr(r17247, r17245, MPFR_RNDN);
        mpfr_mul(r17248, r17241, r17241, MPFR_RNDN);
        mpfr_add(r17249, r17247, r17248, MPFR_RNDN);
        mpfr_div(r17250, r17246, r17249, MPFR_RNDN);
        mpfr_set_d(r17251, a, MPFR_RNDN);
        mpfr_div(r17252, r17247, r17241, MPFR_RNDN);
        mpfr_add(r17253, r17241, r17252, MPFR_RNDN);
        mpfr_div(r17254, r17251, r17253, MPFR_RNDN);
        mpfr_sub(r17255, r17250, r17254, MPFR_RNDN);
        ;
        mpfr_set_si(r17257, mpfr_cmp(r17241, r17256) <= 0, MPFR_RNDN);
        mpfr_div(r17258, r17248, r17245, MPFR_RNDN);
        mpfr_add(r17259, r17258, r17245, MPFR_RNDN);
        mpfr_div(r17260, r17244, r17259, MPFR_RNDN);
        mpfr_mul(r17261, r17251, r17241, MPFR_RNDN);
        mpfr_mul(r17262, r17245, r17245, MPFR_RNDN);
        mpfr_add(r17263, r17248, r17262, MPFR_RNDN);
        mpfr_div(r17264, r17261, r17263, MPFR_RNDN);
        mpfr_sub(r17265, r17260, r17264, MPFR_RNDN);
        if (mpfr_get_si(r17257, MPFR_RNDN)) { mpfr_set(r17266, r17265, MPFR_RNDN); } else { mpfr_set(r17266, r17255, MPFR_RNDN); };
        if (mpfr_get_si(r17243, MPFR_RNDN)) { mpfr_set(r17267, r17255, MPFR_RNDN); } else { mpfr_set(r17267, r17266, MPFR_RNDN); };
        return mpfr_get_d(r17267, MPFR_RNDN);
}