#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 r17151 = b;
        float r17152 = c;
        float r17153 = r17151 * r17152;
        float r17154 = a;
        float r17155 = d;
        float r17156 = r17154 * r17155;
        float r17157 = r17153 - r17156;
        float r17158 = r17152 * r17152;
        float r17159 = r17155 * r17155;
        float r17160 = r17158 + r17159;
        float r17161 = r17157 / r17160;
        return r17161;
}

double f_id(double a, double b, double c, double d) {
        double r17162 = b;
        double r17163 = c;
        double r17164 = r17162 * r17163;
        double r17165 = a;
        double r17166 = d;
        double r17167 = r17165 * r17166;
        double r17168 = r17164 - r17167;
        double r17169 = r17163 * r17163;
        double r17170 = r17166 * r17166;
        double r17171 = r17169 + r17170;
        double r17172 = r17168 / r17171;
        return r17172;
}


double f_of(float a, float b, float c, float d) {
        float r17173 = d;
        float r17174 = -4.347783524512778e+20f;
        bool r17175 = r17173 <= r17174;
        float r17176 = a;
        float r17177 = r17176 / r17173;
        float r17178 = -r17177;
        float r17179 = -8.791165892236865e-22f;
        bool r17180 = r17173 <= r17179;
        float r17181 = b;
        float r17182 = c;
        float r17183 = r17181 * r17182;
        float r17184 = r17176 * r17173;
        float r17185 = r17183 - r17184;
        float r17186 = r17182 * r17182;
        float r17187 = r17173 * r17173;
        float r17188 = r17186 + r17187;
        float r17189 = sqrt(r17188);
        float r17190 = r17189 * r17189;
        float r17191 = r17185 / r17190;
        float r17192 = 10322673664.0f;
        bool r17193 = r17173 <= r17192;
        float r17194 = r17181 / r17182;
        float r17195 = r17173 / r17182;
        float r17196 = r17176 / r17182;
        float r17197 = r17195 * r17196;
        float r17198 = r17194 - r17197;
        float r17199 = r17193 ? r17198 : r17178;
        float r17200 = r17180 ? r17191 : r17199;
        float r17201 = r17175 ? r17178 : r17200;
        return r17201;
}

double f_od(double a, double b, double c, double d) {
        double r17202 = d;
        double r17203 = -4.347783524512778e+20;
        bool r17204 = r17202 <= r17203;
        double r17205 = a;
        double r17206 = r17205 / r17202;
        double r17207 = -r17206;
        double r17208 = -8.791165892236865e-22;
        bool r17209 = r17202 <= r17208;
        double r17210 = b;
        double r17211 = c;
        double r17212 = r17210 * r17211;
        double r17213 = r17205 * r17202;
        double r17214 = r17212 - r17213;
        double r17215 = r17211 * r17211;
        double r17216 = r17202 * r17202;
        double r17217 = r17215 + r17216;
        double r17218 = sqrt(r17217);
        double r17219 = r17218 * r17218;
        double r17220 = r17214 / r17219;
        double r17221 = 10322673664.0;
        bool r17222 = r17202 <= r17221;
        double r17223 = r17210 / r17211;
        double r17224 = r17202 / r17211;
        double r17225 = r17205 / r17211;
        double r17226 = r17224 * r17225;
        double r17227 = r17223 - r17226;
        double r17228 = r17222 ? r17227 : r17207;
        double r17229 = r17209 ? r17220 : r17228;
        double r17230 = r17204 ? r17207 : r17229;
        return r17230;
}

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 r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240, r17241;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        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);
        mpfr_init(r17241);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r17231, b, MPFR_RNDN);
        mpfr_set_d(r17232, c, MPFR_RNDN);
        mpfr_mul(r17233, r17231, r17232, MPFR_RNDN);
        mpfr_set_d(r17234, a, MPFR_RNDN);
        mpfr_set_d(r17235, d, MPFR_RNDN);
        mpfr_mul(r17236, r17234, r17235, MPFR_RNDN);
        mpfr_sub(r17237, r17233, r17236, MPFR_RNDN);
        mpfr_sqr(r17238, r17232, MPFR_RNDN);
        mpfr_sqr(r17239, r17235, MPFR_RNDN);
        mpfr_add(r17240, r17238, r17239, MPFR_RNDN);
        mpfr_div(r17241, r17237, r17240, MPFR_RNDN);
        return mpfr_get_d(r17241, MPFR_RNDN);
}

static mpfr_t 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, r17268, r17269, r17270;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17242);
        mpfr_init_set_str(r17243, "-4.3477835f+20", 10, MPFR_RNDN);
        mpfr_init(r17244);
        mpfr_init(r17245);
        mpfr_init(r17246);
        mpfr_init(r17247);
        mpfr_init_set_str(r17248, "-8.791166f-22", 10, MPFR_RNDN);
        mpfr_init(r17249);
        mpfr_init(r17250);
        mpfr_init(r17251);
        mpfr_init(r17252);
        mpfr_init(r17253);
        mpfr_init(r17254);
        mpfr_init(r17255);
        mpfr_init(r17256);
        mpfr_init(r17257);
        mpfr_init(r17258);
        mpfr_init(r17259);
        mpfr_init(r17260);
        mpfr_init_set_str(r17261, "1.0322674f+10", 10, MPFR_RNDN);
        mpfr_init(r17262);
        mpfr_init(r17263);
        mpfr_init(r17264);
        mpfr_init(r17265);
        mpfr_init(r17266);
        mpfr_init(r17267);
        mpfr_init(r17268);
        mpfr_init(r17269);
        mpfr_init(r17270);
}

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

static mpfr_t r17271, r17272, r17273, r17274, r17275, r17276, r17277, r17278, r17279, r17280, r17281, r17282, r17283, r17284, r17285, r17286, r17287, r17288, r17289, r17290, r17291, r17292, r17293, r17294, r17295, r17296, r17297, r17298, r17299;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17271);
        mpfr_init_set_str(r17272, "-4.3477835f+20", 10, MPFR_RNDN);
        mpfr_init(r17273);
        mpfr_init(r17274);
        mpfr_init(r17275);
        mpfr_init(r17276);
        mpfr_init_set_str(r17277, "-8.791166f-22", 10, MPFR_RNDN);
        mpfr_init(r17278);
        mpfr_init(r17279);
        mpfr_init(r17280);
        mpfr_init(r17281);
        mpfr_init(r17282);
        mpfr_init(r17283);
        mpfr_init(r17284);
        mpfr_init(r17285);
        mpfr_init(r17286);
        mpfr_init(r17287);
        mpfr_init(r17288);
        mpfr_init(r17289);
        mpfr_init_set_str(r17290, "1.0322674f+10", 10, MPFR_RNDN);
        mpfr_init(r17291);
        mpfr_init(r17292);
        mpfr_init(r17293);
        mpfr_init(r17294);
        mpfr_init(r17295);
        mpfr_init(r17296);
        mpfr_init(r17297);
        mpfr_init(r17298);
        mpfr_init(r17299);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r17271, d, MPFR_RNDN);
        ;
        mpfr_set_si(r17273, mpfr_cmp(r17271, r17272) <= 0, MPFR_RNDN);
        mpfr_set_d(r17274, a, MPFR_RNDN);
        mpfr_div(r17275, r17274, r17271, MPFR_RNDN);
        mpfr_neg(r17276, r17275, MPFR_RNDN);
        ;
        mpfr_set_si(r17278, mpfr_cmp(r17271, r17277) <= 0, MPFR_RNDN);
        mpfr_set_d(r17279, b, MPFR_RNDN);
        mpfr_set_d(r17280, c, MPFR_RNDN);
        mpfr_mul(r17281, r17279, r17280, MPFR_RNDN);
        mpfr_mul(r17282, r17274, r17271, MPFR_RNDN);
        mpfr_sub(r17283, r17281, r17282, MPFR_RNDN);
        mpfr_sqr(r17284, r17280, MPFR_RNDN);
        mpfr_sqr(r17285, r17271, MPFR_RNDN);
        mpfr_add(r17286, r17284, r17285, MPFR_RNDN);
        mpfr_sqrt(r17287, r17286, MPFR_RNDN);
        mpfr_sqr(r17288, r17287, MPFR_RNDN);
        mpfr_div(r17289, r17283, r17288, MPFR_RNDN);
        ;
        mpfr_set_si(r17291, mpfr_cmp(r17271, r17290) <= 0, MPFR_RNDN);
        mpfr_div(r17292, r17279, r17280, MPFR_RNDN);
        mpfr_div(r17293, r17271, r17280, MPFR_RNDN);
        mpfr_div(r17294, r17274, r17280, MPFR_RNDN);
        mpfr_mul(r17295, r17293, r17294, MPFR_RNDN);
        mpfr_sub(r17296, r17292, r17295, MPFR_RNDN);
        if (mpfr_get_si(r17291, MPFR_RNDN)) { mpfr_set(r17297, r17296, MPFR_RNDN); } else { mpfr_set(r17297, r17276, MPFR_RNDN); };
        if (mpfr_get_si(r17278, MPFR_RNDN)) { mpfr_set(r17298, r17289, MPFR_RNDN); } else { mpfr_set(r17298, r17297, MPFR_RNDN); };
        if (mpfr_get_si(r17273, MPFR_RNDN)) { mpfr_set(r17299, r17276, MPFR_RNDN); } else { mpfr_set(r17299, r17298, MPFR_RNDN); };
        return mpfr_get_d(r17299, MPFR_RNDN);
}