#include #include #include #include #include char *name = "Complex division, real part"; double f_if(float a, float b, float c, float d) { float r23676 = a; float r23677 = c; float r23678 = r23676 * r23677; float r23679 = b; float r23680 = d; float r23681 = r23679 * r23680; float r23682 = r23678 + r23681; float r23683 = r23677 * r23677; float r23684 = r23680 * r23680; float r23685 = r23683 + r23684; float r23686 = r23682 / r23685; return r23686; } double f_id(double a, double b, double c, double d) { double r23687 = a; double r23688 = c; double r23689 = r23687 * r23688; double r23690 = b; double r23691 = d; double r23692 = r23690 * r23691; double r23693 = r23689 + r23692; double r23694 = r23688 * r23688; double r23695 = r23691 * r23691; double r23696 = r23694 + r23695; double r23697 = r23693 / r23696; return r23697; } double f_of(float a, float b, float c, float d) { float r23698 = c; float r23699 = -2.114794877571543e+145; bool r23700 = r23698 <= r23699; float r23701 = b; float r23702 = r23701 / r23698; float r23703 = d; float r23704 = -r23703; float r23705 = a; float r23706 = -r23705; float r23707 = fma(r23702, r23704, r23706); float r23708 = hypot(r23698, r23703); float r23709 = r23707 / r23708; float r23710 = 1.5767625849600095e+154; bool r23711 = r23698 <= r23710; float r23712 = 1; float r23713 = r23712 / r23708; float r23714 = r23698 * r23705; float r23715 = fma(r23701, r23703, r23714); float r23716 = r23715 / r23708; float r23717 = r23713 * r23716; float r23718 = fma(r23703, r23702, r23705); float r23719 = r23718 / r23708; float r23720 = r23711 ? r23717 : r23719; float r23721 = r23700 ? r23709 : r23720; return r23721; } double f_od(double a, double b, double c, double d) { double r23722 = c; double r23723 = -2.114794877571543e+145; bool r23724 = r23722 <= r23723; double r23725 = b; double r23726 = r23725 / r23722; double r23727 = d; double r23728 = -r23727; double r23729 = a; double r23730 = -r23729; double r23731 = fma(r23726, r23728, r23730); double r23732 = hypot(r23722, r23727); double r23733 = r23731 / r23732; double r23734 = 1.5767625849600095e+154; bool r23735 = r23722 <= r23734; double r23736 = 1; double r23737 = r23736 / r23732; double r23738 = r23722 * r23729; double r23739 = fma(r23725, r23727, r23738); double r23740 = r23739 / r23732; double r23741 = r23737 * r23740; double r23742 = fma(r23727, r23726, r23729); double r23743 = r23742 / r23732; double r23744 = r23735 ? r23741 : r23743; double r23745 = r23724 ? r23733 : r23744; return r23745; } 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 r23746, r23747, r23748, r23749, r23750, r23751, r23752, r23753, r23754, r23755, r23756; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23746); mpfr_init(r23747); mpfr_init(r23748); mpfr_init(r23749); mpfr_init(r23750); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init(r23754); mpfr_init(r23755); mpfr_init(r23756); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r23746, a, MPFR_RNDN); mpfr_set_d(r23747, c, MPFR_RNDN); mpfr_mul(r23748, r23746, r23747, MPFR_RNDN); mpfr_set_d(r23749, b, MPFR_RNDN); mpfr_set_d(r23750, d, MPFR_RNDN); mpfr_mul(r23751, r23749, r23750, MPFR_RNDN); mpfr_add(r23752, r23748, r23751, MPFR_RNDN); mpfr_mul(r23753, r23747, r23747, MPFR_RNDN); mpfr_mul(r23754, r23750, r23750, MPFR_RNDN); mpfr_add(r23755, r23753, r23754, MPFR_RNDN); mpfr_div(r23756, r23752, r23755, MPFR_RNDN); return mpfr_get_d(r23756, MPFR_RNDN); } static mpfr_t r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778, r23779, r23780; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23757); mpfr_init_set_str(r23758, "-2.114794877571543e+145", 10, MPFR_RNDN); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init(r23768); mpfr_init_set_str(r23769, "1.5767625849600095e+154", 10, MPFR_RNDN); mpfr_init(r23770); mpfr_init_set_str(r23771, "1", 10, MPFR_RNDN); mpfr_init(r23772); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); mpfr_init(r23779); mpfr_init(r23780); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r23757, c, MPFR_RNDN); ; mpfr_set_si(r23759, mpfr_cmp(r23757, r23758) <= 0, MPFR_RNDN); mpfr_set_d(r23760, b, MPFR_RNDN); mpfr_div(r23761, r23760, r23757, MPFR_RNDN); mpfr_set_d(r23762, d, MPFR_RNDN); mpfr_neg(r23763, r23762, MPFR_RNDN); mpfr_set_d(r23764, a, MPFR_RNDN); mpfr_neg(r23765, r23764, MPFR_RNDN); mpfr_fma(r23766, r23761, r23763, r23765, MPFR_RNDN); mpfr_hypot(r23767, r23757, r23762, MPFR_RNDN); mpfr_div(r23768, r23766, r23767, MPFR_RNDN); ; mpfr_set_si(r23770, mpfr_cmp(r23757, r23769) <= 0, MPFR_RNDN); ; mpfr_div(r23772, r23771, r23767, MPFR_RNDN); mpfr_mul(r23773, r23757, r23764, MPFR_RNDN); mpfr_fma(r23774, r23760, r23762, r23773, MPFR_RNDN); mpfr_div(r23775, r23774, r23767, MPFR_RNDN); mpfr_mul(r23776, r23772, r23775, MPFR_RNDN); mpfr_fma(r23777, r23762, r23761, r23764, MPFR_RNDN); mpfr_div(r23778, r23777, r23767, MPFR_RNDN); if (mpfr_get_si(r23770, MPFR_RNDN)) { mpfr_set(r23779, r23776, MPFR_RNDN); } else { mpfr_set(r23779, r23778, MPFR_RNDN); }; if (mpfr_get_si(r23759, MPFR_RNDN)) { mpfr_set(r23780, r23768, MPFR_RNDN); } else { mpfr_set(r23780, r23779, MPFR_RNDN); }; return mpfr_get_d(r23780, MPFR_RNDN); } static mpfr_t r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23781); mpfr_init_set_str(r23782, "-2.114794877571543e+145", 10, MPFR_RNDN); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init(r23791); mpfr_init(r23792); mpfr_init_set_str(r23793, "1.5767625849600095e+154", 10, MPFR_RNDN); mpfr_init(r23794); mpfr_init_set_str(r23795, "1", 10, MPFR_RNDN); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init(r23801); mpfr_init(r23802); mpfr_init(r23803); mpfr_init(r23804); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r23781, c, MPFR_RNDN); ; mpfr_set_si(r23783, mpfr_cmp(r23781, r23782) <= 0, MPFR_RNDN); mpfr_set_d(r23784, b, MPFR_RNDN); mpfr_div(r23785, r23784, r23781, MPFR_RNDN); mpfr_set_d(r23786, d, MPFR_RNDN); mpfr_neg(r23787, r23786, MPFR_RNDN); mpfr_set_d(r23788, a, MPFR_RNDN); mpfr_neg(r23789, r23788, MPFR_RNDN); mpfr_fma(r23790, r23785, r23787, r23789, MPFR_RNDN); mpfr_hypot(r23791, r23781, r23786, MPFR_RNDN); mpfr_div(r23792, r23790, r23791, MPFR_RNDN); ; mpfr_set_si(r23794, mpfr_cmp(r23781, r23793) <= 0, MPFR_RNDN); ; mpfr_div(r23796, r23795, r23791, MPFR_RNDN); mpfr_mul(r23797, r23781, r23788, MPFR_RNDN); mpfr_fma(r23798, r23784, r23786, r23797, MPFR_RNDN); mpfr_div(r23799, r23798, r23791, MPFR_RNDN); mpfr_mul(r23800, r23796, r23799, MPFR_RNDN); mpfr_fma(r23801, r23786, r23785, r23788, MPFR_RNDN); mpfr_div(r23802, r23801, r23791, MPFR_RNDN); if (mpfr_get_si(r23794, MPFR_RNDN)) { mpfr_set(r23803, r23800, MPFR_RNDN); } else { mpfr_set(r23803, r23802, MPFR_RNDN); }; if (mpfr_get_si(r23783, MPFR_RNDN)) { mpfr_set(r23804, r23792, MPFR_RNDN); } else { mpfr_set(r23804, r23803, MPFR_RNDN); }; return mpfr_get_d(r23804, MPFR_RNDN); }