#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 r7868 = 0.5f;
        float r7869 = 2.0f;
        float r7870 = re;
        float r7871 = r7870 * r7870;
        float r7872 = im;
        float r7873 = r7872 * r7872;
        float r7874 = r7871 + r7873;
        float r7875 = sqrt(r7874);
        float r7876 = r7875 - r7870;
        float r7877 = r7869 * r7876;
        float r7878 = sqrt(r7877);
        float r7879 = r7868 * r7878;
        return r7879;
}

double f_id(double re, double im) {
        double r7880 = 0.5;
        double r7881 = 2.0;
        double r7882 = re;
        double r7883 = r7882 * r7882;
        double r7884 = im;
        double r7885 = r7884 * r7884;
        double r7886 = r7883 + r7885;
        double r7887 = sqrt(r7886);
        double r7888 = r7887 - r7882;
        double r7889 = r7881 * r7888;
        double r7890 = sqrt(r7889);
        double r7891 = r7880 * r7890;
        return r7891;
}


double f_of(float re, float im) {
        float r7892 = re;
        float r7893 = im;
        float r7894 = hypot(r7892, r7893);
        float r7895 = r7894 - r7892;
        float r7896 = 2.0f;
        float r7897 = r7895 * r7896;
        float r7898 = sqrt(r7897);
        float r7899 = 0.5f;
        float r7900 = r7898 * r7899;
        return r7900;
}

double f_od(double re, double im) {
        double r7901 = re;
        double r7902 = im;
        double r7903 = hypot(r7901, r7902);
        double r7904 = r7903 - r7901;
        double r7905 = 2.0;
        double r7906 = r7904 * r7905;
        double r7907 = sqrt(r7906);
        double r7908 = 0.5;
        double r7909 = r7907 * r7908;
        return r7909;
}

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 r7910, r7911, r7912, r7913, r7914, r7915, r7916, r7917, r7918, r7919, r7920, r7921;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3664);
        mpfr_init_set_str(r7910, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r7911, "2.0", 10, MPFR_RNDN);
        mpfr_init(r7912);
        mpfr_init(r7913);
        mpfr_init(r7914);
        mpfr_init(r7915);
        mpfr_init(r7916);
        mpfr_init(r7917);
        mpfr_init(r7918);
        mpfr_init(r7919);
        mpfr_init(r7920);
        mpfr_init(r7921);
}

double f_im(double re, double im) {
        ;
        ;
        mpfr_set_d(r7912, re, MPFR_RNDN);
        mpfr_mul(r7913, r7912, r7912, MPFR_RNDN);
        mpfr_set_d(r7914, im, MPFR_RNDN);
        mpfr_mul(r7915, r7914, r7914, MPFR_RNDN);
        mpfr_add(r7916, r7913, r7915, MPFR_RNDN);
        mpfr_sqrt(r7917, r7916, MPFR_RNDN);
        mpfr_sub(r7918, r7917, r7912, MPFR_RNDN);
        mpfr_mul(r7919, r7911, r7918, MPFR_RNDN);
        mpfr_sqrt(r7920, r7919, MPFR_RNDN);
        mpfr_mul(r7921, r7910, r7920, MPFR_RNDN);
        return mpfr_get_d(r7921, MPFR_RNDN);
}

static mpfr_t r7922, r7923, r7924, r7925, r7926, r7927, r7928, r7929, r7930;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3664);
        mpfr_init(r7922);
        mpfr_init(r7923);
        mpfr_init(r7924);
        mpfr_init(r7925);
        mpfr_init_set_str(r7926, "2.0", 10, MPFR_RNDN);
        mpfr_init(r7927);
        mpfr_init(r7928);
        mpfr_init_set_str(r7929, "0.5", 10, MPFR_RNDN);
        mpfr_init(r7930);
}

double f_fm(double re, double im) {
        mpfr_set_d(r7922, re, MPFR_RNDN);
        mpfr_set_d(r7923, im, MPFR_RNDN);
        mpfr_hypot(r7924, r7922, r7923, MPFR_RNDN);
        mpfr_sub(r7925, r7924, r7922, MPFR_RNDN);
        ;
        mpfr_mul(r7927, r7925, r7926, MPFR_RNDN);
        mpfr_sqrt(r7928, r7927, MPFR_RNDN);
        ;
        mpfr_mul(r7930, r7928, r7929, MPFR_RNDN);
        return mpfr_get_d(r7930, MPFR_RNDN);
}

static mpfr_t r7931, r7932, r7933, r7934, r7935, r7936, r7937, r7938, r7939;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3664);
        mpfr_init(r7931);
        mpfr_init(r7932);
        mpfr_init(r7933);
        mpfr_init(r7934);
        mpfr_init_set_str(r7935, "2.0", 10, MPFR_RNDN);
        mpfr_init(r7936);
        mpfr_init(r7937);
        mpfr_init_set_str(r7938, "0.5", 10, MPFR_RNDN);
        mpfr_init(r7939);
}

double f_dm(double re, double im) {
        mpfr_set_d(r7931, re, MPFR_RNDN);
        mpfr_set_d(r7932, im, MPFR_RNDN);
        mpfr_hypot(r7933, r7931, r7932, MPFR_RNDN);
        mpfr_sub(r7934, r7933, r7931, MPFR_RNDN);
        ;
        mpfr_mul(r7936, r7934, r7935, MPFR_RNDN);
        mpfr_sqrt(r7937, r7936, MPFR_RNDN);
        ;
        mpfr_mul(r7939, r7937, r7938, MPFR_RNDN);
        return mpfr_get_d(r7939, MPFR_RNDN);
}