From 3b92b642eb3f7e3ed1b68bda86d073d310c5d131 Mon Sep 17 00:00:00 2001 From: "Yukihiro \"Matz\" Matsumoto" Date: Thu, 28 Mar 2019 14:02:30 +0900 Subject: Break loop whem `sum==0` to avoid zero division; ref #4345 --- mrbgems/mruby-math/src/math.c | 1 + 1 file changed, 1 insertion(+) (limited to 'mrbgems/mruby-math/src/math.c') diff --git a/mrbgems/mruby-math/src/math.c b/mrbgems/mruby-math/src/math.c index 39fed6cd3..ef0dc1857 100644 --- a/mrbgems/mruby-math/src/math.c +++ b/mrbgems/mruby-math/src/math.c @@ -120,6 +120,7 @@ erf(double x) term *= xsqr/j; sum += term/(2*j+1); ++j; + if (sum == 0) break; } while (fabs(term/sum) > DBL_EPSILON); return two_sqrtpi*sum; } -- cgit v1.2.3