diff options
author | nobu <nobu@b2dd03c8-39d4-4d8f-98ff-823fe69b080e> | 2018-06-06 14:01:43 +0000 |
---|---|---|
committer | nobu <nobu@b2dd03c8-39d4-4d8f-98ff-823fe69b080e> | 2018-06-06 14:01:43 +0000 |
commit | c995315153f7403636ec0d6119b7885668b32039 (patch) | |
tree | f73f0aee5ae6ff2b100f76df7f2049014bf2bf98 /math.c | |
parent | 439e50789b2b49676c5108894d324c5d1ccebf1f (diff) |
math.c: adjust cbrt
* math.c (math_cbrt): refine the approximation result on boundary
values by an iteration of Newton-Raphson method.
git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@63592 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
Diffstat (limited to 'math.c')
-rw-r--r-- | math.c | 9 |
1 files changed, 8 insertions, 1 deletions
@@ -689,7 +689,14 @@ rb_math_sqrt(VALUE x) static VALUE math_cbrt(VALUE unused_obj, VALUE x) { - return DBL2NUM(cbrt(Get_Double(x))); + double f = Get_Double(x); + double r = cbrt(f); +#if defined __GLIBC__ + if (isfinite(r)) { + r = (2.0 * r + (f / r / r)) / 3.0; + } +#endif + return DBL2NUM(r); } /* |