#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 r22956 = x_im;
        float r22957 = y_re;
        float r22958 = r22956 * r22957;
        float r22959 = x_re;
        float r22960 = y_im;
        float r22961 = r22959 * r22960;
        float r22962 = r22958 - r22961;
        float r22963 = r22957 * r22957;
        float r22964 = r22960 * r22960;
        float r22965 = r22963 + r22964;
        float r22966 = r22962 / r22965;
        return r22966;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r22967 = x_im;
        double r22968 = y_re;
        double r22969 = r22967 * r22968;
        double r22970 = x_re;
        double r22971 = y_im;
        double r22972 = r22970 * r22971;
        double r22973 = r22969 - r22972;
        double r22974 = r22968 * r22968;
        double r22975 = r22971 * r22971;
        double r22976 = r22974 + r22975;
        double r22977 = r22973 / r22976;
        return r22977;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r22978 = y_re;
        float r22979 = -2.1991094302601127e+162;
        bool r22980 = r22978 <= r22979;
        float r22981 = x_im;
        float r22982 = -r22981;
        float r22983 = y_im;
        float r22984 = hypot(r22978, r22983);
        float r22985 = r22982 / r22984;
        float r22986 = 3.9488752170732723e+114;
        bool r22987 = r22978 <= r22986;
        float r22988 = 1;
        float r22989 = r22988 / r22984;
        float r22990 = r22978 * r22981;
        float r22991 = x_re;
        float r22992 = r22991 * r22983;
        float r22993 = r22990 - r22992;
        float r22994 = r22984 / r22993;
        float r22995 = r22988 / r22994;
        float r22996 = r22989 * r22995;
        float r22997 = r22981 / r22984;
        float r22998 = r22987 ? r22996 : r22997;
        float r22999 = r22980 ? r22985 : r22998;
        return r22999;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r23000 = y_re;
        double r23001 = -2.1991094302601127e+162;
        bool r23002 = r23000 <= r23001;
        double r23003 = x_im;
        double r23004 = -r23003;
        double r23005 = y_im;
        double r23006 = hypot(r23000, r23005);
        double r23007 = r23004 / r23006;
        double r23008 = 3.9488752170732723e+114;
        bool r23009 = r23000 <= r23008;
        double r23010 = 1;
        double r23011 = r23010 / r23006;
        double r23012 = r23000 * r23003;
        double r23013 = x_re;
        double r23014 = r23013 * r23005;
        double r23015 = r23012 - r23014;
        double r23016 = r23006 / r23015;
        double r23017 = r23010 / r23016;
        double r23018 = r23011 * r23017;
        double r23019 = r23003 / r23006;
        double r23020 = r23009 ? r23018 : r23019;
        double r23021 = r23002 ? r23007 : r23020;
        return r23021;
}

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 r23022, r23023, r23024, r23025, r23026, r23027, r23028, r23029, r23030, r23031, r23032;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r23022);
        mpfr_init(r23023);
        mpfr_init(r23024);
        mpfr_init(r23025);
        mpfr_init(r23026);
        mpfr_init(r23027);
        mpfr_init(r23028);
        mpfr_init(r23029);
        mpfr_init(r23030);
        mpfr_init(r23031);
        mpfr_init(r23032);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r23022, x_im, MPFR_RNDN);
        mpfr_set_d(r23023, y_re, MPFR_RNDN);
        mpfr_mul(r23024, r23022, r23023, MPFR_RNDN);
        mpfr_set_d(r23025, x_re, MPFR_RNDN);
        mpfr_set_d(r23026, y_im, MPFR_RNDN);
        mpfr_mul(r23027, r23025, r23026, MPFR_RNDN);
        mpfr_sub(r23028, r23024, r23027, MPFR_RNDN);
        mpfr_mul(r23029, r23023, r23023, MPFR_RNDN);
        mpfr_mul(r23030, r23026, r23026, MPFR_RNDN);
        mpfr_add(r23031, r23029, r23030, MPFR_RNDN);
        mpfr_div(r23032, r23028, r23031, MPFR_RNDN);
        return mpfr_get_d(r23032, MPFR_RNDN);
}

static mpfr_t r23033, r23034, r23035, r23036, r23037, r23038, r23039, r23040, r23041, r23042, r23043, r23044, r23045, r23046, r23047, r23048, r23049, r23050, r23051, r23052, r23053, r23054;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r23033);
        mpfr_init_set_str(r23034, "-2.1991094302601127e+162", 10, MPFR_RNDN);
        mpfr_init(r23035);
        mpfr_init(r23036);
        mpfr_init(r23037);
        mpfr_init(r23038);
        mpfr_init(r23039);
        mpfr_init(r23040);
        mpfr_init_set_str(r23041, "3.9488752170732723e+114", 10, MPFR_RNDN);
        mpfr_init(r23042);
        mpfr_init_set_str(r23043, "1", 10, MPFR_RNDN);
        mpfr_init(r23044);
        mpfr_init(r23045);
        mpfr_init(r23046);
        mpfr_init(r23047);
        mpfr_init(r23048);
        mpfr_init(r23049);
        mpfr_init(r23050);
        mpfr_init(r23051);
        mpfr_init(r23052);
        mpfr_init(r23053);
        mpfr_init(r23054);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r23033, y_re, MPFR_RNDN);
        ;
        mpfr_set_si(r23035, mpfr_cmp(r23033, r23034) <= 0, MPFR_RNDN);
        mpfr_set_d(r23036, x_im, MPFR_RNDN);
        mpfr_neg(r23037, r23036, MPFR_RNDN);
        mpfr_set_d(r23038, y_im, MPFR_RNDN);
        mpfr_hypot(r23039, r23033, r23038, MPFR_RNDN);
        mpfr_div(r23040, r23037, r23039, MPFR_RNDN);
        ;
        mpfr_set_si(r23042, mpfr_cmp(r23033, r23041) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r23044, r23043, r23039, MPFR_RNDN);
        mpfr_mul(r23045, r23033, r23036, MPFR_RNDN);
        mpfr_set_d(r23046, x_re, MPFR_RNDN);
        mpfr_mul(r23047, r23046, r23038, MPFR_RNDN);
        mpfr_sub(r23048, r23045, r23047, MPFR_RNDN);
        mpfr_div(r23049, r23039, r23048, MPFR_RNDN);
        mpfr_div(r23050, r23043, r23049, MPFR_RNDN);
        mpfr_mul(r23051, r23044, r23050, MPFR_RNDN);
        mpfr_div(r23052, r23036, r23039, MPFR_RNDN);
        if (mpfr_get_si(r23042, MPFR_RNDN)) { mpfr_set(r23053, r23051, MPFR_RNDN); } else { mpfr_set(r23053, r23052, MPFR_RNDN); };
        if (mpfr_get_si(r23035, MPFR_RNDN)) { mpfr_set(r23054, r23040, MPFR_RNDN); } else { mpfr_set(r23054, r23053, MPFR_RNDN); };
        return mpfr_get_d(r23054, MPFR_RNDN);
}

static mpfr_t r23055, r23056, r23057, r23058, r23059, r23060, r23061, r23062, r23063, r23064, r23065, r23066, r23067, r23068, r23069, r23070, r23071, r23072, r23073, r23074, r23075, r23076;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r23055);
        mpfr_init_set_str(r23056, "-2.1991094302601127e+162", 10, MPFR_RNDN);
        mpfr_init(r23057);
        mpfr_init(r23058);
        mpfr_init(r23059);
        mpfr_init(r23060);
        mpfr_init(r23061);
        mpfr_init(r23062);
        mpfr_init_set_str(r23063, "3.9488752170732723e+114", 10, MPFR_RNDN);
        mpfr_init(r23064);
        mpfr_init_set_str(r23065, "1", 10, MPFR_RNDN);
        mpfr_init(r23066);
        mpfr_init(r23067);
        mpfr_init(r23068);
        mpfr_init(r23069);
        mpfr_init(r23070);
        mpfr_init(r23071);
        mpfr_init(r23072);
        mpfr_init(r23073);
        mpfr_init(r23074);
        mpfr_init(r23075);
        mpfr_init(r23076);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r23055, y_re, MPFR_RNDN);
        ;
        mpfr_set_si(r23057, mpfr_cmp(r23055, r23056) <= 0, MPFR_RNDN);
        mpfr_set_d(r23058, x_im, MPFR_RNDN);
        mpfr_neg(r23059, r23058, MPFR_RNDN);
        mpfr_set_d(r23060, y_im, MPFR_RNDN);
        mpfr_hypot(r23061, r23055, r23060, MPFR_RNDN);
        mpfr_div(r23062, r23059, r23061, MPFR_RNDN);
        ;
        mpfr_set_si(r23064, mpfr_cmp(r23055, r23063) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r23066, r23065, r23061, MPFR_RNDN);
        mpfr_mul(r23067, r23055, r23058, MPFR_RNDN);
        mpfr_set_d(r23068, x_re, MPFR_RNDN);
        mpfr_mul(r23069, r23068, r23060, MPFR_RNDN);
        mpfr_sub(r23070, r23067, r23069, MPFR_RNDN);
        mpfr_div(r23071, r23061, r23070, MPFR_RNDN);
        mpfr_div(r23072, r23065, r23071, MPFR_RNDN);
        mpfr_mul(r23073, r23066, r23072, MPFR_RNDN);
        mpfr_div(r23074, r23058, r23061, MPFR_RNDN);
        if (mpfr_get_si(r23064, MPFR_RNDN)) { mpfr_set(r23075, r23073, MPFR_RNDN); } else { mpfr_set(r23075, r23074, MPFR_RNDN); };
        if (mpfr_get_si(r23057, MPFR_RNDN)) { mpfr_set(r23076, r23062, MPFR_RNDN); } else { mpfr_set(r23076, r23075, MPFR_RNDN); };
        return mpfr_get_d(r23076, MPFR_RNDN);
}