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

char *name = "_divideComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r25998 = x_im;
        float r25999 = y_re;
        float r26000 = r25998 * r25999;
        float r26001 = x_re;
        float r26002 = y_im;
        float r26003 = r26001 * r26002;
        float r26004 = r26000 - r26003;
        float r26005 = r25999 * r25999;
        float r26006 = r26002 * r26002;
        float r26007 = r26005 + r26006;
        float r26008 = r26004 / r26007;
        return r26008;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r26009 = x_im;
        double r26010 = y_re;
        double r26011 = r26009 * r26010;
        double r26012 = x_re;
        double r26013 = y_im;
        double r26014 = r26012 * r26013;
        double r26015 = r26011 - r26014;
        double r26016 = r26010 * r26010;
        double r26017 = r26013 * r26013;
        double r26018 = r26016 + r26017;
        double r26019 = r26015 / r26018;
        return r26019;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r26020 = y_re;
        float r26021 = -2.2149249625222455e+181;
        bool r26022 = r26020 <= r26021;
        float r26023 = x_im;
        float r26024 = -r26023;
        float r26025 = y_im;
        float r26026 = hypot(r26020, r26025);
        float r26027 = r26024 / r26026;
        float r26028 = 2.143387036101042e+85;
        bool r26029 = r26020 <= r26028;
        float r26030 = 1;
        float r26031 = r26030 / r26026;
        float r26032 = r26020 * r26023;
        float r26033 = x_re;
        float r26034 = r26033 * r26025;
        float r26035 = r26032 - r26034;
        float r26036 = r26035 * r26031;
        float r26037 = r26031 * r26036;
        float r26038 = r26023 / r26026;
        float r26039 = r26029 ? r26037 : r26038;
        float r26040 = r26022 ? r26027 : r26039;
        return r26040;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r26041 = y_re;
        double r26042 = -2.2149249625222455e+181;
        bool r26043 = r26041 <= r26042;
        double r26044 = x_im;
        double r26045 = -r26044;
        double r26046 = y_im;
        double r26047 = hypot(r26041, r26046);
        double r26048 = r26045 / r26047;
        double r26049 = 2.143387036101042e+85;
        bool r26050 = r26041 <= r26049;
        double r26051 = 1;
        double r26052 = r26051 / r26047;
        double r26053 = r26041 * r26044;
        double r26054 = x_re;
        double r26055 = r26054 * r26046;
        double r26056 = r26053 - r26055;
        double r26057 = r26056 * r26052;
        double r26058 = r26052 * r26057;
        double r26059 = r26044 / r26047;
        double r26060 = r26050 ? r26058 : r26059;
        double r26061 = r26043 ? r26048 : r26060;
        return r26061;
}

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 r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26062);
        mpfr_init(r26063);
        mpfr_init(r26064);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init(r26068);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r26062, x_im, MPFR_RNDN);
        mpfr_set_d(r26063, y_re, MPFR_RNDN);
        mpfr_mul(r26064, r26062, r26063, MPFR_RNDN);
        mpfr_set_d(r26065, x_re, MPFR_RNDN);
        mpfr_set_d(r26066, y_im, MPFR_RNDN);
        mpfr_mul(r26067, r26065, r26066, MPFR_RNDN);
        mpfr_sub(r26068, r26064, r26067, MPFR_RNDN);
        mpfr_mul(r26069, r26063, r26063, MPFR_RNDN);
        mpfr_mul(r26070, r26066, r26066, MPFR_RNDN);
        mpfr_add(r26071, r26069, r26070, MPFR_RNDN);
        mpfr_div(r26072, r26068, r26071, MPFR_RNDN);
        return mpfr_get_d(r26072, MPFR_RNDN);
}

static mpfr_t r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26073);
        mpfr_init_set_str(r26074, "-2.2149249625222455e+181", 10, MPFR_RNDN);
        mpfr_init(r26075);
        mpfr_init(r26076);
        mpfr_init(r26077);
        mpfr_init(r26078);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init_set_str(r26081, "2.143387036101042e+85", 10, MPFR_RNDN);
        mpfr_init(r26082);
        mpfr_init_set_str(r26083, "1", 10, MPFR_RNDN);
        mpfr_init(r26084);
        mpfr_init(r26085);
        mpfr_init(r26086);
        mpfr_init(r26087);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
        mpfr_init(r26091);
        mpfr_init(r26092);
        mpfr_init(r26093);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r26073, y_re, MPFR_RNDN);
        ;
        mpfr_set_si(r26075, mpfr_cmp(r26073, r26074) <= 0, MPFR_RNDN);
        mpfr_set_d(r26076, x_im, MPFR_RNDN);
        mpfr_neg(r26077, r26076, MPFR_RNDN);
        mpfr_set_d(r26078, y_im, MPFR_RNDN);
        mpfr_hypot(r26079, r26073, r26078, MPFR_RNDN);
        mpfr_div(r26080, r26077, r26079, MPFR_RNDN);
        ;
        mpfr_set_si(r26082, mpfr_cmp(r26073, r26081) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26084, r26083, r26079, MPFR_RNDN);
        mpfr_mul(r26085, r26073, r26076, MPFR_RNDN);
        mpfr_set_d(r26086, x_re, MPFR_RNDN);
        mpfr_mul(r26087, r26086, r26078, MPFR_RNDN);
        mpfr_sub(r26088, r26085, r26087, MPFR_RNDN);
        mpfr_mul(r26089, r26088, r26084, MPFR_RNDN);
        mpfr_mul(r26090, r26084, r26089, MPFR_RNDN);
        mpfr_div(r26091, r26076, r26079, MPFR_RNDN);
        if (mpfr_get_si(r26082, MPFR_RNDN)) { mpfr_set(r26092, r26090, MPFR_RNDN); } else { mpfr_set(r26092, r26091, MPFR_RNDN); };
        if (mpfr_get_si(r26075, MPFR_RNDN)) { mpfr_set(r26093, r26080, MPFR_RNDN); } else { mpfr_set(r26093, r26092, MPFR_RNDN); };
        return mpfr_get_d(r26093, MPFR_RNDN);
}

static mpfr_t r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26094);
        mpfr_init_set_str(r26095, "-2.2149249625222455e+181", 10, MPFR_RNDN);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init(r26098);
        mpfr_init(r26099);
        mpfr_init(r26100);
        mpfr_init(r26101);
        mpfr_init_set_str(r26102, "2.143387036101042e+85", 10, MPFR_RNDN);
        mpfr_init(r26103);
        mpfr_init_set_str(r26104, "1", 10, MPFR_RNDN);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init(r26107);
        mpfr_init(r26108);
        mpfr_init(r26109);
        mpfr_init(r26110);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init(r26113);
        mpfr_init(r26114);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r26094, y_re, MPFR_RNDN);
        ;
        mpfr_set_si(r26096, mpfr_cmp(r26094, r26095) <= 0, MPFR_RNDN);
        mpfr_set_d(r26097, x_im, MPFR_RNDN);
        mpfr_neg(r26098, r26097, MPFR_RNDN);
        mpfr_set_d(r26099, y_im, MPFR_RNDN);
        mpfr_hypot(r26100, r26094, r26099, MPFR_RNDN);
        mpfr_div(r26101, r26098, r26100, MPFR_RNDN);
        ;
        mpfr_set_si(r26103, mpfr_cmp(r26094, r26102) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26105, r26104, r26100, MPFR_RNDN);
        mpfr_mul(r26106, r26094, r26097, MPFR_RNDN);
        mpfr_set_d(r26107, x_re, MPFR_RNDN);
        mpfr_mul(r26108, r26107, r26099, MPFR_RNDN);
        mpfr_sub(r26109, r26106, r26108, MPFR_RNDN);
        mpfr_mul(r26110, r26109, r26105, MPFR_RNDN);
        mpfr_mul(r26111, r26105, r26110, MPFR_RNDN);
        mpfr_div(r26112, r26097, r26100, MPFR_RNDN);
        if (mpfr_get_si(r26103, MPFR_RNDN)) { mpfr_set(r26113, r26111, MPFR_RNDN); } else { mpfr_set(r26113, r26112, MPFR_RNDN); };
        if (mpfr_get_si(r26096, MPFR_RNDN)) { mpfr_set(r26114, r26101, MPFR_RNDN); } else { mpfr_set(r26114, r26113, MPFR_RNDN); };
        return mpfr_get_d(r26114, MPFR_RNDN);
}