#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 r19144 = 0.5f;
        float r19145 = 2.0f;
        float r19146 = re;
        float r19147 = r19146 * r19146;
        float r19148 = im;
        float r19149 = r19148 * r19148;
        float r19150 = r19147 + r19149;
        float r19151 = sqrt(r19150);
        float r19152 = r19151 + r19146;
        float r19153 = r19145 * r19152;
        float r19154 = sqrt(r19153);
        float r19155 = r19144 * r19154;
        return r19155;
}

double f_id(double re, double im) {
        double r19156 = 0.5;
        double r19157 = 2.0;
        double r19158 = re;
        double r19159 = r19158 * r19158;
        double r19160 = im;
        double r19161 = r19160 * r19160;
        double r19162 = r19159 + r19161;
        double r19163 = sqrt(r19162);
        double r19164 = r19163 + r19158;
        double r19165 = r19157 * r19164;
        double r19166 = sqrt(r19165);
        double r19167 = r19156 * r19166;
        return r19167;
}


double f_of(float re, float im) {
        float r19168 = re;
        float r19169 = -6.6386501330895076e-12f;
        bool r19170 = r19168 <= r19169;
        float r19171 = 0.5f;
        float r19172 = im;
        float r19173 = r19172 * r19172;
        float r19174 = 2.0f;
        float r19175 = r19173 * r19174;
        float r19176 = sqrt(r19175);
        float r19177 = r19171 * r19176;
        float r19178 = -r19168;
        float r19179 = r19178 - r19168;
        float r19180 = sqrt(r19179);
        float r19181 = r19177 / r19180;
        float r19182 = 43888257531904.0f;
        bool r19183 = r19168 <= r19182;
        float r19184 = r19168 * r19168;
        float r19185 = r19184 + r19173;
        float r19186 = sqrt(r19185);
        float r19187 = r19186 + r19168;
        float r19188 = r19174 * r19187;
        float r19189 = sqrt(r19188);
        float r19190 = r19171 * r19189;
        float r19191 = 2.0f;
        float r19192 = r19191 * r19168;
        float r19193 = 0.5f;
        float r19194 = r19172 * r19193;
        float r19195 = r19168 / r19172;
        float r19196 = r19194 / r19195;
        float r19197 = r19192 + r19196;
        float r19198 = r19197 * r19174;
        float r19199 = sqrt(r19198);
        float r19200 = r19171 * r19199;
        float r19201 = r19183 ? r19190 : r19200;
        float r19202 = r19170 ? r19181 : r19201;
        return r19202;
}

double f_od(double re, double im) {
        double r19203 = re;
        double r19204 = -6.6386501330895076e-12;
        bool r19205 = r19203 <= r19204;
        double r19206 = 0.5;
        double r19207 = im;
        double r19208 = r19207 * r19207;
        double r19209 = 2.0;
        double r19210 = r19208 * r19209;
        double r19211 = sqrt(r19210);
        double r19212 = r19206 * r19211;
        double r19213 = -r19203;
        double r19214 = r19213 - r19203;
        double r19215 = sqrt(r19214);
        double r19216 = r19212 / r19215;
        double r19217 = 43888257531904.0;
        bool r19218 = r19203 <= r19217;
        double r19219 = r19203 * r19203;
        double r19220 = r19219 + r19208;
        double r19221 = sqrt(r19220);
        double r19222 = r19221 + r19203;
        double r19223 = r19209 * r19222;
        double r19224 = sqrt(r19223);
        double r19225 = r19206 * r19224;
        double r19226 = 2.0;
        double r19227 = r19226 * r19203;
        double r19228 = 0.5;
        double r19229 = r19207 * r19228;
        double r19230 = r19203 / r19207;
        double r19231 = r19229 / r19230;
        double r19232 = r19227 + r19231;
        double r19233 = r19232 * r19209;
        double r19234 = sqrt(r19233);
        double r19235 = r19206 * r19234;
        double r19236 = r19218 ? r19225 : r19235;
        double r19237 = r19205 ? r19216 : r19236;
        return r19237;
}

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 r19238, r19239, r19240, r19241, r19242, r19243, r19244, r19245, r19246, r19247, r19248, r19249;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19238, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19239, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19240);
        mpfr_init(r19241);
        mpfr_init(r19242);
        mpfr_init(r19243);
        mpfr_init(r19244);
        mpfr_init(r19245);
        mpfr_init(r19246);
        mpfr_init(r19247);
        mpfr_init(r19248);
        mpfr_init(r19249);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r19240, re, MPFR_RNDN);
        mpfr_mul(r19241, r19240, r19240, MPFR_RNDN);
        mpfr_set_d(r19242, im, MPFR_RNDN);
        mpfr_mul(r19243, r19242, r19242, MPFR_RNDN);
        mpfr_add(r19244, r19241, r19243, MPFR_RNDN);
        mpfr_sqrt(r19245, r19244, MPFR_RNDN);
        mpfr_add(r19246, r19245, r19240, MPFR_RNDN);
        mpfr_mul(r19247, r19239, r19246, MPFR_RNDN);
        mpfr_sqrt(r19248, r19247, MPFR_RNDN);
        mpfr_mul(r19249, r19238, r19248, MPFR_RNDN);
        return mpfr_get_d(r19249, MPFR_RNDN);
}

static mpfr_t r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265, r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19250);
        mpfr_init_set_str(r19251, "-6.63865f-12", 10, MPFR_RNDN);
        mpfr_init(r19252);
        mpfr_init_set_str(r19253, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19254);
        mpfr_init(r19255);
        mpfr_init_set_str(r19256, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19257);
        mpfr_init(r19258);
        mpfr_init(r19259);
        mpfr_init(r19260);
        mpfr_init(r19261);
        mpfr_init(r19262);
        mpfr_init(r19263);
        mpfr_init_set_str(r19264, "4.3888258f+13", 10, MPFR_RNDN);
        mpfr_init(r19265);
        mpfr_init(r19266);
        mpfr_init(r19267);
        mpfr_init(r19268);
        mpfr_init(r19269);
        mpfr_init(r19270);
        mpfr_init(r19271);
        mpfr_init(r19272);
        mpfr_init_set_str(r19273, "2", 10, MPFR_RNDN);
        mpfr_init(r19274);
        mpfr_init_set_str(r19275, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19276);
        mpfr_init(r19277);
        mpfr_init(r19278);
        mpfr_init(r19279);
        mpfr_init(r19280);
        mpfr_init(r19281);
        mpfr_init(r19282);
        mpfr_init(r19283);
        mpfr_init(r19284);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19250, re, MPFR_RNDN);
        ;
        mpfr_set_si(r19252, mpfr_cmp(r19250, r19251) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19254, im, MPFR_RNDN);
        mpfr_mul(r19255, r19254, r19254, MPFR_RNDN);
        ;
        mpfr_mul(r19257, r19255, r19256, MPFR_RNDN);
        mpfr_sqrt(r19258, r19257, MPFR_RNDN);
        mpfr_mul(r19259, r19253, r19258, MPFR_RNDN);
        mpfr_neg(r19260, r19250, MPFR_RNDN);
        mpfr_sub(r19261, r19260, r19250, MPFR_RNDN);
        mpfr_sqrt(r19262, r19261, MPFR_RNDN);
        mpfr_div(r19263, r19259, r19262, MPFR_RNDN);
        ;
        mpfr_set_si(r19265, mpfr_cmp(r19250, r19264) <= 0, MPFR_RNDN);
        mpfr_mul(r19266, r19250, r19250, MPFR_RNDN);
        mpfr_add(r19267, r19266, r19255, MPFR_RNDN);
        mpfr_sqrt(r19268, r19267, MPFR_RNDN);
        mpfr_add(r19269, r19268, r19250, MPFR_RNDN);
        mpfr_mul(r19270, r19256, r19269, MPFR_RNDN);
        mpfr_sqrt(r19271, r19270, MPFR_RNDN);
        mpfr_mul(r19272, r19253, r19271, MPFR_RNDN);
        ;
        mpfr_mul(r19274, r19273, r19250, MPFR_RNDN);
        ;
        mpfr_mul(r19276, r19254, r19275, MPFR_RNDN);
        mpfr_div(r19277, r19250, r19254, MPFR_RNDN);
        mpfr_div(r19278, r19276, r19277, MPFR_RNDN);
        mpfr_add(r19279, r19274, r19278, MPFR_RNDN);
        mpfr_mul(r19280, r19279, r19256, MPFR_RNDN);
        mpfr_sqrt(r19281, r19280, MPFR_RNDN);
        mpfr_mul(r19282, r19253, r19281, MPFR_RNDN);
        if (mpfr_get_si(r19265, MPFR_RNDN)) { mpfr_set(r19283, r19272, MPFR_RNDN); } else { mpfr_set(r19283, r19282, MPFR_RNDN); };
        if (mpfr_get_si(r19252, MPFR_RNDN)) { mpfr_set(r19284, r19263, MPFR_RNDN); } else { mpfr_set(r19284, r19283, MPFR_RNDN); };
        return mpfr_get_d(r19284, MPFR_RNDN);
}

static mpfr_t r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301, r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317, r19318, r19319;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19285);
        mpfr_init_set_str(r19286, "-6.63865f-12", 10, MPFR_RNDN);
        mpfr_init(r19287);
        mpfr_init_set_str(r19288, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19289);
        mpfr_init(r19290);
        mpfr_init_set_str(r19291, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19292);
        mpfr_init(r19293);
        mpfr_init(r19294);
        mpfr_init(r19295);
        mpfr_init(r19296);
        mpfr_init(r19297);
        mpfr_init(r19298);
        mpfr_init_set_str(r19299, "4.3888258f+13", 10, MPFR_RNDN);
        mpfr_init(r19300);
        mpfr_init(r19301);
        mpfr_init(r19302);
        mpfr_init(r19303);
        mpfr_init(r19304);
        mpfr_init(r19305);
        mpfr_init(r19306);
        mpfr_init(r19307);
        mpfr_init_set_str(r19308, "2", 10, MPFR_RNDN);
        mpfr_init(r19309);
        mpfr_init_set_str(r19310, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19311);
        mpfr_init(r19312);
        mpfr_init(r19313);
        mpfr_init(r19314);
        mpfr_init(r19315);
        mpfr_init(r19316);
        mpfr_init(r19317);
        mpfr_init(r19318);
        mpfr_init(r19319);
}

double f_dm(double re, double im) {
        mpfr_set_d(r19285, re, MPFR_RNDN);
        ;
        mpfr_set_si(r19287, mpfr_cmp(r19285, r19286) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19289, im, MPFR_RNDN);
        mpfr_mul(r19290, r19289, r19289, MPFR_RNDN);
        ;
        mpfr_mul(r19292, r19290, r19291, MPFR_RNDN);
        mpfr_sqrt(r19293, r19292, MPFR_RNDN);
        mpfr_mul(r19294, r19288, r19293, MPFR_RNDN);
        mpfr_neg(r19295, r19285, MPFR_RNDN);
        mpfr_sub(r19296, r19295, r19285, MPFR_RNDN);
        mpfr_sqrt(r19297, r19296, MPFR_RNDN);
        mpfr_div(r19298, r19294, r19297, MPFR_RNDN);
        ;
        mpfr_set_si(r19300, mpfr_cmp(r19285, r19299) <= 0, MPFR_RNDN);
        mpfr_mul(r19301, r19285, r19285, MPFR_RNDN);
        mpfr_add(r19302, r19301, r19290, MPFR_RNDN);
        mpfr_sqrt(r19303, r19302, MPFR_RNDN);
        mpfr_add(r19304, r19303, r19285, MPFR_RNDN);
        mpfr_mul(r19305, r19291, r19304, MPFR_RNDN);
        mpfr_sqrt(r19306, r19305, MPFR_RNDN);
        mpfr_mul(r19307, r19288, r19306, MPFR_RNDN);
        ;
        mpfr_mul(r19309, r19308, r19285, MPFR_RNDN);
        ;
        mpfr_mul(r19311, r19289, r19310, MPFR_RNDN);
        mpfr_div(r19312, r19285, r19289, MPFR_RNDN);
        mpfr_div(r19313, r19311, r19312, MPFR_RNDN);
        mpfr_add(r19314, r19309, r19313, MPFR_RNDN);
        mpfr_mul(r19315, r19314, r19291, MPFR_RNDN);
        mpfr_sqrt(r19316, r19315, MPFR_RNDN);
        mpfr_mul(r19317, r19288, r19316, MPFR_RNDN);
        if (mpfr_get_si(r19300, MPFR_RNDN)) { mpfr_set(r19318, r19307, MPFR_RNDN); } else { mpfr_set(r19318, r19317, MPFR_RNDN); };
        if (mpfr_get_si(r19287, MPFR_RNDN)) { mpfr_set(r19319, r19298, MPFR_RNDN); } else { mpfr_set(r19319, r19318, MPFR_RNDN); };
        return mpfr_get_d(r19319, MPFR_RNDN);
}