#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "math.sqrt on complex, imaginary part, im greater than 0 branch";

double f_if(float re, float im) {
        float r20823 = 0.5;
        float r20824 = 2.0;
        float r20825 = re;
        float r20826 = r20825 * r20825;
        float r20827 = im;
        float r20828 = r20827 * r20827;
        float r20829 = r20826 - r20828;
        float r20830 = sqrt(r20829);
        float r20831 = r20830 + r20825;
        float r20832 = r20824 * r20831;
        float r20833 = sqrt(r20832);
        float r20834 = r20823 * r20833;
        return r20834;
}

double f_id(double re, double im) {
        double r20835 = 0.5;
        double r20836 = 2.0;
        double r20837 = re;
        double r20838 = r20837 * r20837;
        double r20839 = im;
        double r20840 = r20839 * r20839;
        double r20841 = r20838 - r20840;
        double r20842 = sqrt(r20841);
        double r20843 = r20842 + r20837;
        double r20844 = r20836 * r20843;
        double r20845 = sqrt(r20844);
        double r20846 = r20835 * r20845;
        return r20846;
}


double f_of(float re, float im) {
        float r20847 = 0.5;
        float r20848 = 2.0;
        float r20849 = re;
        float r20850 = im;
        float r20851 = r20849 + r20850;
        float r20852 = sqrt(r20851);
        float r20853 = r20849 - r20850;
        float r20854 = sqrt(r20853);
        float r20855 = r20852 * r20854;
        float r20856 = r20855 + r20849;
        float r20857 = r20848 * r20856;
        float r20858 = sqrt(r20857);
        float r20859 = r20847 * r20858;
        return r20859;
}

double f_od(double re, double im) {
        double r20860 = 0.5;
        double r20861 = 2.0;
        double r20862 = re;
        double r20863 = im;
        double r20864 = r20862 + r20863;
        double r20865 = sqrt(r20864);
        double r20866 = r20862 - r20863;
        double r20867 = sqrt(r20866);
        double r20868 = r20865 * r20867;
        double r20869 = r20868 + r20862;
        double r20870 = r20861 * r20869;
        double r20871 = sqrt(r20870);
        double r20872 = r20860 * r20871;
        return r20872;
}

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 r20873, r20874, r20875, r20876, r20877, r20878, r20879, r20880, r20881, r20882, r20883, r20884;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(4240);
        mpfr_init_set_str(r20873, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20874, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20875);
        mpfr_init(r20876);
        mpfr_init(r20877);
        mpfr_init(r20878);
        mpfr_init(r20879);
        mpfr_init(r20880);
        mpfr_init(r20881);
        mpfr_init(r20882);
        mpfr_init(r20883);
        mpfr_init(r20884);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r20875, re, MPFR_RNDN);
        mpfr_mul(r20876, r20875, r20875, MPFR_RNDN);
        mpfr_set_d(r20877, im, MPFR_RNDN);
        mpfr_mul(r20878, r20877, r20877, MPFR_RNDN);
        mpfr_sub(r20879, r20876, r20878, MPFR_RNDN);
        mpfr_sqrt(r20880, r20879, MPFR_RNDN);
        mpfr_add(r20881, r20880, r20875, MPFR_RNDN);
        mpfr_mul(r20882, r20874, r20881, MPFR_RNDN);
        mpfr_sqrt(r20883, r20882, MPFR_RNDN);
        mpfr_mul(r20884, r20873, r20883, MPFR_RNDN);
        return mpfr_get_d(r20884, MPFR_RNDN);
}

static mpfr_t r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(4240);
        mpfr_init_set_str(r20885, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20886, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20887);
        mpfr_init(r20888);
        mpfr_init(r20889);
        mpfr_init(r20890);
        mpfr_init(r20891);
        mpfr_init(r20892);
        mpfr_init(r20893);
        mpfr_init(r20894);
        mpfr_init(r20895);
        mpfr_init(r20896);
        mpfr_init(r20897);
}

double f_fm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20887, re, MPFR_RNDN);
        mpfr_set_d(r20888, im, MPFR_RNDN);
        mpfr_add(r20889, r20887, r20888, MPFR_RNDN);
        mpfr_sqrt(r20890, r20889, MPFR_RNDN);
        mpfr_sub(r20891, r20887, r20888, MPFR_RNDN);
        mpfr_sqrt(r20892, r20891, MPFR_RNDN);
        mpfr_mul(r20893, r20890, r20892, MPFR_RNDN);
        mpfr_add(r20894, r20893, r20887, MPFR_RNDN);
        mpfr_mul(r20895, r20886, r20894, MPFR_RNDN);
        mpfr_sqrt(r20896, r20895, MPFR_RNDN);
        mpfr_mul(r20897, r20885, r20896, MPFR_RNDN);
        return mpfr_get_d(r20897, MPFR_RNDN);
}

static mpfr_t r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(4240);
        mpfr_init_set_str(r20898, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r20899, "2.0", 10, MPFR_RNDN);
        mpfr_init(r20900);
        mpfr_init(r20901);
        mpfr_init(r20902);
        mpfr_init(r20903);
        mpfr_init(r20904);
        mpfr_init(r20905);
        mpfr_init(r20906);
        mpfr_init(r20907);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init(r20910);
}

double f_dm(double re, double im) {
        ;
        ;
        mpfr_set_d(r20900, re, MPFR_RNDN);
        mpfr_set_d(r20901, im, MPFR_RNDN);
        mpfr_add(r20902, r20900, r20901, MPFR_RNDN);
        mpfr_sqrt(r20903, r20902, MPFR_RNDN);
        mpfr_sub(r20904, r20900, r20901, MPFR_RNDN);
        mpfr_sqrt(r20905, r20904, MPFR_RNDN);
        mpfr_mul(r20906, r20903, r20905, MPFR_RNDN);
        mpfr_add(r20907, r20906, r20900, MPFR_RNDN);
        mpfr_mul(r20908, r20899, r20907, MPFR_RNDN);
        mpfr_sqrt(r20909, r20908, MPFR_RNDN);
        mpfr_mul(r20910, r20898, r20909, MPFR_RNDN);
        return mpfr_get_d(r20910, MPFR_RNDN);
}