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

char *name = "math.sqrt on complex, imaginary part, im greater than 0 branch";

double f_if(float re, float im) {
        float r20196 = 0.5;
        float r20197 = 2.0;
        float r20198 = re;
        float r20199 = r20198 * r20198;
        float r20200 = im;
        float r20201 = r20200 * r20200;
        float r20202 = r20199 - r20201;
        float r20203 = sqrt(r20202);
        float r20204 = r20203 + r20198;
        float r20205 = r20197 * r20204;
        float r20206 = sqrt(r20205);
        float r20207 = r20196 * r20206;
        return r20207;
}

double f_id(double re, double im) {
        double r20208 = 0.5;
        double r20209 = 2.0;
        double r20210 = re;
        double r20211 = r20210 * r20210;
        double r20212 = im;
        double r20213 = r20212 * r20212;
        double r20214 = r20211 - r20213;
        double r20215 = sqrt(r20214);
        double r20216 = r20215 + r20210;
        double r20217 = r20209 * r20216;
        double r20218 = sqrt(r20217);
        double r20219 = r20208 * r20218;
        return r20219;
}


double f_of(float re, float im) {
        float r20220 = 0.5;
        float r20221 = 2.0;
        float r20222 = re;
        float r20223 = im;
        float r20224 = r20222 + r20223;
        float r20225 = sqrt(r20224);
        float r20226 = r20222 - r20223;
        float r20227 = sqrt(r20226);
        float r20228 = r20225 * r20227;
        float r20229 = r20228 + r20222;
        float r20230 = r20221 * r20229;
        float r20231 = sqrt(r20230);
        float r20232 = r20220 * r20231;
        return r20232;
}

double f_od(double re, double im) {
        double r20233 = 0.5;
        double r20234 = 2.0;
        double r20235 = re;
        double r20236 = im;
        double r20237 = r20235 + r20236;
        double r20238 = sqrt(r20237);
        double r20239 = r20235 - r20236;
        double r20240 = sqrt(r20239);
        double r20241 = r20238 * r20240;
        double r20242 = r20241 + r20235;
        double r20243 = r20234 * r20242;
        double r20244 = sqrt(r20243);
        double r20245 = r20233 * r20244;
        return r20245;
}

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 r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20246, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20247, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r20248, re, MPFR_RNDN);
        mpfr_mul(r20249, r20248, r20248, MPFR_RNDN);
        mpfr_set_d(r20250, im, MPFR_RNDN);
        mpfr_mul(r20251, r20250, r20250, MPFR_RNDN);
        mpfr_sub(r20252, r20249, r20251, MPFR_RNDN);
        mpfr_sqrt(r20253, r20252, MPFR_RNDN);
        mpfr_add(r20254, r20253, r20248, MPFR_RNDN);
        mpfr_mul(r20255, r20247, r20254, MPFR_RNDN);
        mpfr_sqrt(r20256, r20255, MPFR_RNDN);
        mpfr_mul(r20257, r20246, r20256, MPFR_RNDN);
        return mpfr_get_d(r20257, MPFR_RNDN);
}

static mpfr_t r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267, r20268, r20269, r20270;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20258, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20259, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20260);
        mpfr_init(r20261);
        mpfr_init(r20262);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init(r20265);
        mpfr_init(r20266);
        mpfr_init(r20267);
        mpfr_init(r20268);
        mpfr_init(r20269);
        mpfr_init(r20270);
}

double f_fm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20260, re, MPFR_RNDN);
        mpfr_set_d(r20261, im, MPFR_RNDN);
        mpfr_add(r20262, r20260, r20261, MPFR_RNDN);
        mpfr_sqrt(r20263, r20262, MPFR_RNDN);
        mpfr_sub(r20264, r20260, r20261, MPFR_RNDN);
        mpfr_sqrt(r20265, r20264, MPFR_RNDN);
        mpfr_mul(r20266, r20263, r20265, MPFR_RNDN);
        mpfr_add(r20267, r20266, r20260, MPFR_RNDN);
        mpfr_mul(r20268, r20259, r20267, MPFR_RNDN);
        mpfr_sqrt(r20269, r20268, MPFR_RNDN);
        mpfr_mul(r20270, r20258, r20269, MPFR_RNDN);
        return mpfr_get_d(r20270, MPFR_RNDN);
}

static mpfr_t r20271, r20272, r20273, r20274, r20275, r20276, r20277, r20278, r20279, r20280, r20281, r20282, r20283;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20271, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20272, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20273);
        mpfr_init(r20274);
        mpfr_init(r20275);
        mpfr_init(r20276);
        mpfr_init(r20277);
        mpfr_init(r20278);
        mpfr_init(r20279);
        mpfr_init(r20280);
        mpfr_init(r20281);
        mpfr_init(r20282);
        mpfr_init(r20283);
}

double f_dm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20273, re, MPFR_RNDN);
        mpfr_set_d(r20274, im, MPFR_RNDN);
        mpfr_add(r20275, r20273, r20274, MPFR_RNDN);
        mpfr_sqrt(r20276, r20275, MPFR_RNDN);
        mpfr_sub(r20277, r20273, r20274, MPFR_RNDN);
        mpfr_sqrt(r20278, r20277, MPFR_RNDN);
        mpfr_mul(r20279, r20276, r20278, MPFR_RNDN);
        mpfr_add(r20280, r20279, r20273, MPFR_RNDN);
        mpfr_mul(r20281, r20272, r20280, MPFR_RNDN);
        mpfr_sqrt(r20282, r20281, MPFR_RNDN);
        mpfr_mul(r20283, r20271, r20282, MPFR_RNDN);
        return mpfr_get_d(r20283, MPFR_RNDN);
}