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

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

double f_if(float a, float b, float c, float d) {
        float r10319 = a;
        float r10320 = c;
        float r10321 = r10319 * r10320;
        float r10322 = b;
        float r10323 = d;
        float r10324 = r10322 * r10323;
        float r10325 = r10321 + r10324;
        float r10326 = r10320 * r10320;
        float r10327 = r10323 * r10323;
        float r10328 = r10326 + r10327;
        float r10329 = r10325 / r10328;
        return r10329;
}

double f_id(double a, double b, double c, double d) {
        double r10330 = a;
        double r10331 = c;
        double r10332 = r10330 * r10331;
        double r10333 = b;
        double r10334 = d;
        double r10335 = r10333 * r10334;
        double r10336 = r10332 + r10335;
        double r10337 = r10331 * r10331;
        double r10338 = r10334 * r10334;
        double r10339 = r10337 + r10338;
        double r10340 = r10336 / r10339;
        return r10340;
}


double f_of(float a, float b, float c, float d) {
        float r10341 = 1;
        float r10342 = c;
        float r10343 = d;
        float r10344 = hypot(r10342, r10343);
        float r10345 = r10341 / r10344;
        float r10346 = b;
        float r10347 = a;
        float r10348 = r10342 * r10347;
        float r10349 = fma(r10346, r10343, r10348);
        float r10350 = r10349 / r10344;
        float r10351 = r10345 * r10350;
        float r10352 = -5.508858057345239e+291;
        bool r10353 = r10351 <= r10352;
        float r10354 = r10347 / r10344;
        float r10355 = 3.768562974861941e+300;
        bool r10356 = r10351 <= r10355;
        float r10357 = r10346 / r10344;
        float r10358 = r10356 ? r10351 : r10357;
        float r10359 = r10353 ? r10354 : r10358;
        return r10359;
}

double f_od(double a, double b, double c, double d) {
        double r10360 = 1;
        double r10361 = c;
        double r10362 = d;
        double r10363 = hypot(r10361, r10362);
        double r10364 = r10360 / r10363;
        double r10365 = b;
        double r10366 = a;
        double r10367 = r10361 * r10366;
        double r10368 = fma(r10365, r10362, r10367);
        double r10369 = r10368 / r10363;
        double r10370 = r10364 * r10369;
        double r10371 = -5.508858057345239e+291;
        bool r10372 = r10370 <= r10371;
        double r10373 = r10366 / r10363;
        double r10374 = 3.768562974861941e+300;
        bool r10375 = r10370 <= r10374;
        double r10376 = r10365 / r10363;
        double r10377 = r10375 ? r10370 : r10376;
        double r10378 = r10372 ? r10373 : r10377;
        return r10378;
}

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 r10379, r10380, r10381, r10382, r10383, r10384, r10385, r10386, r10387, r10388, r10389;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r10379);
        mpfr_init(r10380);
        mpfr_init(r10381);
        mpfr_init(r10382);
        mpfr_init(r10383);
        mpfr_init(r10384);
        mpfr_init(r10385);
        mpfr_init(r10386);
        mpfr_init(r10387);
        mpfr_init(r10388);
        mpfr_init(r10389);
}

double f_im(double a, double b, double c, double d) {
        mpfr_set_d(r10379, a, MPFR_RNDN);
        mpfr_set_d(r10380, c, MPFR_RNDN);
        mpfr_mul(r10381, r10379, r10380, MPFR_RNDN);
        mpfr_set_d(r10382, b, MPFR_RNDN);
        mpfr_set_d(r10383, d, MPFR_RNDN);
        mpfr_mul(r10384, r10382, r10383, MPFR_RNDN);
        mpfr_add(r10385, r10381, r10384, MPFR_RNDN);
        mpfr_mul(r10386, r10380, r10380, MPFR_RNDN);
        mpfr_mul(r10387, r10383, r10383, MPFR_RNDN);
        mpfr_add(r10388, r10386, r10387, MPFR_RNDN);
        mpfr_div(r10389, r10385, r10388, MPFR_RNDN);
        return mpfr_get_d(r10389, MPFR_RNDN);
}

static mpfr_t r10390, r10391, r10392, r10393, r10394, r10395, r10396, r10397, r10398, r10399, r10400, r10401, r10402, r10403, r10404, r10405, r10406, r10407, r10408;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r10390, "1", 10, MPFR_RNDN);
        mpfr_init(r10391);
        mpfr_init(r10392);
        mpfr_init(r10393);
        mpfr_init(r10394);
        mpfr_init(r10395);
        mpfr_init(r10396);
        mpfr_init(r10397);
        mpfr_init(r10398);
        mpfr_init(r10399);
        mpfr_init(r10400);
        mpfr_init_set_str(r10401, "-5.508858057345239e+291", 10, MPFR_RNDN);
        mpfr_init(r10402);
        mpfr_init(r10403);
        mpfr_init_set_str(r10404, "3.768562974861941e+300", 10, MPFR_RNDN);
        mpfr_init(r10405);
        mpfr_init(r10406);
        mpfr_init(r10407);
        mpfr_init(r10408);
}

double f_fm(double a, double b, double c, double d) {
        ;
        mpfr_set_d(r10391, c, MPFR_RNDN);
        mpfr_set_d(r10392, d, MPFR_RNDN);
        mpfr_hypot(r10393, r10391, r10392, MPFR_RNDN);
        mpfr_div(r10394, r10390, r10393, MPFR_RNDN);
        mpfr_set_d(r10395, b, MPFR_RNDN);
        mpfr_set_d(r10396, a, MPFR_RNDN);
        mpfr_mul(r10397, r10391, r10396, MPFR_RNDN);
        mpfr_fma(r10398, r10395, r10392, r10397, MPFR_RNDN);
        mpfr_div(r10399, r10398, r10393, MPFR_RNDN);
        mpfr_mul(r10400, r10394, r10399, MPFR_RNDN);
        ;
        mpfr_set_si(r10402, mpfr_cmp(r10400, r10401) <= 0, MPFR_RNDN);
        mpfr_div(r10403, r10396, r10393, MPFR_RNDN);
        ;
        mpfr_set_si(r10405, mpfr_cmp(r10400, r10404) <= 0, MPFR_RNDN);
        mpfr_div(r10406, r10395, r10393, MPFR_RNDN);
        if (mpfr_get_si(r10405, MPFR_RNDN)) { mpfr_set(r10407, r10400, MPFR_RNDN); } else { mpfr_set(r10407, r10406, MPFR_RNDN); };
        if (mpfr_get_si(r10402, MPFR_RNDN)) { mpfr_set(r10408, r10403, MPFR_RNDN); } else { mpfr_set(r10408, r10407, MPFR_RNDN); };
        return mpfr_get_d(r10408, MPFR_RNDN);
}

static mpfr_t r10409, r10410, r10411, r10412, r10413, r10414, r10415, r10416, r10417, r10418, r10419, r10420, r10421, r10422, r10423, r10424, r10425, r10426, r10427;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r10409, "1", 10, MPFR_RNDN);
        mpfr_init(r10410);
        mpfr_init(r10411);
        mpfr_init(r10412);
        mpfr_init(r10413);
        mpfr_init(r10414);
        mpfr_init(r10415);
        mpfr_init(r10416);
        mpfr_init(r10417);
        mpfr_init(r10418);
        mpfr_init(r10419);
        mpfr_init_set_str(r10420, "-5.508858057345239e+291", 10, MPFR_RNDN);
        mpfr_init(r10421);
        mpfr_init(r10422);
        mpfr_init_set_str(r10423, "3.768562974861941e+300", 10, MPFR_RNDN);
        mpfr_init(r10424);
        mpfr_init(r10425);
        mpfr_init(r10426);
        mpfr_init(r10427);
}

double f_dm(double a, double b, double c, double d) {
        ;
        mpfr_set_d(r10410, c, MPFR_RNDN);
        mpfr_set_d(r10411, d, MPFR_RNDN);
        mpfr_hypot(r10412, r10410, r10411, MPFR_RNDN);
        mpfr_div(r10413, r10409, r10412, MPFR_RNDN);
        mpfr_set_d(r10414, b, MPFR_RNDN);
        mpfr_set_d(r10415, a, MPFR_RNDN);
        mpfr_mul(r10416, r10410, r10415, MPFR_RNDN);
        mpfr_fma(r10417, r10414, r10411, r10416, MPFR_RNDN);
        mpfr_div(r10418, r10417, r10412, MPFR_RNDN);
        mpfr_mul(r10419, r10413, r10418, MPFR_RNDN);
        ;
        mpfr_set_si(r10421, mpfr_cmp(r10419, r10420) <= 0, MPFR_RNDN);
        mpfr_div(r10422, r10415, r10412, MPFR_RNDN);
        ;
        mpfr_set_si(r10424, mpfr_cmp(r10419, r10423) <= 0, MPFR_RNDN);
        mpfr_div(r10425, r10414, r10412, MPFR_RNDN);
        if (mpfr_get_si(r10424, MPFR_RNDN)) { mpfr_set(r10426, r10419, MPFR_RNDN); } else { mpfr_set(r10426, r10425, MPFR_RNDN); };
        if (mpfr_get_si(r10421, MPFR_RNDN)) { mpfr_set(r10427, r10422, MPFR_RNDN); } else { mpfr_set(r10427, r10426, MPFR_RNDN); };
        return mpfr_get_d(r10427, MPFR_RNDN);
}