#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 r18427 = 0.5f;
        float r18428 = 2.0f;
        float r18429 = re;
        float r18430 = r18429 * r18429;
        float r18431 = im;
        float r18432 = r18431 * r18431;
        float r18433 = r18430 + r18432;
        float r18434 = sqrt(r18433);
        float r18435 = r18434 + r18429;
        float r18436 = r18428 * r18435;
        float r18437 = sqrt(r18436);
        float r18438 = r18427 * r18437;
        return r18438;
}

double f_id(double re, double im) {
        double r18439 = 0.5;
        double r18440 = 2.0;
        double r18441 = re;
        double r18442 = r18441 * r18441;
        double r18443 = im;
        double r18444 = r18443 * r18443;
        double r18445 = r18442 + r18444;
        double r18446 = sqrt(r18445);
        double r18447 = r18446 + r18441;
        double r18448 = r18440 * r18447;
        double r18449 = sqrt(r18448);
        double r18450 = r18439 * r18449;
        return r18450;
}


double f_of(float re, float im) {
        float r18451 = 0.5f;
        float r18452 = re;
        float r18453 = im;
        float r18454 = hypot(r18452, r18453);
        float r18455 = 2.0f;
        float r18456 = r18455 * r18452;
        float r18457 = fma(r18454, r18455, r18456);
        float r18458 = sqrt(r18457);
        float r18459 = r18451 * r18458;
        return r18459;
}

double f_od(double re, double im) {
        double r18460 = 0.5;
        double r18461 = re;
        double r18462 = im;
        double r18463 = hypot(r18461, r18462);
        double r18464 = 2.0;
        double r18465 = r18464 * r18461;
        double r18466 = fma(r18463, r18464, r18465);
        double r18467 = sqrt(r18466);
        double r18468 = r18460 * r18467;
        return r18468;
}

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 r18469, r18470, r18471, r18472, r18473, r18474, r18475, r18476, r18477, r18478, r18479, r18480;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18469, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r18470, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18471);
        mpfr_init(r18472);
        mpfr_init(r18473);
        mpfr_init(r18474);
        mpfr_init(r18475);
        mpfr_init(r18476);
        mpfr_init(r18477);
        mpfr_init(r18478);
        mpfr_init(r18479);
        mpfr_init(r18480);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r18471, re, MPFR_RNDN);
        mpfr_mul(r18472, r18471, r18471, MPFR_RNDN);
        mpfr_set_d(r18473, im, MPFR_RNDN);
        mpfr_mul(r18474, r18473, r18473, MPFR_RNDN);
        mpfr_add(r18475, r18472, r18474, MPFR_RNDN);
        mpfr_sqrt(r18476, r18475, MPFR_RNDN);
        mpfr_add(r18477, r18476, r18471, MPFR_RNDN);
        mpfr_mul(r18478, r18470, r18477, MPFR_RNDN);
        mpfr_sqrt(r18479, r18478, MPFR_RNDN);
        mpfr_mul(r18480, r18469, r18479, MPFR_RNDN);
        return mpfr_get_d(r18480, MPFR_RNDN);
}

static mpfr_t r18481, r18482, r18483, r18484, r18485, r18486, r18487, r18488, r18489;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18481, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18482);
        mpfr_init(r18483);
        mpfr_init(r18484);
        mpfr_init_set_str(r18485, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18486);
        mpfr_init(r18487);
        mpfr_init(r18488);
        mpfr_init(r18489);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r18482, re, MPFR_RNDN);
        mpfr_set_d(r18483, im, MPFR_RNDN);
        mpfr_hypot(r18484, r18482, r18483, MPFR_RNDN);
        ;
        mpfr_mul(r18486, r18485, r18482, MPFR_RNDN);
        mpfr_fma(r18487, r18484, r18485, r18486, MPFR_RNDN);
        mpfr_sqrt(r18488, r18487, MPFR_RNDN);
        mpfr_mul(r18489, r18481, r18488, MPFR_RNDN);
        return mpfr_get_d(r18489, MPFR_RNDN);
}

static mpfr_t r18490, r18491, r18492, r18493, r18494, r18495, r18496, r18497, r18498;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18490, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18491);
        mpfr_init(r18492);
        mpfr_init(r18493);
        mpfr_init_set_str(r18494, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18495);
        mpfr_init(r18496);
        mpfr_init(r18497);
        mpfr_init(r18498);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r18491, re, MPFR_RNDN);
        mpfr_set_d(r18492, im, MPFR_RNDN);
        mpfr_hypot(r18493, r18491, r18492, MPFR_RNDN);
        ;
        mpfr_mul(r18495, r18494, r18491, MPFR_RNDN);
        mpfr_fma(r18496, r18493, r18494, r18495, MPFR_RNDN);
        mpfr_sqrt(r18497, r18496, MPFR_RNDN);
        mpfr_mul(r18498, r18490, r18497, MPFR_RNDN);
        return mpfr_get_d(r18498, MPFR_RNDN);
}