/* Assemble desired bits into floating point variable */
for (y=i=0; i<LD_B1B_DIG; i++) {
- if ((a+i & MASK)==z) x[z=(z+1 & MASK)] = 0;
+ if ((a+i & MASK)==z) x[(z=(z+1 & MASK))-1] = 0;
y = 1000000000.0L * y + x[a+i & MASK];
}
} else {
shlim(f, 0);
}
- return 0;
+ return sign * 0.0;
}
if (!gotrad) rp = dc;
while (dc<8) x *= 16, dc++;