#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 r18463 = 0.5f;
        float r18464 = 2.0f;
        float r18465 = re;
        float r18466 = r18465 * r18465;
        float r18467 = im;
        float r18468 = r18467 * r18467;
        float r18469 = r18466 + r18468;
        float r18470 = sqrt(r18469);
        float r18471 = r18470 + r18465;
        float r18472 = r18464 * r18471;
        float r18473 = sqrt(r18472);
        float r18474 = r18463 * r18473;
        return r18474;
}

double f_id(double re, double im) {
        double r18475 = 0.5;
        double r18476 = 2.0;
        double r18477 = re;
        double r18478 = r18477 * r18477;
        double r18479 = im;
        double r18480 = r18479 * r18479;
        double r18481 = r18478 + r18480;
        double r18482 = sqrt(r18481);
        double r18483 = r18482 + r18477;
        double r18484 = r18476 * r18483;
        double r18485 = sqrt(r18484);
        double r18486 = r18475 * r18485;
        return r18486;
}


double f_of(float re, float im) {
        float r18487 = 0.5f;
        float r18488 = re;
        float r18489 = im;
        float r18490 = hypot(r18488, r18489);
        float r18491 = 2.0f;
        float r18492 = r18491 * r18488;
        float r18493 = fma(r18490, r18491, r18492);
        float r18494 = sqrt(r18493);
        float r18495 = r18487 * r18494;
        return r18495;
}

double f_od(double re, double im) {
        double r18496 = 0.5;
        double r18497 = re;
        double r18498 = im;
        double r18499 = hypot(r18497, r18498);
        double r18500 = 2.0;
        double r18501 = r18500 * r18497;
        double r18502 = fma(r18499, r18500, r18501);
        double r18503 = sqrt(r18502);
        double r18504 = r18496 * r18503;
        return r18504;
}

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 r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513, r18514, r18515, r18516;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18505, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r18506, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18507);
        mpfr_init(r18508);
        mpfr_init(r18509);
        mpfr_init(r18510);
        mpfr_init(r18511);
        mpfr_init(r18512);
        mpfr_init(r18513);
        mpfr_init(r18514);
        mpfr_init(r18515);
        mpfr_init(r18516);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r18507, re, MPFR_RNDN);
        mpfr_mul(r18508, r18507, r18507, MPFR_RNDN);
        mpfr_set_d(r18509, im, MPFR_RNDN);
        mpfr_mul(r18510, r18509, r18509, MPFR_RNDN);
        mpfr_add(r18511, r18508, r18510, MPFR_RNDN);
        mpfr_sqrt(r18512, r18511, MPFR_RNDN);
        mpfr_add(r18513, r18512, r18507, MPFR_RNDN);
        mpfr_mul(r18514, r18506, r18513, MPFR_RNDN);
        mpfr_sqrt(r18515, r18514, MPFR_RNDN);
        mpfr_mul(r18516, r18505, r18515, MPFR_RNDN);
        return mpfr_get_d(r18516, MPFR_RNDN);
}

static mpfr_t r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18517, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18518);
        mpfr_init(r18519);
        mpfr_init(r18520);
        mpfr_init_set_str(r18521, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18522);
        mpfr_init(r18523);
        mpfr_init(r18524);
        mpfr_init(r18525);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r18518, re, MPFR_RNDN);
        mpfr_set_d(r18519, im, MPFR_RNDN);
        mpfr_hypot(r18520, r18518, r18519, MPFR_RNDN);
        ;
        mpfr_mul(r18522, r18521, r18518, MPFR_RNDN);
        mpfr_fma(r18523, r18520, r18521, r18522, MPFR_RNDN);
        mpfr_sqrt(r18524, r18523, MPFR_RNDN);
        mpfr_mul(r18525, r18517, r18524, MPFR_RNDN);
        return mpfr_get_d(r18525, MPFR_RNDN);
}

static mpfr_t r18526, r18527, r18528, r18529, r18530, r18531, r18532, r18533, r18534;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18526, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18527);
        mpfr_init(r18528);
        mpfr_init(r18529);
        mpfr_init_set_str(r18530, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18531);
        mpfr_init(r18532);
        mpfr_init(r18533);
        mpfr_init(r18534);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r18527, re, MPFR_RNDN);
        mpfr_set_d(r18528, im, MPFR_RNDN);
        mpfr_hypot(r18529, r18527, r18528, MPFR_RNDN);
        ;
        mpfr_mul(r18531, r18530, r18527, MPFR_RNDN);
        mpfr_fma(r18532, r18529, r18530, r18531, MPFR_RNDN);
        mpfr_sqrt(r18533, r18532, MPFR_RNDN);
        mpfr_mul(r18534, r18526, r18533, MPFR_RNDN);
        return mpfr_get_d(r18534, MPFR_RNDN);
}