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

char *name = "math.sqrt on complex, real part";

double f_if(float re, float im) {
        float r26229 = 0.5;
        float r26230 = 2.0;
        float r26231 = re;
        float r26232 = r26231 * r26231;
        float r26233 = im;
        float r26234 = r26233 * r26233;
        float r26235 = r26232 + r26234;
        float r26236 = sqrt(r26235);
        float r26237 = r26236 + r26231;
        float r26238 = r26230 * r26237;
        float r26239 = sqrt(r26238);
        float r26240 = r26229 * r26239;
        return r26240;
}

double f_id(double re, double im) {
        double r26241 = 0.5;
        double r26242 = 2.0;
        double r26243 = re;
        double r26244 = r26243 * r26243;
        double r26245 = im;
        double r26246 = r26245 * r26245;
        double r26247 = r26244 + r26246;
        double r26248 = sqrt(r26247);
        double r26249 = r26248 + r26243;
        double r26250 = r26242 * r26249;
        double r26251 = sqrt(r26250);
        double r26252 = r26241 * r26251;
        return r26252;
}


double f_of(float re, float im) {
        float r26253 = 2.0;
        float r26254 = im;
        float r26255 = re;
        float r26256 = r26254 + r26255;
        float r26257 = r26253 * r26256;
        float r26258 = -1.847819979812116e-147;
        bool r26259 = r26257 <= r26258;
        float r26260 = 0.5;
        float r26261 = r26254 * r26254;
        float r26262 = r26253 * r26261;
        float r26263 = sqrt(r26262);
        float r26264 = r26255 * r26255;
        float r26265 = r26264 + r26261;
        float r26266 = sqrt(r26265);
        float r26267 = r26266 - r26255;
        float r26268 = sqrt(r26267);
        float r26269 = r26263 / r26268;
        float r26270 = r26260 * r26269;
        float r26271 = 1.842616673093658e+147;
        bool r26272 = r26257 <= r26271;
        float r26273 = r26266 + r26255;
        float r26274 = r26253 * r26273;
        float r26275 = sqrt(r26274);
        float r26276 = r26260 * r26275;
        float r26277 = sqrt(r26257);
        float r26278 = r26260 * r26277;
        float r26279 = r26272 ? r26276 : r26278;
        float r26280 = r26259 ? r26270 : r26279;
        return r26280;
}

double f_od(double re, double im) {
        double r26281 = 2.0;
        double r26282 = im;
        double r26283 = re;
        double r26284 = r26282 + r26283;
        double r26285 = r26281 * r26284;
        double r26286 = -1.847819979812116e-147;
        bool r26287 = r26285 <= r26286;
        double r26288 = 0.5;
        double r26289 = r26282 * r26282;
        double r26290 = r26281 * r26289;
        double r26291 = sqrt(r26290);
        double r26292 = r26283 * r26283;
        double r26293 = r26292 + r26289;
        double r26294 = sqrt(r26293);
        double r26295 = r26294 - r26283;
        double r26296 = sqrt(r26295);
        double r26297 = r26291 / r26296;
        double r26298 = r26288 * r26297;
        double r26299 = 1.842616673093658e+147;
        bool r26300 = r26285 <= r26299;
        double r26301 = r26294 + r26283;
        double r26302 = r26281 * r26301;
        double r26303 = sqrt(r26302);
        double r26304 = r26288 * r26303;
        double r26305 = sqrt(r26285);
        double r26306 = r26288 * r26305;
        double r26307 = r26300 ? r26304 : r26306;
        double r26308 = r26287 ? r26298 : r26307;
        return r26308;
}

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 r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3472);
        mpfr_init_set_str(r26309, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r26310, "2.0", 10, MPFR_RNDN);
        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);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r26311, re, MPFR_RNDN);
        mpfr_mul(r26312, r26311, r26311, MPFR_RNDN);
        mpfr_set_d(r26313, im, MPFR_RNDN);
        mpfr_mul(r26314, r26313, r26313, MPFR_RNDN);
        mpfr_add(r26315, r26312, r26314, MPFR_RNDN);
        mpfr_sqrt(r26316, r26315, MPFR_RNDN);
        mpfr_add(r26317, r26316, r26311, MPFR_RNDN);
        mpfr_mul(r26318, r26310, r26317, MPFR_RNDN);
        mpfr_sqrt(r26319, r26318, MPFR_RNDN);
        mpfr_mul(r26320, r26309, r26319, MPFR_RNDN);
        return mpfr_get_d(r26320, MPFR_RNDN);
}

static mpfr_t r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3472);
        mpfr_init_set_str(r26321, "2.0", 10, MPFR_RNDN);
        mpfr_init(r26322);
        mpfr_init(r26323);
        mpfr_init(r26324);
        mpfr_init(r26325);
        mpfr_init_set_str(r26326, "-1.847819979812116e-147", 10, MPFR_RNDN);
        mpfr_init(r26327);
        mpfr_init_set_str(r26328, "0.5", 10, MPFR_RNDN);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init_set_str(r26339, "1.842616673093658e+147", 10, MPFR_RNDN);
        mpfr_init(r26340);
        mpfr_init(r26341);
        mpfr_init(r26342);
        mpfr_init(r26343);
        mpfr_init(r26344);
        mpfr_init(r26345);
        mpfr_init(r26346);
        mpfr_init(r26347);
        mpfr_init(r26348);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r26322, im, MPFR_RNDN);
        mpfr_set_d(r26323, re, MPFR_RNDN);
        mpfr_add(r26324, r26322, r26323, MPFR_RNDN);
        mpfr_mul(r26325, r26321, r26324, MPFR_RNDN);
        ;
        mpfr_set_si(r26327, mpfr_cmp(r26325, r26326) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r26329, r26322, r26322, MPFR_RNDN);
        mpfr_mul(r26330, r26321, r26329, MPFR_RNDN);
        mpfr_sqrt(r26331, r26330, MPFR_RNDN);
        mpfr_mul(r26332, r26323, r26323, MPFR_RNDN);
        mpfr_add(r26333, r26332, r26329, MPFR_RNDN);
        mpfr_sqrt(r26334, r26333, MPFR_RNDN);
        mpfr_sub(r26335, r26334, r26323, MPFR_RNDN);
        mpfr_sqrt(r26336, r26335, MPFR_RNDN);
        mpfr_div(r26337, r26331, r26336, MPFR_RNDN);
        mpfr_mul(r26338, r26328, r26337, MPFR_RNDN);
        ;
        mpfr_set_si(r26340, mpfr_cmp(r26325, r26339) <= 0, MPFR_RNDN);
        mpfr_add(r26341, r26334, r26323, MPFR_RNDN);
        mpfr_mul(r26342, r26321, r26341, MPFR_RNDN);
        mpfr_sqrt(r26343, r26342, MPFR_RNDN);
        mpfr_mul(r26344, r26328, r26343, MPFR_RNDN);
        mpfr_sqrt(r26345, r26325, MPFR_RNDN);
        mpfr_mul(r26346, r26328, r26345, MPFR_RNDN);
        if (mpfr_get_si(r26340, MPFR_RNDN)) { mpfr_set(r26347, r26344, MPFR_RNDN); } else { mpfr_set(r26347, r26346, MPFR_RNDN); };
        if (mpfr_get_si(r26327, MPFR_RNDN)) { mpfr_set(r26348, r26338, MPFR_RNDN); } else { mpfr_set(r26348, r26347, MPFR_RNDN); };
        return mpfr_get_d(r26348, MPFR_RNDN);
}

static mpfr_t r26349, r26350, r26351, r26352, r26353, r26354, r26355, r26356, r26357, r26358, r26359, r26360, r26361, r26362, r26363, r26364, r26365, r26366, r26367, r26368, r26369, r26370, r26371, r26372, r26373, r26374, r26375, r26376;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3472);
        mpfr_init_set_str(r26349, "2.0", 10, MPFR_RNDN);
        mpfr_init(r26350);
        mpfr_init(r26351);
        mpfr_init(r26352);
        mpfr_init(r26353);
        mpfr_init_set_str(r26354, "-1.847819979812116e-147", 10, MPFR_RNDN);
        mpfr_init(r26355);
        mpfr_init_set_str(r26356, "0.5", 10, MPFR_RNDN);
        mpfr_init(r26357);
        mpfr_init(r26358);
        mpfr_init(r26359);
        mpfr_init(r26360);
        mpfr_init(r26361);
        mpfr_init(r26362);
        mpfr_init(r26363);
        mpfr_init(r26364);
        mpfr_init(r26365);
        mpfr_init(r26366);
        mpfr_init_set_str(r26367, "1.842616673093658e+147", 10, MPFR_RNDN);
        mpfr_init(r26368);
        mpfr_init(r26369);
        mpfr_init(r26370);
        mpfr_init(r26371);
        mpfr_init(r26372);
        mpfr_init(r26373);
        mpfr_init(r26374);
        mpfr_init(r26375);
        mpfr_init(r26376);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r26350, im, MPFR_RNDN);
        mpfr_set_d(r26351, re, MPFR_RNDN);
        mpfr_add(r26352, r26350, r26351, MPFR_RNDN);
        mpfr_mul(r26353, r26349, r26352, MPFR_RNDN);
        ;
        mpfr_set_si(r26355, mpfr_cmp(r26353, r26354) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r26357, r26350, r26350, MPFR_RNDN);
        mpfr_mul(r26358, r26349, r26357, MPFR_RNDN);
        mpfr_sqrt(r26359, r26358, MPFR_RNDN);
        mpfr_mul(r26360, r26351, r26351, MPFR_RNDN);
        mpfr_add(r26361, r26360, r26357, MPFR_RNDN);
        mpfr_sqrt(r26362, r26361, MPFR_RNDN);
        mpfr_sub(r26363, r26362, r26351, MPFR_RNDN);
        mpfr_sqrt(r26364, r26363, MPFR_RNDN);
        mpfr_div(r26365, r26359, r26364, MPFR_RNDN);
        mpfr_mul(r26366, r26356, r26365, MPFR_RNDN);
        ;
        mpfr_set_si(r26368, mpfr_cmp(r26353, r26367) <= 0, MPFR_RNDN);
        mpfr_add(r26369, r26362, r26351, MPFR_RNDN);
        mpfr_mul(r26370, r26349, r26369, MPFR_RNDN);
        mpfr_sqrt(r26371, r26370, MPFR_RNDN);
        mpfr_mul(r26372, r26356, r26371, MPFR_RNDN);
        mpfr_sqrt(r26373, r26353, MPFR_RNDN);
        mpfr_mul(r26374, r26356, r26373, MPFR_RNDN);
        if (mpfr_get_si(r26368, MPFR_RNDN)) { mpfr_set(r26375, r26372, MPFR_RNDN); } else { mpfr_set(r26375, r26374, MPFR_RNDN); };
        if (mpfr_get_si(r26355, MPFR_RNDN)) { mpfr_set(r26376, r26366, MPFR_RNDN); } else { mpfr_set(r26376, r26375, MPFR_RNDN); };
        return mpfr_get_d(r26376, MPFR_RNDN);
}