diff options
author | Nobuyoshi Nakada <nobu@ruby-lang.org> | 2019-07-14 16:31:36 +0900 |
---|---|---|
committer | Nobuyoshi Nakada <nobu@ruby-lang.org> | 2019-07-14 16:56:25 +0900 |
commit | 10de5f149acc620e12be5943e4e29f4f555b7551 (patch) | |
tree | 77a72e57694c9dca3920463f338e4a16ed64dfed | |
parent | 83153bbb149f5c827c4ce537722d3169d24b9de9 (diff) |
Calculate float complex division per each part
Arguments to f_complex_new2 should not be Complex, or violate the
assertion.
-rw-r--r-- | complex.c | 14 |
1 files changed, 4 insertions, 10 deletions
@@ -817,25 +817,19 @@ f_divide(VALUE self, VALUE other, if (f_gt_p(f_abs(bdat->real), f_abs(bdat->imag))) { r = (*func)(bdat->imag, bdat->real); n = f_mul(bdat->real, f_add(ONE, f_mul(r, r))); - if (flo) - return f_complex_new2(CLASS_OF(self), - (*func)(self, n), - (*func)(f_negate(f_mul(self, r)), n)); x = (*func)(f_add(adat->real, f_mul(adat->imag, r)), n); y = (*func)(f_sub(adat->imag, f_mul(adat->real, r)), n); } else { r = (*func)(bdat->real, bdat->imag); n = f_mul(bdat->imag, f_add(ONE, f_mul(r, r))); - if (flo) - return f_complex_new2(CLASS_OF(self), - (*func)(f_mul(self, r), n), - (*func)(f_negate(self), n)); x = (*func)(f_add(f_mul(adat->real, r), adat->imag), n); y = (*func)(f_sub(f_mul(adat->imag, r), adat->real), n); } - x = rb_rational_canonicalize(x); - y = rb_rational_canonicalize(y); + if (!flo) { + x = rb_rational_canonicalize(x); + y = rb_rational_canonicalize(y); + } return f_complex_new2(CLASS_OF(self), x, y); } if (k_numeric_p(other) && f_real_p(other)) { |