#include #include #include #include #include char *name = "Complex division, imag part"; double f_if(float a, float b, float c, float d) { float r23635 = b; float r23636 = c; float r23637 = r23635 * r23636; float r23638 = a; float r23639 = d; float r23640 = r23638 * r23639; float r23641 = r23637 - r23640; float r23642 = r23636 * r23636; float r23643 = r23639 * r23639; float r23644 = r23642 + r23643; float r23645 = r23641 / r23644; return r23645; } double f_id(double a, double b, double c, double d) { double r23646 = b; double r23647 = c; double r23648 = r23646 * r23647; double r23649 = a; double r23650 = d; double r23651 = r23649 * r23650; double r23652 = r23648 - r23651; double r23653 = r23647 * r23647; double r23654 = r23650 * r23650; double r23655 = r23653 + r23654; double r23656 = r23652 / r23655; return r23656; } double f_of(float a, float b, float c, float d) { float r23657 = c; float r23658 = -3.520039360321946e+135; bool r23659 = r23657 <= r23658; float r23660 = b; float r23661 = -r23660; float r23662 = d; float r23663 = hypot(r23657, r23662); float r23664 = r23661 / r23663; float r23665 = 1.0151773538233026e+121; bool r23666 = r23657 <= r23665; float r23667 = 1; float r23668 = r23667 / r23663; float r23669 = a; float r23670 = r23669 * r23662; float r23671 = -r23670; float r23672 = fma(r23657, r23660, r23671); float r23673 = r23672 / r23663; float r23674 = r23668 * r23673; float r23675 = r23669 / r23657; float r23676 = -r23662; float r23677 = fma(r23675, r23676, r23660); float r23678 = r23677 / r23663; float r23679 = r23666 ? r23674 : r23678; float r23680 = r23659 ? r23664 : r23679; return r23680; } double f_od(double a, double b, double c, double d) { double r23681 = c; double r23682 = -3.520039360321946e+135; bool r23683 = r23681 <= r23682; double r23684 = b; double r23685 = -r23684; double r23686 = d; double r23687 = hypot(r23681, r23686); double r23688 = r23685 / r23687; double r23689 = 1.0151773538233026e+121; bool r23690 = r23681 <= r23689; double r23691 = 1; double r23692 = r23691 / r23687; double r23693 = a; double r23694 = r23693 * r23686; double r23695 = -r23694; double r23696 = fma(r23681, r23684, r23695); double r23697 = r23696 / r23687; double r23698 = r23692 * r23697; double r23699 = r23693 / r23681; double r23700 = -r23686; double r23701 = fma(r23699, r23700, r23684); double r23702 = r23701 / r23687; double r23703 = r23690 ? r23698 : r23702; double r23704 = r23683 ? r23688 : r23703; return r23704; } 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 r23705, r23706, r23707, r23708, r23709, r23710, r23711, r23712, r23713, r23714, r23715; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23705); mpfr_init(r23706); mpfr_init(r23707); mpfr_init(r23708); mpfr_init(r23709); mpfr_init(r23710); mpfr_init(r23711); mpfr_init(r23712); mpfr_init(r23713); mpfr_init(r23714); mpfr_init(r23715); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r23705, b, MPFR_RNDN); mpfr_set_d(r23706, c, MPFR_RNDN); mpfr_mul(r23707, r23705, r23706, MPFR_RNDN); mpfr_set_d(r23708, a, MPFR_RNDN); mpfr_set_d(r23709, d, MPFR_RNDN); mpfr_mul(r23710, r23708, r23709, MPFR_RNDN); mpfr_sub(r23711, r23707, r23710, MPFR_RNDN); mpfr_mul(r23712, r23706, r23706, MPFR_RNDN); mpfr_mul(r23713, r23709, r23709, MPFR_RNDN); mpfr_add(r23714, r23712, r23713, MPFR_RNDN); mpfr_div(r23715, r23711, r23714, MPFR_RNDN); return mpfr_get_d(r23715, MPFR_RNDN); } static mpfr_t r23716, r23717, r23718, r23719, r23720, r23721, r23722, r23723, r23724, r23725, r23726, r23727, r23728, r23729, r23730, r23731, r23732, r23733, r23734, r23735, r23736, r23737, r23738, r23739; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23716); mpfr_init_set_str(r23717, "-3.520039360321946e+135", 10, MPFR_RNDN); mpfr_init(r23718); mpfr_init(r23719); mpfr_init(r23720); mpfr_init(r23721); mpfr_init(r23722); mpfr_init(r23723); mpfr_init_set_str(r23724, "1.0151773538233026e+121", 10, MPFR_RNDN); mpfr_init(r23725); mpfr_init_set_str(r23726, "1", 10, MPFR_RNDN); mpfr_init(r23727); mpfr_init(r23728); mpfr_init(r23729); mpfr_init(r23730); mpfr_init(r23731); mpfr_init(r23732); mpfr_init(r23733); mpfr_init(r23734); mpfr_init(r23735); mpfr_init(r23736); mpfr_init(r23737); mpfr_init(r23738); mpfr_init(r23739); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r23716, c, MPFR_RNDN); ; mpfr_set_si(r23718, mpfr_cmp(r23716, r23717) <= 0, MPFR_RNDN); mpfr_set_d(r23719, b, MPFR_RNDN); mpfr_neg(r23720, r23719, MPFR_RNDN); mpfr_set_d(r23721, d, MPFR_RNDN); mpfr_hypot(r23722, r23716, r23721, MPFR_RNDN); mpfr_div(r23723, r23720, r23722, MPFR_RNDN); ; mpfr_set_si(r23725, mpfr_cmp(r23716, r23724) <= 0, MPFR_RNDN); ; mpfr_div(r23727, r23726, r23722, MPFR_RNDN); mpfr_set_d(r23728, a, MPFR_RNDN); mpfr_mul(r23729, r23728, r23721, MPFR_RNDN); mpfr_neg(r23730, r23729, MPFR_RNDN); mpfr_fma(r23731, r23716, r23719, r23730, MPFR_RNDN); mpfr_div(r23732, r23731, r23722, MPFR_RNDN); mpfr_mul(r23733, r23727, r23732, MPFR_RNDN); mpfr_div(r23734, r23728, r23716, MPFR_RNDN); mpfr_neg(r23735, r23721, MPFR_RNDN); mpfr_fma(r23736, r23734, r23735, r23719, MPFR_RNDN); mpfr_div(r23737, r23736, r23722, MPFR_RNDN); if (mpfr_get_si(r23725, MPFR_RNDN)) { mpfr_set(r23738, r23733, MPFR_RNDN); } else { mpfr_set(r23738, r23737, MPFR_RNDN); }; if (mpfr_get_si(r23718, MPFR_RNDN)) { mpfr_set(r23739, r23723, MPFR_RNDN); } else { mpfr_set(r23739, r23738, MPFR_RNDN); }; return mpfr_get_d(r23739, MPFR_RNDN); } static mpfr_t r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747, r23748, r23749, r23750, r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23740); mpfr_init_set_str(r23741, "-3.520039360321946e+135", 10, MPFR_RNDN); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init(r23745); mpfr_init(r23746); mpfr_init(r23747); mpfr_init_set_str(r23748, "1.0151773538233026e+121", 10, MPFR_RNDN); mpfr_init(r23749); mpfr_init_set_str(r23750, "1", 10, MPFR_RNDN); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init(r23754); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r23740, c, MPFR_RNDN); ; mpfr_set_si(r23742, mpfr_cmp(r23740, r23741) <= 0, MPFR_RNDN); mpfr_set_d(r23743, b, MPFR_RNDN); mpfr_neg(r23744, r23743, MPFR_RNDN); mpfr_set_d(r23745, d, MPFR_RNDN); mpfr_hypot(r23746, r23740, r23745, MPFR_RNDN); mpfr_div(r23747, r23744, r23746, MPFR_RNDN); ; mpfr_set_si(r23749, mpfr_cmp(r23740, r23748) <= 0, MPFR_RNDN); ; mpfr_div(r23751, r23750, r23746, MPFR_RNDN); mpfr_set_d(r23752, a, MPFR_RNDN); mpfr_mul(r23753, r23752, r23745, MPFR_RNDN); mpfr_neg(r23754, r23753, MPFR_RNDN); mpfr_fma(r23755, r23740, r23743, r23754, MPFR_RNDN); mpfr_div(r23756, r23755, r23746, MPFR_RNDN); mpfr_mul(r23757, r23751, r23756, MPFR_RNDN); mpfr_div(r23758, r23752, r23740, MPFR_RNDN); mpfr_neg(r23759, r23745, MPFR_RNDN); mpfr_fma(r23760, r23758, r23759, r23743, MPFR_RNDN); mpfr_div(r23761, r23760, r23746, MPFR_RNDN); if (mpfr_get_si(r23749, MPFR_RNDN)) { mpfr_set(r23762, r23757, MPFR_RNDN); } else { mpfr_set(r23762, r23761, MPFR_RNDN); }; if (mpfr_get_si(r23742, MPFR_RNDN)) { mpfr_set(r23763, r23747, MPFR_RNDN); } else { mpfr_set(r23763, r23762, MPFR_RNDN); }; return mpfr_get_d(r23763, MPFR_RNDN); }