#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 r27387 = 0.5;
        float r27388 = 2.0;
        float r27389 = re;
        float r27390 = r27389 * r27389;
        float r27391 = im;
        float r27392 = r27391 * r27391;
        float r27393 = r27390 + r27392;
        float r27394 = sqrt(r27393);
        float r27395 = r27394 + r27389;
        float r27396 = r27388 * r27395;
        float r27397 = sqrt(r27396);
        float r27398 = r27387 * r27397;
        return r27398;
}

double f_id(double re, double im) {
        double r27399 = 0.5;
        double r27400 = 2.0;
        double r27401 = re;
        double r27402 = r27401 * r27401;
        double r27403 = im;
        double r27404 = r27403 * r27403;
        double r27405 = r27402 + r27404;
        double r27406 = sqrt(r27405);
        double r27407 = r27406 + r27401;
        double r27408 = r27400 * r27407;
        double r27409 = sqrt(r27408);
        double r27410 = r27399 * r27409;
        return r27410;
}


double f_of(float re, float im) {
        float r27411 = 0.5;
        float r27412 = re;
        float r27413 = im;
        float r27414 = hypot(r27412, r27413);
        float r27415 = 2.0;
        float r27416 = r27412 * r27415;
        float r27417 = fma(r27414, r27415, r27416);
        float r27418 = sqrt(r27417);
        float r27419 = r27411 * r27418;
        return r27419;
}

double f_od(double re, double im) {
        double r27420 = 0.5;
        double r27421 = re;
        double r27422 = im;
        double r27423 = hypot(r27421, r27422);
        double r27424 = 2.0;
        double r27425 = r27421 * r27424;
        double r27426 = fma(r27423, r27424, r27425);
        double r27427 = sqrt(r27426);
        double r27428 = r27420 * r27427;
        return r27428;
}

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 r27429, r27430, r27431, r27432, r27433, r27434, r27435, r27436, r27437, r27438, r27439, r27440;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init_set_str(r27429, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r27430, "2.0", 10, MPFR_RNDN);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
        mpfr_init(r27434);
        mpfr_init(r27435);
        mpfr_init(r27436);
        mpfr_init(r27437);
        mpfr_init(r27438);
        mpfr_init(r27439);
        mpfr_init(r27440);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r27431, re, MPFR_RNDN);
        mpfr_mul(r27432, r27431, r27431, MPFR_RNDN);
        mpfr_set_d(r27433, im, MPFR_RNDN);
        mpfr_mul(r27434, r27433, r27433, MPFR_RNDN);
        mpfr_add(r27435, r27432, r27434, MPFR_RNDN);
        mpfr_sqrt(r27436, r27435, MPFR_RNDN);
        mpfr_add(r27437, r27436, r27431, MPFR_RNDN);
        mpfr_mul(r27438, r27430, r27437, MPFR_RNDN);
        mpfr_sqrt(r27439, r27438, MPFR_RNDN);
        mpfr_mul(r27440, r27429, r27439, MPFR_RNDN);
        return mpfr_get_d(r27440, MPFR_RNDN);
}

static mpfr_t r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init_set_str(r27441, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27442);
        mpfr_init(r27443);
        mpfr_init(r27444);
        mpfr_init_set_str(r27445, "2.0", 10, MPFR_RNDN);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init(r27449);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r27442, re, MPFR_RNDN);
        mpfr_set_d(r27443, im, MPFR_RNDN);
        mpfr_hypot(r27444, r27442, r27443, MPFR_RNDN);
        ;
        mpfr_mul(r27446, r27442, r27445, MPFR_RNDN);
        mpfr_fma(r27447, r27444, r27445, r27446, MPFR_RNDN);
        mpfr_sqrt(r27448, r27447, MPFR_RNDN);
        mpfr_mul(r27449, r27441, r27448, MPFR_RNDN);
        return mpfr_get_d(r27449, MPFR_RNDN);
}

static mpfr_t r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init_set_str(r27450, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27451);
        mpfr_init(r27452);
        mpfr_init(r27453);
        mpfr_init_set_str(r27454, "2.0", 10, MPFR_RNDN);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init(r27457);
        mpfr_init(r27458);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r27451, re, MPFR_RNDN);
        mpfr_set_d(r27452, im, MPFR_RNDN);
        mpfr_hypot(r27453, r27451, r27452, MPFR_RNDN);
        ;
        mpfr_mul(r27455, r27451, r27454, MPFR_RNDN);
        mpfr_fma(r27456, r27453, r27454, r27455, MPFR_RNDN);
        mpfr_sqrt(r27457, r27456, MPFR_RNDN);
        mpfr_mul(r27458, r27450, r27457, MPFR_RNDN);
        return mpfr_get_d(r27458, MPFR_RNDN);
}