We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 1d670b1 commit 0bf5b2fCopy full SHA for 0bf5b2f
src/nmath/bd0.c
@@ -52,7 +52,7 @@ attribute_hidden double bd0(double x, double np)
52
if (fabs(x-np) < 0.1*(x+np)) {
53
double d = x - np,
54
v = d/(x+np);
55
- if(d && !v) { // v has underflown to 0 (as x+np = inf)
+ if((d != 0.0) && (v == 0.0)) { // v has underflown to 0 (as x+np = inf)
56
double
57
x_ = ldexp(x, -2),
58
n_ = ldexp(np,-2);
0 commit comments