In the file eig.c, function jacobi, there are these two lines (265-266):
tol = pow(0.01,(2*iter));
tol = 0.0;
In Numerical Recipes, I find that the lines are embedded in an IF statement.
A tolerance of zero seems a bit overkill.
Also, one of the input parameters is tol.
If indeed wrong, what should it be?