#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 r20180 = 0.5;
        float r20181 = 2.0;
        float r20182 = re;
        float r20183 = r20182 * r20182;
        float r20184 = im;
        float r20185 = r20184 * r20184;
        float r20186 = r20183 - r20185;
        float r20187 = sqrt(r20186);
        float r20188 = r20187 + r20182;
        float r20189 = r20181 * r20188;
        float r20190 = sqrt(r20189);
        float r20191 = r20180 * r20190;
        return r20191;
}

double f_id(double re, double im) {
        double r20192 = 0.5;
        double r20193 = 2.0;
        double r20194 = re;
        double r20195 = r20194 * r20194;
        double r20196 = im;
        double r20197 = r20196 * r20196;
        double r20198 = r20195 - r20197;
        double r20199 = sqrt(r20198);
        double r20200 = r20199 + r20194;
        double r20201 = r20193 * r20200;
        double r20202 = sqrt(r20201);
        double r20203 = r20192 * r20202;
        return r20203;
}


double f_of(float re, float im) {
        float r20204 = 0.5;
        float r20205 = 2.0;
        float r20206 = re;
        float r20207 = im;
        float r20208 = r20206 + r20207;
        float r20209 = sqrt(r20208);
        float r20210 = r20206 - r20207;
        float r20211 = sqrt(r20210);
        float r20212 = r20209 * r20211;
        float r20213 = r20212 + r20206;
        float r20214 = r20205 * r20213;
        float r20215 = sqrt(r20214);
        float r20216 = r20204 * r20215;
        return r20216;
}

double f_od(double re, double im) {
        double r20217 = 0.5;
        double r20218 = 2.0;
        double r20219 = re;
        double r20220 = im;
        double r20221 = r20219 + r20220;
        double r20222 = sqrt(r20221);
        double r20223 = r20219 - r20220;
        double r20224 = sqrt(r20223);
        double r20225 = r20222 * r20224;
        double r20226 = r20225 + r20219;
        double r20227 = r20218 * r20226;
        double r20228 = sqrt(r20227);
        double r20229 = r20217 * r20228;
        return r20229;
}

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 r20230, r20231, r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20230, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20231, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20232);
        mpfr_init(r20233);
        mpfr_init(r20234);
        mpfr_init(r20235);
        mpfr_init(r20236);
        mpfr_init(r20237);
        mpfr_init(r20238);
        mpfr_init(r20239);
        mpfr_init(r20240);
        mpfr_init(r20241);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r20232, re, MPFR_RNDN);
        mpfr_mul(r20233, r20232, r20232, MPFR_RNDN);
        mpfr_set_d(r20234, im, MPFR_RNDN);
        mpfr_mul(r20235, r20234, r20234, MPFR_RNDN);
        mpfr_sub(r20236, r20233, r20235, MPFR_RNDN);
        mpfr_sqrt(r20237, r20236, MPFR_RNDN);
        mpfr_add(r20238, r20237, r20232, MPFR_RNDN);
        mpfr_mul(r20239, r20231, r20238, MPFR_RNDN);
        mpfr_sqrt(r20240, r20239, MPFR_RNDN);
        mpfr_mul(r20241, r20230, r20240, MPFR_RNDN);
        return mpfr_get_d(r20241, MPFR_RNDN);
}

static mpfr_t r20242, r20243, r20244, r20245, r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20242, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20243, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20244);
        mpfr_init(r20245);
        mpfr_init(r20246);
        mpfr_init(r20247);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
}

double f_fm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20244, re, MPFR_RNDN);
        mpfr_set_d(r20245, im, MPFR_RNDN);
        mpfr_add(r20246, r20244, r20245, MPFR_RNDN);
        mpfr_sqrt(r20247, r20246, MPFR_RNDN);
        mpfr_sub(r20248, r20244, r20245, MPFR_RNDN);
        mpfr_sqrt(r20249, r20248, MPFR_RNDN);
        mpfr_mul(r20250, r20247, r20249, MPFR_RNDN);
        mpfr_add(r20251, r20250, r20244, MPFR_RNDN);
        mpfr_mul(r20252, r20243, r20251, MPFR_RNDN);
        mpfr_sqrt(r20253, r20252, MPFR_RNDN);
        mpfr_mul(r20254, r20242, r20253, MPFR_RNDN);
        return mpfr_get_d(r20254, MPFR_RNDN);
}

static mpfr_t r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(4496);
        mpfr_init_set_str(r20255, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20256, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20257);
        mpfr_init(r20258);
        mpfr_init(r20259);
        mpfr_init(r20260);
        mpfr_init(r20261);
        mpfr_init(r20262);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init(r20265);
        mpfr_init(r20266);
        mpfr_init(r20267);
}

double f_dm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20257, re, MPFR_RNDN);
        mpfr_set_d(r20258, im, MPFR_RNDN);
        mpfr_add(r20259, r20257, r20258, MPFR_RNDN);
        mpfr_sqrt(r20260, r20259, MPFR_RNDN);
        mpfr_sub(r20261, r20257, r20258, MPFR_RNDN);
        mpfr_sqrt(r20262, r20261, MPFR_RNDN);
        mpfr_mul(r20263, r20260, r20262, MPFR_RNDN);
        mpfr_add(r20264, r20263, r20257, MPFR_RNDN);
        mpfr_mul(r20265, r20256, r20264, MPFR_RNDN);
        mpfr_sqrt(r20266, r20265, MPFR_RNDN);
        mpfr_mul(r20267, r20255, r20266, MPFR_RNDN);
        return mpfr_get_d(r20267, MPFR_RNDN);
}