kotovalexarian-likes-github/ruby--ruby
1
0
Fork 0
mirror of https://github.com/ruby/ruby.git synced 2022-11-09 12:17:21 -05:00
Code Releases Activity

Calculate float complex division per each part

Browse source 
Arguments to f_complex_new2 should not be Complex, or violate the
assertion.
This commit is contained in:
Nobuyoshi Nakada 2019-07-14 16:31:36 +09:00
parent 83153bbb14
commit 10de5f149a
No known key found for this signature in database
GPG key ID: 4BC7D6DF58D8DF60
1 changed files with 4 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
complex.c
View file

@ -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)) {