Menu

#1424 tolerance in jacobi confusing

open
nobody
None
5
2017-07-30
2017-07-30
Erik
No

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?

Discussion


Log in to post a comment.

MongoDB Logo MongoDB