GNOME Bugzilla – Bug 697293
=r.psnorm(-0.1,1111,0,1) < 0
Last modified: 2013-04-07 20:58:12 UTC
=r.psnorm(-0.1,1111,0,1) is not supposed to be -0.01627980883793 And that's a little far for a rounding error.
It's really funky there: -0.1 ---> -0.016280 0.0 ---> 0.437235 0.1 ---> 0.063376 I.e., it is not even close to monotonic.
This problem has been fixed in our software repository. The fix will go into the next software release. Thank you for your bug report.
Created attachment 240797 [details] sample Gnumeric file I see some curious blips in B17 and B29. Obviously they are well within rounding problem size. I just find it curious that they show up when the surrounding values are listed as 0.
There are 18 different parameter areas covering six different evaluation methods, some used with different parameters. Minor bumps when moving from one area to another are to be expected. I don't think the accuracy claim is better than 10ulp. I'd have to look to be sure.
On second thought, what we have here is probably a case of catastrophic cancellation in psnorm. The reason you see zeros is probably because I added a clamp to [0;1].
Should be fixed now.