Commit 82bfe40f authored by Pavel Kus's avatar Pavel Kus

test all eigenvalues in test analytic

even if only part of eigenvectors are calculated, all eigenvalues should
be checked
parent e2d9b343
......@@ -115,10 +115,12 @@ module test_analytic
max_z_diff = ZERO
max_ev_diff = ZERO
do globJ = 1, nev
do globJ = 1, na
diff = abs(ev(globJ) - analytic_eigenvalues(na, globJ))
max_ev_diff = max(diff, max_ev_diff)
end do
do globJ = 1, nev
! calculated eigenvector can be in opposite direction
max_curr_z_diff_minus = ZERO
max_curr_z_diff_plus = ZERO
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment