At Mon, 22 Oct 2007 16:16:22 -0400, Benjamin Redelings wrote: > Q1: I'm very much interested in advice on how to handle this kind of > situation. Should I expect this to work in GSL, or is there some reason > why GSL does not want this to work (e.g. speed issues)? Should I write > some code myself to revert to a gaussian when a > 10?
Hello, Thanks for the bug report, I have added a test case for that value. You should expect this sort of thing to work in GSL -- definitely. I added the hard failures in 1.10 as the problem was reported when I was making the release and I didn't have time to fix it fully. Some of the other inverse functions were also affected so I thought it best to report an error. I would welcome any patches or other suggestions on the best way to calculate this. High precision test values to add to test.c are also useful. -- Brian Gough GNU Scientific Library - http://www.gnu.org/software/gsl/
