Hello, I believe that intrinsic to the algorithm in gsl_fcmp is that it is not a reasonable thing to write something like gsl_fcmp(1e-16,0.0,1e-13) and expect to receive a 0 indicating that 1e-16 is approximately 0.0.
(|u-v|<=|u|*e and |u-v|<=|v|*e )cannot be true if either u is 0.0 and e and v are not. I was hung up for a moment on this point and think that the gsl documentation for gsl_fcmp should address this point so that novice gsl programmers don't make the same mistake. If I am still confused please correct me, Karl Edler _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
