#include #include #include #include #include char *name = "_divideComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r25487 = x_im; float r25488 = y_re; float r25489 = r25487 * r25488; float r25490 = x_re; float r25491 = y_im; float r25492 = r25490 * r25491; float r25493 = r25489 - r25492; float r25494 = r25488 * r25488; float r25495 = r25491 * r25491; float r25496 = r25494 + r25495; float r25497 = r25493 / r25496; return r25497; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r25498 = x_im; double r25499 = y_re; double r25500 = r25498 * r25499; double r25501 = x_re; double r25502 = y_im; double r25503 = r25501 * r25502; double r25504 = r25500 - r25503; double r25505 = r25499 * r25499; double r25506 = r25502 * r25502; double r25507 = r25505 + r25506; double r25508 = r25504 / r25507; return r25508; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r25509 = y_re; float r25510 = y_im; float r25511 = hypot(r25509, r25510); float r25512 = sqrt(r25511); float r25513 = r25509 / r25512; float r25514 = x_im; float r25515 = r25514 / r25512; float r25516 = -r25510; float r25517 = r25516 / r25511; float r25518 = x_re; float r25519 = r25517 * r25518; float r25520 = fma(r25513, r25515, r25519); float r25521 = r25520 / r25511; float r25522 = 0; float r25523 = r25521 + r25522; return r25523; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r25524 = y_re; double r25525 = y_im; double r25526 = hypot(r25524, r25525); double r25527 = sqrt(r25526); double r25528 = r25524 / r25527; double r25529 = x_im; double r25530 = r25529 / r25527; double r25531 = -r25525; double r25532 = r25531 / r25526; double r25533 = x_re; double r25534 = r25532 * r25533; double r25535 = fma(r25528, r25530, r25534); double r25536 = r25535 / r25526; double r25537 = 0; double r25538 = r25536 + r25537; return r25538; } 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 r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25539, x_im, MPFR_RNDN); mpfr_set_d(r25540, y_re, MPFR_RNDN); mpfr_mul(r25541, r25539, r25540, MPFR_RNDN); mpfr_set_d(r25542, x_re, MPFR_RNDN); mpfr_set_d(r25543, y_im, MPFR_RNDN); mpfr_mul(r25544, r25542, r25543, MPFR_RNDN); mpfr_sub(r25545, r25541, r25544, MPFR_RNDN); mpfr_mul(r25546, r25540, r25540, MPFR_RNDN); mpfr_mul(r25547, r25543, r25543, MPFR_RNDN); mpfr_add(r25548, r25546, r25547, MPFR_RNDN); mpfr_div(r25549, r25545, r25548, MPFR_RNDN); return mpfr_get_d(r25549, MPFR_RNDN); } static mpfr_t r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init(r25562); mpfr_init_set_str(r25563, "0", 10, MPFR_RNDN); mpfr_init(r25564); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25550, y_re, MPFR_RNDN); mpfr_set_d(r25551, y_im, MPFR_RNDN); mpfr_hypot(r25552, r25550, r25551, MPFR_RNDN); mpfr_sqrt(r25553, r25552, MPFR_RNDN); mpfr_div(r25554, r25550, r25553, MPFR_RNDN); mpfr_set_d(r25555, x_im, MPFR_RNDN); mpfr_div(r25556, r25555, r25553, MPFR_RNDN); mpfr_neg(r25557, r25551, MPFR_RNDN); mpfr_div(r25558, r25557, r25552, MPFR_RNDN); mpfr_set_d(r25559, x_re, MPFR_RNDN); mpfr_mul(r25560, r25558, r25559, MPFR_RNDN); mpfr_fma(r25561, r25554, r25556, r25560, MPFR_RNDN); mpfr_div(r25562, r25561, r25552, MPFR_RNDN); ; mpfr_add(r25564, r25562, r25563, MPFR_RNDN); return mpfr_get_d(r25564, MPFR_RNDN); } static mpfr_t r25565, r25566, r25567, r25568, r25569, r25570, r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); mpfr_init(r25574); mpfr_init(r25575); mpfr_init(r25576); mpfr_init(r25577); mpfr_init_set_str(r25578, "0", 10, MPFR_RNDN); mpfr_init(r25579); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25565, y_re, MPFR_RNDN); mpfr_set_d(r25566, y_im, MPFR_RNDN); mpfr_hypot(r25567, r25565, r25566, MPFR_RNDN); mpfr_sqrt(r25568, r25567, MPFR_RNDN); mpfr_div(r25569, r25565, r25568, MPFR_RNDN); mpfr_set_d(r25570, x_im, MPFR_RNDN); mpfr_div(r25571, r25570, r25568, MPFR_RNDN); mpfr_neg(r25572, r25566, MPFR_RNDN); mpfr_div(r25573, r25572, r25567, MPFR_RNDN); mpfr_set_d(r25574, x_re, MPFR_RNDN); mpfr_mul(r25575, r25573, r25574, MPFR_RNDN); mpfr_fma(r25576, r25569, r25571, r25575, MPFR_RNDN); mpfr_div(r25577, r25576, r25567, MPFR_RNDN); ; mpfr_add(r25579, r25577, r25578, MPFR_RNDN); return mpfr_get_d(r25579, MPFR_RNDN); }