(minor) bug in KSdist::pks (fixed in V3.04)


Bill Press
08-17-2010, 03:47 PM
The line in pks
if (z == 0.) return 0.;
should be replaced by
if (z < 0.042) return 0.;
This fixes a potential underflow in the exponential, two lines later. Thanks to user tbegg for pointing this out.

emilpohl
09-06-2010, 04:10 PM
Hi,

KSdist::invqks warning "The local variable yp is never read".

Is yp necessary?

Bill Press
09-07-2010, 09:26 AM
Yes, you're right. It's a debug fossil. Remove the declaration of yp, and the line
yp = y;
Thanks for pointing this out.