#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 r19579 = 0.5f;
        float r19580 = 2.0f;
        float r19581 = re;
        float r19582 = r19581 * r19581;
        float r19583 = im;
        float r19584 = r19583 * r19583;
        float r19585 = r19582 + r19584;
        float r19586 = sqrt(r19585);
        float r19587 = r19586 + r19581;
        float r19588 = r19580 * r19587;
        float r19589 = sqrt(r19588);
        float r19590 = r19579 * r19589;
        return r19590;
}

double f_id(double re, double im) {
        double r19591 = 0.5;
        double r19592 = 2.0;
        double r19593 = re;
        double r19594 = r19593 * r19593;
        double r19595 = im;
        double r19596 = r19595 * r19595;
        double r19597 = r19594 + r19596;
        double r19598 = sqrt(r19597);
        double r19599 = r19598 + r19593;
        double r19600 = r19592 * r19599;
        double r19601 = sqrt(r19600);
        double r19602 = r19591 * r19601;
        return r19602;
}


double f_of(float re, float im) {
        float r19603 = re;
        float r19604 = -1.5608708420993077e-235f;
        bool r19605 = r19603 <= r19604;
        float r19606 = 0.5f;
        float r19607 = 2.0f;
        float r19608 = im;
        float r19609 = r19607 * r19608;
        float r19610 = r19609 * r19608;
        float r19611 = sqrt(r19610);
        float r19612 = r19603 * r19603;
        float r19613 = r19608 * r19608;
        float r19614 = r19612 + r19613;
        float r19615 = sqrt(r19614);
        float r19616 = r19615 - r19603;
        float r19617 = sqrt(r19616);
        float r19618 = r19611 / r19617;
        float r19619 = r19606 * r19618;
        float r19620 = -2.3665634646026902e-293f;
        bool r19621 = r19603 <= r19620;
        float r19622 = r19608 + r19603;
        float r19623 = r19607 * r19622;
        float r19624 = sqrt(r19623);
        float r19625 = r19606 * r19624;
        float r19626 = 1.245134237263001e-160f;
        bool r19627 = r19603 <= r19626;
        float r19628 = r19613 / r19616;
        float r19629 = r19607 * r19628;
        float r19630 = sqrt(r19629);
        float r19631 = r19606 * r19630;
        float r19632 = 2.56976209479494e+109f;
        bool r19633 = r19603 <= r19632;
        float r19634 = r19603 * r19603;
        float r19635 = r19634 + r19613;
        float r19636 = sqrt(r19635);
        float r19637 = r19636 * (r19636 * r19636);
        float r19638 = cbrt(r19637);
        float r19639 = r19638 + r19603;
        float r19640 = r19607 * r19639;
        float r19641 = sqrt(r19640);
        float r19642 = r19606 * r19641;
        float r19643 = r19603 + r19603;
        float r19644 = r19607 * r19643;
        float r19645 = sqrt(r19644);
        float r19646 = r19606 * r19645;
        float r19647 = r19633 ? r19642 : r19646;
        float r19648 = r19627 ? r19631 : r19647;
        float r19649 = r19621 ? r19625 : r19648;
        float r19650 = r19605 ? r19619 : r19649;
        return r19650;
}

double f_od(double re, double im) {
        double r19651 = re;
        double r19652 = -1.5608708420993077e-235;
        bool r19653 = r19651 <= r19652;
        double r19654 = 0.5;
        double r19655 = 2.0;
        double r19656 = im;
        double r19657 = r19655 * r19656;
        double r19658 = r19657 * r19656;
        double r19659 = sqrt(r19658);
        double r19660 = r19651 * r19651;
        double r19661 = r19656 * r19656;
        double r19662 = r19660 + r19661;
        double r19663 = sqrt(r19662);
        double r19664 = r19663 - r19651;
        double r19665 = sqrt(r19664);
        double r19666 = r19659 / r19665;
        double r19667 = r19654 * r19666;
        double r19668 = -2.3665634646026902e-293;
        bool r19669 = r19651 <= r19668;
        double r19670 = r19656 + r19651;
        double r19671 = r19655 * r19670;
        double r19672 = sqrt(r19671);
        double r19673 = r19654 * r19672;
        double r19674 = 1.245134237263001e-160;
        bool r19675 = r19651 <= r19674;
        double r19676 = r19661 / r19664;
        double r19677 = r19655 * r19676;
        double r19678 = sqrt(r19677);
        double r19679 = r19654 * r19678;
        double r19680 = 2.56976209479494e+109;
        bool r19681 = r19651 <= r19680;
        double r19682 = r19651 * r19651;
        double r19683 = r19682 + r19661;
        double r19684 = sqrt(r19683);
        double r19685 = r19684 * (r19684 * r19684);
        double r19686 = cbrt(r19685);
        double r19687 = r19686 + r19651;
        double r19688 = r19655 * r19687;
        double r19689 = sqrt(r19688);
        double r19690 = r19654 * r19689;
        double r19691 = r19651 + r19651;
        double r19692 = r19655 * r19691;
        double r19693 = sqrt(r19692);
        double r19694 = r19654 * r19693;
        double r19695 = r19681 ? r19690 : r19694;
        double r19696 = r19675 ? r19679 : r19695;
        double r19697 = r19669 ? r19673 : r19696;
        double r19698 = r19653 ? r19667 : r19697;
        return r19698;
}

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 r19699, r19700, r19701, r19702, r19703, r19704, r19705, r19706, r19707, r19708, r19709, r19710;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19699, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19700, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19701);
        mpfr_init(r19702);
        mpfr_init(r19703);
        mpfr_init(r19704);
        mpfr_init(r19705);
        mpfr_init(r19706);
        mpfr_init(r19707);
        mpfr_init(r19708);
        mpfr_init(r19709);
        mpfr_init(r19710);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r19701, re, MPFR_RNDN);
        mpfr_mul(r19702, r19701, r19701, MPFR_RNDN);
        mpfr_set_d(r19703, im, MPFR_RNDN);
        mpfr_mul(r19704, r19703, r19703, MPFR_RNDN);
        mpfr_add(r19705, r19702, r19704, MPFR_RNDN);
        mpfr_sqrt(r19706, r19705, MPFR_RNDN);
        mpfr_add(r19707, r19706, r19701, MPFR_RNDN);
        mpfr_mul(r19708, r19700, r19707, MPFR_RNDN);
        mpfr_sqrt(r19709, r19708, MPFR_RNDN);
        mpfr_mul(r19710, r19699, r19709, MPFR_RNDN);
        return mpfr_get_d(r19710, MPFR_RNDN);
}

static mpfr_t r19711, r19712, r19713, r19714, r19715, r19716, r19717, r19718, r19719, r19720, r19721, r19722, r19723, r19724, r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19711);
        mpfr_init_set_str(r19712, "-1.5608708420993077e-235", 10, MPFR_RNDN);
        mpfr_init(r19713);
        mpfr_init_set_str(r19714, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19715, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19716);
        mpfr_init(r19717);
        mpfr_init(r19718);
        mpfr_init(r19719);
        mpfr_init(r19720);
        mpfr_init(r19721);
        mpfr_init(r19722);
        mpfr_init(r19723);
        mpfr_init(r19724);
        mpfr_init(r19725);
        mpfr_init(r19726);
        mpfr_init(r19727);
        mpfr_init_set_str(r19728, "-2.3665634646026902e-293", 10, MPFR_RNDN);
        mpfr_init(r19729);
        mpfr_init(r19730);
        mpfr_init(r19731);
        mpfr_init(r19732);
        mpfr_init(r19733);
        mpfr_init_set_str(r19734, "1.245134237263001e-160", 10, MPFR_RNDN);
        mpfr_init(r19735);
        mpfr_init(r19736);
        mpfr_init(r19737);
        mpfr_init(r19738);
        mpfr_init(r19739);
        mpfr_init_set_str(r19740, "2.56976209479494e+109", 10, MPFR_RNDN);
        mpfr_init(r19741);
        mpfr_init(r19742);
        mpfr_init(r19743);
        mpfr_init(r19744);
        mpfr_init(r19745);
        mpfr_init(r19746);
        mpfr_init(r19747);
        mpfr_init(r19748);
        mpfr_init(r19749);
        mpfr_init(r19750);
        mpfr_init(r19751);
        mpfr_init(r19752);
        mpfr_init(r19753);
        mpfr_init(r19754);
        mpfr_init(r19755);
        mpfr_init(r19756);
        mpfr_init(r19757);
        mpfr_init(r19758);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19711, re, MPFR_RNDN);
        ;
        mpfr_set_si(r19713, mpfr_cmp(r19711, r19712) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r19716, im, MPFR_RNDN);
        mpfr_mul(r19717, r19715, r19716, MPFR_RNDN);
        mpfr_mul(r19718, r19717, r19716, MPFR_RNDN);
        mpfr_sqrt(r19719, r19718, MPFR_RNDN);
        mpfr_sqr(r19720, r19711, MPFR_RNDN);
        mpfr_mul(r19721, r19716, r19716, MPFR_RNDN);
        mpfr_add(r19722, r19720, r19721, MPFR_RNDN);
        mpfr_sqrt(r19723, r19722, MPFR_RNDN);
        mpfr_sub(r19724, r19723, r19711, MPFR_RNDN);
        mpfr_sqrt(r19725, r19724, MPFR_RNDN);
        mpfr_div(r19726, r19719, r19725, MPFR_RNDN);
        mpfr_mul(r19727, r19714, r19726, MPFR_RNDN);
        ;
        mpfr_set_si(r19729, mpfr_cmp(r19711, r19728) <= 0, MPFR_RNDN);
        mpfr_add(r19730, r19716, r19711, MPFR_RNDN);
        mpfr_mul(r19731, r19715, r19730, MPFR_RNDN);
        mpfr_sqrt(r19732, r19731, MPFR_RNDN);
        mpfr_mul(r19733, r19714, r19732, MPFR_RNDN);
        ;
        mpfr_set_si(r19735, mpfr_cmp(r19711, r19734) <= 0, MPFR_RNDN);
        mpfr_div(r19736, r19721, r19724, MPFR_RNDN);
        mpfr_mul(r19737, r19715, r19736, MPFR_RNDN);
        mpfr_sqrt(r19738, r19737, MPFR_RNDN);
        mpfr_mul(r19739, r19714, r19738, MPFR_RNDN);
        ;
        mpfr_set_si(r19741, mpfr_cmp(r19711, r19740) <= 0, MPFR_RNDN);
        mpfr_mul(r19742, r19711, r19711, MPFR_RNDN);
        mpfr_add(r19743, r19742, r19721, MPFR_RNDN);
        mpfr_sqrt(r19744, r19743, MPFR_RNDN);
        mpfr_mul(r19745, r19744, r19744, MPFR_RNDN); mpfr_mul(r19745, r19745, r19744, MPFR_RNDN);
        mpfr_cbrt(r19746, r19745, MPFR_RNDN);
        mpfr_add(r19747, r19746, r19711, MPFR_RNDN);
        mpfr_mul(r19748, r19715, r19747, MPFR_RNDN);
        mpfr_sqrt(r19749, r19748, MPFR_RNDN);
        mpfr_mul(r19750, r19714, r19749, MPFR_RNDN);
        mpfr_add(r19751, r19711, r19711, MPFR_RNDN);
        mpfr_mul(r19752, r19715, r19751, MPFR_RNDN);
        mpfr_sqrt(r19753, r19752, MPFR_RNDN);
        mpfr_mul(r19754, r19714, r19753, MPFR_RNDN);
        if (mpfr_get_si(r19741, MPFR_RNDN)) { mpfr_set(r19755, r19750, MPFR_RNDN); } else { mpfr_set(r19755, r19754, MPFR_RNDN); };
        if (mpfr_get_si(r19735, MPFR_RNDN)) { mpfr_set(r19756, r19739, MPFR_RNDN); } else { mpfr_set(r19756, r19755, MPFR_RNDN); };
        if (mpfr_get_si(r19729, MPFR_RNDN)) { mpfr_set(r19757, r19733, MPFR_RNDN); } else { mpfr_set(r19757, r19756, MPFR_RNDN); };
        if (mpfr_get_si(r19713, MPFR_RNDN)) { mpfr_set(r19758, r19727, MPFR_RNDN); } else { mpfr_set(r19758, r19757, MPFR_RNDN); };
        return mpfr_get_d(r19758, MPFR_RNDN);
}

static mpfr_t r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803, r19804, r19805, r19806;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19759);
        mpfr_init_set_str(r19760, "-1.5608708420993077e-235", 10, MPFR_RNDN);
        mpfr_init(r19761);
        mpfr_init_set_str(r19762, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r19763, "2.0", 10, MPFR_RNDN);
        mpfr_init(r19764);
        mpfr_init(r19765);
        mpfr_init(r19766);
        mpfr_init(r19767);
        mpfr_init(r19768);
        mpfr_init(r19769);
        mpfr_init(r19770);
        mpfr_init(r19771);
        mpfr_init(r19772);
        mpfr_init(r19773);
        mpfr_init(r19774);
        mpfr_init(r19775);
        mpfr_init_set_str(r19776, "-2.3665634646026902e-293", 10, MPFR_RNDN);
        mpfr_init(r19777);
        mpfr_init(r19778);
        mpfr_init(r19779);
        mpfr_init(r19780);
        mpfr_init(r19781);
        mpfr_init_set_str(r19782, "1.245134237263001e-160", 10, MPFR_RNDN);
        mpfr_init(r19783);
        mpfr_init(r19784);
        mpfr_init(r19785);
        mpfr_init(r19786);
        mpfr_init(r19787);
        mpfr_init_set_str(r19788, "2.56976209479494e+109", 10, MPFR_RNDN);
        mpfr_init(r19789);
        mpfr_init(r19790);
        mpfr_init(r19791);
        mpfr_init(r19792);
        mpfr_init(r19793);
        mpfr_init(r19794);
        mpfr_init(r19795);
        mpfr_init(r19796);
        mpfr_init(r19797);
        mpfr_init(r19798);
        mpfr_init(r19799);
        mpfr_init(r19800);
        mpfr_init(r19801);
        mpfr_init(r19802);
        mpfr_init(r19803);
        mpfr_init(r19804);
        mpfr_init(r19805);
        mpfr_init(r19806);
}

double f_dm(double re, double im) {
        mpfr_set_d(r19759, re, MPFR_RNDN);
        ;
        mpfr_set_si(r19761, mpfr_cmp(r19759, r19760) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_set_d(r19764, im, MPFR_RNDN);
        mpfr_mul(r19765, r19763, r19764, MPFR_RNDN);
        mpfr_mul(r19766, r19765, r19764, MPFR_RNDN);
        mpfr_sqrt(r19767, r19766, MPFR_RNDN);
        mpfr_sqr(r19768, r19759, MPFR_RNDN);
        mpfr_mul(r19769, r19764, r19764, MPFR_RNDN);
        mpfr_add(r19770, r19768, r19769, MPFR_RNDN);
        mpfr_sqrt(r19771, r19770, MPFR_RNDN);
        mpfr_sub(r19772, r19771, r19759, MPFR_RNDN);
        mpfr_sqrt(r19773, r19772, MPFR_RNDN);
        mpfr_div(r19774, r19767, r19773, MPFR_RNDN);
        mpfr_mul(r19775, r19762, r19774, MPFR_RNDN);
        ;
        mpfr_set_si(r19777, mpfr_cmp(r19759, r19776) <= 0, MPFR_RNDN);
        mpfr_add(r19778, r19764, r19759, MPFR_RNDN);
        mpfr_mul(r19779, r19763, r19778, MPFR_RNDN);
        mpfr_sqrt(r19780, r19779, MPFR_RNDN);
        mpfr_mul(r19781, r19762, r19780, MPFR_RNDN);
        ;
        mpfr_set_si(r19783, mpfr_cmp(r19759, r19782) <= 0, MPFR_RNDN);
        mpfr_div(r19784, r19769, r19772, MPFR_RNDN);
        mpfr_mul(r19785, r19763, r19784, MPFR_RNDN);
        mpfr_sqrt(r19786, r19785, MPFR_RNDN);
        mpfr_mul(r19787, r19762, r19786, MPFR_RNDN);
        ;
        mpfr_set_si(r19789, mpfr_cmp(r19759, r19788) <= 0, MPFR_RNDN);
        mpfr_mul(r19790, r19759, r19759, MPFR_RNDN);
        mpfr_add(r19791, r19790, r19769, MPFR_RNDN);
        mpfr_sqrt(r19792, r19791, MPFR_RNDN);
        mpfr_mul(r19793, r19792, r19792, MPFR_RNDN); mpfr_mul(r19793, r19793, r19792, MPFR_RNDN);
        mpfr_cbrt(r19794, r19793, MPFR_RNDN);
        mpfr_add(r19795, r19794, r19759, MPFR_RNDN);
        mpfr_mul(r19796, r19763, r19795, MPFR_RNDN);
        mpfr_sqrt(r19797, r19796, MPFR_RNDN);
        mpfr_mul(r19798, r19762, r19797, MPFR_RNDN);
        mpfr_add(r19799, r19759, r19759, MPFR_RNDN);
        mpfr_mul(r19800, r19763, r19799, MPFR_RNDN);
        mpfr_sqrt(r19801, r19800, MPFR_RNDN);
        mpfr_mul(r19802, r19762, r19801, MPFR_RNDN);
        if (mpfr_get_si(r19789, MPFR_RNDN)) { mpfr_set(r19803, r19798, MPFR_RNDN); } else { mpfr_set(r19803, r19802, MPFR_RNDN); };
        if (mpfr_get_si(r19783, MPFR_RNDN)) { mpfr_set(r19804, r19787, MPFR_RNDN); } else { mpfr_set(r19804, r19803, MPFR_RNDN); };
        if (mpfr_get_si(r19777, MPFR_RNDN)) { mpfr_set(r19805, r19781, MPFR_RNDN); } else { mpfr_set(r19805, r19804, MPFR_RNDN); };
        if (mpfr_get_si(r19761, MPFR_RNDN)) { mpfr_set(r19806, r19775, MPFR_RNDN); } else { mpfr_set(r19806, r19805, MPFR_RNDN); };
        return mpfr_get_d(r19806, MPFR_RNDN);
}