#include #include #include #include #include char *name = "Complex division, real part"; double f_if(float a, float b, float c, float d) { float r23524 = a; float r23525 = c; float r23526 = r23524 * r23525; float r23527 = b; float r23528 = d; float r23529 = r23527 * r23528; float r23530 = r23526 + r23529; float r23531 = r23525 * r23525; float r23532 = r23528 * r23528; float r23533 = r23531 + r23532; float r23534 = r23530 / r23533; return r23534; } double f_id(double a, double b, double c, double d) { double r23535 = a; double r23536 = c; double r23537 = r23535 * r23536; double r23538 = b; double r23539 = d; double r23540 = r23538 * r23539; double r23541 = r23537 + r23540; double r23542 = r23536 * r23536; double r23543 = r23539 * r23539; double r23544 = r23542 + r23543; double r23545 = r23541 / r23544; return r23545; } double f_of(float a, float b, float c, float d) { float r23546 = d; float r23547 = -3.2769180631198613e+229; bool r23548 = r23546 <= r23547; float r23549 = b; float r23550 = -r23549; float r23551 = c; float r23552 = hypot(r23546, r23551); float r23553 = r23550 / r23552; float r23554 = 7.015217207138249e+139; bool r23555 = r23546 <= r23554; float r23556 = 1; float r23557 = r23556 / r23552; float r23558 = sqrt(r23552); float r23559 = r23556 / r23558; float r23560 = a; float r23561 = r23551 * r23560; float r23562 = fma(r23549, r23546, r23561); float r23563 = r23562 / r23558; float r23564 = r23559 * r23563; float r23565 = r23557 * r23564; float r23566 = r23549 / r23552; float r23567 = r23555 ? r23565 : r23566; float r23568 = r23548 ? r23553 : r23567; return r23568; } double f_od(double a, double b, double c, double d) { double r23569 = d; double r23570 = -3.2769180631198613e+229; bool r23571 = r23569 <= r23570; double r23572 = b; double r23573 = -r23572; double r23574 = c; double r23575 = hypot(r23569, r23574); double r23576 = r23573 / r23575; double r23577 = 7.015217207138249e+139; bool r23578 = r23569 <= r23577; double r23579 = 1; double r23580 = r23579 / r23575; double r23581 = sqrt(r23575); double r23582 = r23579 / r23581; double r23583 = a; double r23584 = r23574 * r23583; double r23585 = fma(r23572, r23569, r23584); double r23586 = r23585 / r23581; double r23587 = r23582 * r23586; double r23588 = r23580 * r23587; double r23589 = r23572 / r23575; double r23590 = r23578 ? r23588 : r23589; double r23591 = r23571 ? r23576 : r23590; return r23591; } 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 r23592, r23593, r23594, r23595, r23596, r23597, r23598, r23599, r23600, r23601, r23602; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23592); mpfr_init(r23593); mpfr_init(r23594); mpfr_init(r23595); mpfr_init(r23596); mpfr_init(r23597); mpfr_init(r23598); mpfr_init(r23599); mpfr_init(r23600); mpfr_init(r23601); mpfr_init(r23602); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r23592, a, MPFR_RNDN); mpfr_set_d(r23593, c, MPFR_RNDN); mpfr_mul(r23594, r23592, r23593, MPFR_RNDN); mpfr_set_d(r23595, b, MPFR_RNDN); mpfr_set_d(r23596, d, MPFR_RNDN); mpfr_mul(r23597, r23595, r23596, MPFR_RNDN); mpfr_add(r23598, r23594, r23597, MPFR_RNDN); mpfr_mul(r23599, r23593, r23593, MPFR_RNDN); mpfr_mul(r23600, r23596, r23596, MPFR_RNDN); mpfr_add(r23601, r23599, r23600, MPFR_RNDN); mpfr_div(r23602, r23598, r23601, MPFR_RNDN); return mpfr_get_d(r23602, MPFR_RNDN); } static mpfr_t r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23603); mpfr_init_set_str(r23604, "-3.2769180631198613e+229", 10, MPFR_RNDN); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init_set_str(r23611, "7.015217207138249e+139", 10, MPFR_RNDN); mpfr_init(r23612); mpfr_init_set_str(r23613, "1", 10, MPFR_RNDN); mpfr_init(r23614); mpfr_init(r23615); mpfr_init(r23616); mpfr_init(r23617); mpfr_init(r23618); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r23603, d, MPFR_RNDN); ; mpfr_set_si(r23605, mpfr_cmp(r23603, r23604) <= 0, MPFR_RNDN); mpfr_set_d(r23606, b, MPFR_RNDN); mpfr_neg(r23607, r23606, MPFR_RNDN); mpfr_set_d(r23608, c, MPFR_RNDN); mpfr_hypot(r23609, r23603, r23608, MPFR_RNDN); mpfr_div(r23610, r23607, r23609, MPFR_RNDN); ; mpfr_set_si(r23612, mpfr_cmp(r23603, r23611) <= 0, MPFR_RNDN); ; mpfr_div(r23614, r23613, r23609, MPFR_RNDN); mpfr_sqrt(r23615, r23609, MPFR_RNDN); mpfr_div(r23616, r23613, r23615, MPFR_RNDN); mpfr_set_d(r23617, a, MPFR_RNDN); mpfr_mul(r23618, r23608, r23617, MPFR_RNDN); mpfr_fma(r23619, r23606, r23603, r23618, MPFR_RNDN); mpfr_div(r23620, r23619, r23615, MPFR_RNDN); mpfr_mul(r23621, r23616, r23620, MPFR_RNDN); mpfr_mul(r23622, r23614, r23621, MPFR_RNDN); mpfr_div(r23623, r23606, r23609, MPFR_RNDN); if (mpfr_get_si(r23612, MPFR_RNDN)) { mpfr_set(r23624, r23622, MPFR_RNDN); } else { mpfr_set(r23624, r23623, MPFR_RNDN); }; if (mpfr_get_si(r23605, MPFR_RNDN)) { mpfr_set(r23625, r23610, MPFR_RNDN); } else { mpfr_set(r23625, r23624, MPFR_RNDN); }; return mpfr_get_d(r23625, MPFR_RNDN); } static mpfr_t r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23626); mpfr_init_set_str(r23627, "-3.2769180631198613e+229", 10, MPFR_RNDN); mpfr_init(r23628); mpfr_init(r23629); mpfr_init(r23630); mpfr_init(r23631); mpfr_init(r23632); mpfr_init(r23633); mpfr_init_set_str(r23634, "7.015217207138249e+139", 10, MPFR_RNDN); mpfr_init(r23635); mpfr_init_set_str(r23636, "1", 10, MPFR_RNDN); mpfr_init(r23637); mpfr_init(r23638); mpfr_init(r23639); mpfr_init(r23640); mpfr_init(r23641); mpfr_init(r23642); mpfr_init(r23643); mpfr_init(r23644); mpfr_init(r23645); mpfr_init(r23646); mpfr_init(r23647); mpfr_init(r23648); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r23626, d, MPFR_RNDN); ; mpfr_set_si(r23628, mpfr_cmp(r23626, r23627) <= 0, MPFR_RNDN); mpfr_set_d(r23629, b, MPFR_RNDN); mpfr_neg(r23630, r23629, MPFR_RNDN); mpfr_set_d(r23631, c, MPFR_RNDN); mpfr_hypot(r23632, r23626, r23631, MPFR_RNDN); mpfr_div(r23633, r23630, r23632, MPFR_RNDN); ; mpfr_set_si(r23635, mpfr_cmp(r23626, r23634) <= 0, MPFR_RNDN); ; mpfr_div(r23637, r23636, r23632, MPFR_RNDN); mpfr_sqrt(r23638, r23632, MPFR_RNDN); mpfr_div(r23639, r23636, r23638, MPFR_RNDN); mpfr_set_d(r23640, a, MPFR_RNDN); mpfr_mul(r23641, r23631, r23640, MPFR_RNDN); mpfr_fma(r23642, r23629, r23626, r23641, MPFR_RNDN); mpfr_div(r23643, r23642, r23638, MPFR_RNDN); mpfr_mul(r23644, r23639, r23643, MPFR_RNDN); mpfr_mul(r23645, r23637, r23644, MPFR_RNDN); mpfr_div(r23646, r23629, r23632, MPFR_RNDN); if (mpfr_get_si(r23635, MPFR_RNDN)) { mpfr_set(r23647, r23645, MPFR_RNDN); } else { mpfr_set(r23647, r23646, MPFR_RNDN); }; if (mpfr_get_si(r23628, MPFR_RNDN)) { mpfr_set(r23648, r23633, MPFR_RNDN); } else { mpfr_set(r23648, r23647, MPFR_RNDN); }; return mpfr_get_d(r23648, MPFR_RNDN); }