#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 r26231 = b;
        float r26232 = c;
        float r26233 = r26231 * r26232;
        float r26234 = a;
        float r26235 = d;
        float r26236 = r26234 * r26235;
        float r26237 = r26233 - r26236;
        float r26238 = r26232 * r26232;
        float r26239 = r26235 * r26235;
        float r26240 = r26238 + r26239;
        float r26241 = r26237 / r26240;
        return r26241;
}

double f_id(double a, double b, double c, double d) {
        double r26242 = b;
        double r26243 = c;
        double r26244 = r26242 * r26243;
        double r26245 = a;
        double r26246 = d;
        double r26247 = r26245 * r26246;
        double r26248 = r26244 - r26247;
        double r26249 = r26243 * r26243;
        double r26250 = r26246 * r26246;
        double r26251 = r26249 + r26250;
        double r26252 = r26248 / r26251;
        return r26252;
}


double f_of(float a, float b, float c, float d) {
        float r26253 = b;
        float r26254 = c;
        float r26255 = r26253 * r26254;
        float r26256 = a;
        float r26257 = d;
        float r26258 = r26256 * r26257;
        float r26259 = r26255 - r26258;
        float r26260 = hypot(r26254, r26257);
        float r26261 = r26259 / r26260;
        float r26262 = log(r26260);
        float r26263 = exp(r26262);
        float r26264 = r26261 / r26263;
        float r26265 = -1.7770805725111145e+308;
        bool r26266 = r26264 <= r26265;
        float r26267 = r26254 / r26257;
        float r26268 = -r26256;
        float r26269 = fma(r26267, r26253, r26268);
        float r26270 = r26269 / r26260;
        float r26271 = 6.20009197993359e+307;
        bool r26272 = r26264 <= r26271;
        float r26273 = r26261 / r26260;
        float r26274 = r26272 ? r26273 : r26270;
        float r26275 = r26266 ? r26270 : r26274;
        return r26275;
}

double f_od(double a, double b, double c, double d) {
        double r26276 = b;
        double r26277 = c;
        double r26278 = r26276 * r26277;
        double r26279 = a;
        double r26280 = d;
        double r26281 = r26279 * r26280;
        double r26282 = r26278 - r26281;
        double r26283 = hypot(r26277, r26280);
        double r26284 = r26282 / r26283;
        double r26285 = log(r26283);
        double r26286 = exp(r26285);
        double r26287 = r26284 / r26286;
        double r26288 = -1.7770805725111145e+308;
        bool r26289 = r26287 <= r26288;
        double r26290 = r26277 / r26280;
        double r26291 = -r26279;
        double r26292 = fma(r26290, r26276, r26291);
        double r26293 = r26292 / r26283;
        double r26294 = 6.20009197993359e+307;
        bool r26295 = r26287 <= r26294;
        double r26296 = r26284 / r26283;
        double r26297 = r26295 ? r26296 : r26293;
        double r26298 = r26289 ? r26293 : r26297;
        return r26298;
}

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 r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init(r26301);
        mpfr_init(r26302);
        mpfr_init(r26303);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r26299, b, MPFR_RNDN);
        mpfr_set_d(r26300, c, MPFR_RNDN);
        mpfr_mul(r26301, r26299, r26300, MPFR_RNDN);
        mpfr_set_d(r26302, a, MPFR_RNDN);
        mpfr_set_d(r26303, d, MPFR_RNDN);
        mpfr_mul(r26304, r26302, r26303, MPFR_RNDN);
        mpfr_sub(r26305, r26301, r26304, MPFR_RNDN);
        mpfr_mul(r26306, r26300, r26300, MPFR_RNDN);
        mpfr_mul(r26307, r26303, r26303, MPFR_RNDN);
        mpfr_add(r26308, r26306, r26307, MPFR_RNDN);
        mpfr_div(r26309, r26305, r26308, MPFR_RNDN);
        return mpfr_get_d(r26309, MPFR_RNDN);
}

static mpfr_t r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26310);
        mpfr_init(r26311);
        mpfr_init(r26312);
        mpfr_init(r26313);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init(r26317);
        mpfr_init(r26318);
        mpfr_init(r26319);
        mpfr_init(r26320);
        mpfr_init(r26321);
        mpfr_init_set_str(r26322, "-1.7770805725111145e+308", 10, MPFR_RNDN);
        mpfr_init(r26323);
        mpfr_init(r26324);
        mpfr_init(r26325);
        mpfr_init(r26326);
        mpfr_init(r26327);
        mpfr_init_set_str(r26328, "6.20009197993359e+307", 10, MPFR_RNDN);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
}

double f_fm(double a, double b, double c, double d) {
        mpfr_set_d(r26310, b, MPFR_RNDN);
        mpfr_set_d(r26311, c, MPFR_RNDN);
        mpfr_mul(r26312, r26310, r26311, MPFR_RNDN);
        mpfr_set_d(r26313, a, MPFR_RNDN);
        mpfr_set_d(r26314, d, MPFR_RNDN);
        mpfr_mul(r26315, r26313, r26314, MPFR_RNDN);
        mpfr_sub(r26316, r26312, r26315, MPFR_RNDN);
        mpfr_hypot(r26317, r26311, r26314, MPFR_RNDN);
        mpfr_div(r26318, r26316, r26317, MPFR_RNDN);
        mpfr_log(r26319, r26317, MPFR_RNDN);
        mpfr_exp(r26320, r26319, MPFR_RNDN);
        mpfr_div(r26321, r26318, r26320, MPFR_RNDN);
        ;
        mpfr_set_si(r26323, mpfr_cmp(r26321, r26322) <= 0, MPFR_RNDN);
        mpfr_div(r26324, r26311, r26314, MPFR_RNDN);
        mpfr_neg(r26325, r26313, MPFR_RNDN);
        mpfr_fma(r26326, r26324, r26310, r26325, MPFR_RNDN);
        mpfr_div(r26327, r26326, r26317, MPFR_RNDN);
        ;
        mpfr_set_si(r26329, mpfr_cmp(r26321, r26328) <= 0, MPFR_RNDN);
        mpfr_div(r26330, r26318, r26317, MPFR_RNDN);
        if (mpfr_get_si(r26329, MPFR_RNDN)) { mpfr_set(r26331, r26330, MPFR_RNDN); } else { mpfr_set(r26331, r26327, MPFR_RNDN); };
        if (mpfr_get_si(r26323, MPFR_RNDN)) { mpfr_set(r26332, r26327, MPFR_RNDN); } else { mpfr_set(r26332, r26331, MPFR_RNDN); };
        return mpfr_get_d(r26332, MPFR_RNDN);
}

static mpfr_t r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351, r26352, r26353, r26354, r26355;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init(r26339);
        mpfr_init(r26340);
        mpfr_init(r26341);
        mpfr_init(r26342);
        mpfr_init(r26343);
        mpfr_init(r26344);
        mpfr_init_set_str(r26345, "-1.7770805725111145e+308", 10, MPFR_RNDN);
        mpfr_init(r26346);
        mpfr_init(r26347);
        mpfr_init(r26348);
        mpfr_init(r26349);
        mpfr_init(r26350);
        mpfr_init_set_str(r26351, "6.20009197993359e+307", 10, MPFR_RNDN);
        mpfr_init(r26352);
        mpfr_init(r26353);
        mpfr_init(r26354);
        mpfr_init(r26355);
}

double f_dm(double a, double b, double c, double d) {
        mpfr_set_d(r26333, b, MPFR_RNDN);
        mpfr_set_d(r26334, c, MPFR_RNDN);
        mpfr_mul(r26335, r26333, r26334, MPFR_RNDN);
        mpfr_set_d(r26336, a, MPFR_RNDN);
        mpfr_set_d(r26337, d, MPFR_RNDN);
        mpfr_mul(r26338, r26336, r26337, MPFR_RNDN);
        mpfr_sub(r26339, r26335, r26338, MPFR_RNDN);
        mpfr_hypot(r26340, r26334, r26337, MPFR_RNDN);
        mpfr_div(r26341, r26339, r26340, MPFR_RNDN);
        mpfr_log(r26342, r26340, MPFR_RNDN);
        mpfr_exp(r26343, r26342, MPFR_RNDN);
        mpfr_div(r26344, r26341, r26343, MPFR_RNDN);
        ;
        mpfr_set_si(r26346, mpfr_cmp(r26344, r26345) <= 0, MPFR_RNDN);
        mpfr_div(r26347, r26334, r26337, MPFR_RNDN);
        mpfr_neg(r26348, r26336, MPFR_RNDN);
        mpfr_fma(r26349, r26347, r26333, r26348, MPFR_RNDN);
        mpfr_div(r26350, r26349, r26340, MPFR_RNDN);
        ;
        mpfr_set_si(r26352, mpfr_cmp(r26344, r26351) <= 0, MPFR_RNDN);
        mpfr_div(r26353, r26341, r26340, MPFR_RNDN);
        if (mpfr_get_si(r26352, MPFR_RNDN)) { mpfr_set(r26354, r26353, MPFR_RNDN); } else { mpfr_set(r26354, r26350, MPFR_RNDN); };
        if (mpfr_get_si(r26346, MPFR_RNDN)) { mpfr_set(r26355, r26350, MPFR_RNDN); } else { mpfr_set(r26355, r26354, MPFR_RNDN); };
        return mpfr_get_d(r26355, MPFR_RNDN);
}

