Fix eigenvalue solver for complex matrices

2017-09-24 19:11:56 -07:00 · 2017-09-24 19:11:56 -07:00 · c14298484d
commit c14298484d
parent 6748181f8f
1 changed files with 1 additions and 1 deletions
--- a/fdfd_tools/waveguide_mode.py
+++ b/fdfd_tools/waveguide_mode.py
@ -64,7 +64,7 @@ def solve_waveguide_mode_2d(mode_number: int,

    eigval = None
    for _ in range(40):
-        eigval = v @ A @ v
+        eigval = v.conj() @ A @ v
        if numpy.linalg.norm(A @ v - eigval * v) < 1e-13:
            break
        w = spalg.spsolve(A - eigval * sparse.eye(A.shape[0]), v)