Filter large but finite eigenvalues appearing with libopenblas 0.3.21

DedalusProject · Aug 12, 2022 · 8f6d212 · 8f6d212
1 parent 90b7b7f
commit 8f6d212
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/dedalus/tests/test_evp.py b/dedalus/tests/test_evp.py
@@ -217,7 +217,7 @@ def test_half_qho_dense_evp(benchmark, x_basis_class, Nx, dtype):
     solver = problem.build_solver()
     solver.solve_dense(solver.pencils[0])
     # Filter infinite/nan eigenmodes
-    finite = np.isfinite(solver.eigenvalues)
+    finite = np.abs(solver.eigenvalues) < 1e10
     solver.eigenvalues = solver.eigenvalues[finite]
     solver.eigenvectors = solver.eigenvectors[:, finite]
     # Sort eigenmodes by eigenvalue