#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 r18723 = 0.5f;
        float r18724 = 2.0f;
        float r18725 = re;
        float r18726 = r18725 * r18725;
        float r18727 = im;
        float r18728 = r18727 * r18727;
        float r18729 = r18726 + r18728;
        float r18730 = sqrt(r18729);
        float r18731 = r18730 + r18725;
        float r18732 = r18724 * r18731;
        float r18733 = sqrt(r18732);
        float r18734 = r18723 * r18733;
        return r18734;
}

double f_id(double re, double im) {
        double r18735 = 0.5;
        double r18736 = 2.0;
        double r18737 = re;
        double r18738 = r18737 * r18737;
        double r18739 = im;
        double r18740 = r18739 * r18739;
        double r18741 = r18738 + r18740;
        double r18742 = sqrt(r18741);
        double r18743 = r18742 + r18737;
        double r18744 = r18736 * r18743;
        double r18745 = sqrt(r18744);
        double r18746 = r18735 * r18745;
        return r18746;
}


double f_of(float re, float im) {
        float r18747 = 0.5f;
        float r18748 = re;
        float r18749 = im;
        float r18750 = hypot(r18748, r18749);
        float r18751 = 2.0f;
        float r18752 = r18751 * r18748;
        float r18753 = fma(r18750, r18751, r18752);
        float r18754 = sqrt(r18753);
        float r18755 = r18747 * r18754;
        return r18755;
}

double f_od(double re, double im) {
        double r18756 = 0.5;
        double r18757 = re;
        double r18758 = im;
        double r18759 = hypot(r18757, r18758);
        double r18760 = 2.0;
        double r18761 = r18760 * r18757;
        double r18762 = fma(r18759, r18760, r18761);
        double r18763 = sqrt(r18762);
        double r18764 = r18756 * r18763;
        return r18764;
}

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 r18765, r18766, r18767, r18768, r18769, r18770, r18771, r18772, r18773, r18774, r18775, r18776;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18765, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r18766, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18767);
        mpfr_init(r18768);
        mpfr_init(r18769);
        mpfr_init(r18770);
        mpfr_init(r18771);
        mpfr_init(r18772);
        mpfr_init(r18773);
        mpfr_init(r18774);
        mpfr_init(r18775);
        mpfr_init(r18776);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r18767, re, MPFR_RNDN);
        mpfr_mul(r18768, r18767, r18767, MPFR_RNDN);
        mpfr_set_d(r18769, im, MPFR_RNDN);
        mpfr_mul(r18770, r18769, r18769, MPFR_RNDN);
        mpfr_add(r18771, r18768, r18770, MPFR_RNDN);
        mpfr_sqrt(r18772, r18771, MPFR_RNDN);
        mpfr_add(r18773, r18772, r18767, MPFR_RNDN);
        mpfr_mul(r18774, r18766, r18773, MPFR_RNDN);
        mpfr_sqrt(r18775, r18774, MPFR_RNDN);
        mpfr_mul(r18776, r18765, r18775, MPFR_RNDN);
        return mpfr_get_d(r18776, MPFR_RNDN);
}

static mpfr_t r18777, r18778, r18779, r18780, r18781, r18782, r18783, r18784, r18785;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18777, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18778);
        mpfr_init(r18779);
        mpfr_init(r18780);
        mpfr_init_set_str(r18781, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18782);
        mpfr_init(r18783);
        mpfr_init(r18784);
        mpfr_init(r18785);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r18778, re, MPFR_RNDN);
        mpfr_set_d(r18779, im, MPFR_RNDN);
        mpfr_hypot(r18780, r18778, r18779, MPFR_RNDN);
        ;
        mpfr_mul(r18782, r18781, r18778, MPFR_RNDN);
        mpfr_fma(r18783, r18780, r18781, r18782, MPFR_RNDN);
        mpfr_sqrt(r18784, r18783, MPFR_RNDN);
        mpfr_mul(r18785, r18777, r18784, MPFR_RNDN);
        return mpfr_get_d(r18785, MPFR_RNDN);
}

static mpfr_t r18786, r18787, r18788, r18789, r18790, r18791, r18792, r18793, r18794;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18786, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18787);
        mpfr_init(r18788);
        mpfr_init(r18789);
        mpfr_init_set_str(r18790, "2.0", 10, MPFR_RNDN);
        mpfr_init(r18791);
        mpfr_init(r18792);
        mpfr_init(r18793);
        mpfr_init(r18794);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r18787, re, MPFR_RNDN);
        mpfr_set_d(r18788, im, MPFR_RNDN);
        mpfr_hypot(r18789, r18787, r18788, MPFR_RNDN);
        ;
        mpfr_mul(r18791, r18790, r18787, MPFR_RNDN);
        mpfr_fma(r18792, r18789, r18790, r18791, MPFR_RNDN);
        mpfr_sqrt(r18793, r18792, MPFR_RNDN);
        mpfr_mul(r18794, r18786, r18793, MPFR_RNDN);
        return mpfr_get_d(r18794, MPFR_RNDN);
}