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

char *name = "_divideComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r22265 = x_im;
        float r22266 = y_re;
        float r22267 = r22265 * r22266;
        float r22268 = x_re;
        float r22269 = y_im;
        float r22270 = r22268 * r22269;
        float r22271 = r22267 - r22270;
        float r22272 = r22266 * r22266;
        float r22273 = r22269 * r22269;
        float r22274 = r22272 + r22273;
        float r22275 = r22271 / r22274;
        return r22275;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r22276 = x_im;
        double r22277 = y_re;
        double r22278 = r22276 * r22277;
        double r22279 = x_re;
        double r22280 = y_im;
        double r22281 = r22279 * r22280;
        double r22282 = r22278 - r22281;
        double r22283 = r22277 * r22277;
        double r22284 = r22280 * r22280;
        double r22285 = r22283 + r22284;
        double r22286 = r22282 / r22285;
        return r22286;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r22287 = y_im;
        float r22288 = -4.426662094485314e-226;
        bool r22289 = r22287 <= r22288;
        float r22290 = y_re;
        float r22291 = hypot(r22290, r22287);
        float r22292 = r22290 / r22291;
        float r22293 = x_im;
        float r22294 = r22291 / r22293;
        float r22295 = r22292 / r22294;
        float r22296 = x_re;
        float r22297 = -r22296;
        float r22298 = r22297 / r22291;
        float r22299 = r22291 / r22287;
        float r22300 = r22298 / r22299;
        float r22301 = r22295 + r22300;
        float r22302 = 9.332805945311764e-40;
        bool r22303 = r22287 <= r22302;
        float r22304 = 1;
        float r22305 = r22304 / r22291;
        float r22306 = r22290 / r22304;
        float r22307 = r22293 / r22291;
        float r22308 = r22306 * r22307;
        float r22309 = r22296 * r22287;
        float r22310 = r22309 / r22291;
        float r22311 = r22308 - r22310;
        float r22312 = r22305 * r22311;
        float r22313 = r22303 ? r22312 : r22301;
        float r22314 = r22289 ? r22301 : r22313;
        return r22314;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r22315 = y_im;
        double r22316 = -4.426662094485314e-226;
        bool r22317 = r22315 <= r22316;
        double r22318 = y_re;
        double r22319 = hypot(r22318, r22315);
        double r22320 = r22318 / r22319;
        double r22321 = x_im;
        double r22322 = r22319 / r22321;
        double r22323 = r22320 / r22322;
        double r22324 = x_re;
        double r22325 = -r22324;
        double r22326 = r22325 / r22319;
        double r22327 = r22319 / r22315;
        double r22328 = r22326 / r22327;
        double r22329 = r22323 + r22328;
        double r22330 = 9.332805945311764e-40;
        bool r22331 = r22315 <= r22330;
        double r22332 = 1;
        double r22333 = r22332 / r22319;
        double r22334 = r22318 / r22332;
        double r22335 = r22321 / r22319;
        double r22336 = r22334 * r22335;
        double r22337 = r22324 * r22315;
        double r22338 = r22337 / r22319;
        double r22339 = r22336 - r22338;
        double r22340 = r22333 * r22339;
        double r22341 = r22331 ? r22340 : r22329;
        double r22342 = r22317 ? r22329 : r22341;
        return r22342;
}

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 r22343, r22344, r22345, r22346, r22347, r22348, r22349, r22350, r22351, r22352, r22353;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r22343);
        mpfr_init(r22344);
        mpfr_init(r22345);
        mpfr_init(r22346);
        mpfr_init(r22347);
        mpfr_init(r22348);
        mpfr_init(r22349);
        mpfr_init(r22350);
        mpfr_init(r22351);
        mpfr_init(r22352);
        mpfr_init(r22353);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22343, x_im, MPFR_RNDN);
        mpfr_set_d(r22344, y_re, MPFR_RNDN);
        mpfr_mul(r22345, r22343, r22344, MPFR_RNDN);
        mpfr_set_d(r22346, x_re, MPFR_RNDN);
        mpfr_set_d(r22347, y_im, MPFR_RNDN);
        mpfr_mul(r22348, r22346, r22347, MPFR_RNDN);
        mpfr_sub(r22349, r22345, r22348, MPFR_RNDN);
        mpfr_mul(r22350, r22344, r22344, MPFR_RNDN);
        mpfr_mul(r22351, r22347, r22347, MPFR_RNDN);
        mpfr_add(r22352, r22350, r22351, MPFR_RNDN);
        mpfr_div(r22353, r22349, r22352, MPFR_RNDN);
        return mpfr_get_d(r22353, MPFR_RNDN);
}

static mpfr_t r22354, r22355, r22356, r22357, r22358, r22359, r22360, r22361, r22362, r22363, r22364, r22365, r22366, r22367, r22368, r22369, r22370, r22371, r22372, r22373, r22374, r22375, r22376, r22377, r22378, r22379, r22380, r22381;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22354);
        mpfr_init_set_str(r22355, "-4.426662094485314e-226", 10, MPFR_RNDN);
        mpfr_init(r22356);
        mpfr_init(r22357);
        mpfr_init(r22358);
        mpfr_init(r22359);
        mpfr_init(r22360);
        mpfr_init(r22361);
        mpfr_init(r22362);
        mpfr_init(r22363);
        mpfr_init(r22364);
        mpfr_init(r22365);
        mpfr_init(r22366);
        mpfr_init(r22367);
        mpfr_init(r22368);
        mpfr_init_set_str(r22369, "9.332805945311764e-40", 10, MPFR_RNDN);
        mpfr_init(r22370);
        mpfr_init_set_str(r22371, "1", 10, MPFR_RNDN);
        mpfr_init(r22372);
        mpfr_init(r22373);
        mpfr_init(r22374);
        mpfr_init(r22375);
        mpfr_init(r22376);
        mpfr_init(r22377);
        mpfr_init(r22378);
        mpfr_init(r22379);
        mpfr_init(r22380);
        mpfr_init(r22381);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22354, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r22356, mpfr_cmp(r22354, r22355) <= 0, MPFR_RNDN);
        mpfr_set_d(r22357, y_re, MPFR_RNDN);
        mpfr_hypot(r22358, r22357, r22354, MPFR_RNDN);
        mpfr_div(r22359, r22357, r22358, MPFR_RNDN);
        mpfr_set_d(r22360, x_im, MPFR_RNDN);
        mpfr_div(r22361, r22358, r22360, MPFR_RNDN);
        mpfr_div(r22362, r22359, r22361, MPFR_RNDN);
        mpfr_set_d(r22363, x_re, MPFR_RNDN);
        mpfr_neg(r22364, r22363, MPFR_RNDN);
        mpfr_div(r22365, r22364, r22358, MPFR_RNDN);
        mpfr_div(r22366, r22358, r22354, MPFR_RNDN);
        mpfr_div(r22367, r22365, r22366, MPFR_RNDN);
        mpfr_add(r22368, r22362, r22367, MPFR_RNDN);
        ;
        mpfr_set_si(r22370, mpfr_cmp(r22354, r22369) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r22372, r22371, r22358, MPFR_RNDN);
        mpfr_div(r22373, r22357, r22371, MPFR_RNDN);
        mpfr_div(r22374, r22360, r22358, MPFR_RNDN);
        mpfr_mul(r22375, r22373, r22374, MPFR_RNDN);
        mpfr_mul(r22376, r22363, r22354, MPFR_RNDN);
        mpfr_div(r22377, r22376, r22358, MPFR_RNDN);
        mpfr_sub(r22378, r22375, r22377, MPFR_RNDN);
        mpfr_mul(r22379, r22372, r22378, MPFR_RNDN);
        if (mpfr_get_si(r22370, MPFR_RNDN)) { mpfr_set(r22380, r22379, MPFR_RNDN); } else { mpfr_set(r22380, r22368, MPFR_RNDN); };
        if (mpfr_get_si(r22356, MPFR_RNDN)) { mpfr_set(r22381, r22368, MPFR_RNDN); } else { mpfr_set(r22381, r22380, MPFR_RNDN); };
        return mpfr_get_d(r22381, MPFR_RNDN);
}

static mpfr_t r22382, r22383, r22384, r22385, r22386, r22387, r22388, r22389, r22390, r22391, r22392, r22393, r22394, r22395, r22396, r22397, r22398, r22399, r22400, r22401, r22402, r22403, r22404, r22405, r22406, r22407, r22408, r22409;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r22382);
        mpfr_init_set_str(r22383, "-4.426662094485314e-226", 10, MPFR_RNDN);
        mpfr_init(r22384);
        mpfr_init(r22385);
        mpfr_init(r22386);
        mpfr_init(r22387);
        mpfr_init(r22388);
        mpfr_init(r22389);
        mpfr_init(r22390);
        mpfr_init(r22391);
        mpfr_init(r22392);
        mpfr_init(r22393);
        mpfr_init(r22394);
        mpfr_init(r22395);
        mpfr_init(r22396);
        mpfr_init_set_str(r22397, "9.332805945311764e-40", 10, MPFR_RNDN);
        mpfr_init(r22398);
        mpfr_init_set_str(r22399, "1", 10, MPFR_RNDN);
        mpfr_init(r22400);
        mpfr_init(r22401);
        mpfr_init(r22402);
        mpfr_init(r22403);
        mpfr_init(r22404);
        mpfr_init(r22405);
        mpfr_init(r22406);
        mpfr_init(r22407);
        mpfr_init(r22408);
        mpfr_init(r22409);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r22382, y_im, MPFR_RNDN);
        ;
        mpfr_set_si(r22384, mpfr_cmp(r22382, r22383) <= 0, MPFR_RNDN);
        mpfr_set_d(r22385, y_re, MPFR_RNDN);
        mpfr_hypot(r22386, r22385, r22382, MPFR_RNDN);
        mpfr_div(r22387, r22385, r22386, MPFR_RNDN);
        mpfr_set_d(r22388, x_im, MPFR_RNDN);
        mpfr_div(r22389, r22386, r22388, MPFR_RNDN);
        mpfr_div(r22390, r22387, r22389, MPFR_RNDN);
        mpfr_set_d(r22391, x_re, MPFR_RNDN);
        mpfr_neg(r22392, r22391, MPFR_RNDN);
        mpfr_div(r22393, r22392, r22386, MPFR_RNDN);
        mpfr_div(r22394, r22386, r22382, MPFR_RNDN);
        mpfr_div(r22395, r22393, r22394, MPFR_RNDN);
        mpfr_add(r22396, r22390, r22395, MPFR_RNDN);
        ;
        mpfr_set_si(r22398, mpfr_cmp(r22382, r22397) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r22400, r22399, r22386, MPFR_RNDN);
        mpfr_div(r22401, r22385, r22399, MPFR_RNDN);
        mpfr_div(r22402, r22388, r22386, MPFR_RNDN);
        mpfr_mul(r22403, r22401, r22402, MPFR_RNDN);
        mpfr_mul(r22404, r22391, r22382, MPFR_RNDN);
        mpfr_div(r22405, r22404, r22386, MPFR_RNDN);
        mpfr_sub(r22406, r22403, r22405, MPFR_RNDN);
        mpfr_mul(r22407, r22400, r22406, MPFR_RNDN);
        if (mpfr_get_si(r22398, MPFR_RNDN)) { mpfr_set(r22408, r22407, MPFR_RNDN); } else { mpfr_set(r22408, r22396, MPFR_RNDN); };
        if (mpfr_get_si(r22384, MPFR_RNDN)) { mpfr_set(r22409, r22396, MPFR_RNDN); } else { mpfr_set(r22409, r22408, MPFR_RNDN); };
        return mpfr_get_d(r22409, MPFR_RNDN);
}