Fix _normalize_left_eigenvectors in EigenvalueSolver by dividing by t…

…he conjugate of the norm.
DedalusProject · May 17, 2023 · 08d20ef · 08d20ef
1 parent 735a666
commit 08d20ef
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/dedalus/core/solvers.py b/dedalus/core/solvers.py
@@ -111,7 +111,7 @@ def solve_dense(self, pencil, rebuild_coeffs=False, normalize_left=False, **kw):
         self.eigenvalue_pencil = pencil
         if len(eig_output) == 3 and normalize_left:
             norms = np.diag(self.modified_left_eigenvectors.T.conj() @ self.eigenvectors)
-            self.left_eigenvectors /= norms
+            self.left_eigenvectors /= np.conj(norms)
             self.modified_left_eigenvectors = np.conjugate(np.transpose(np.conjugate(self.left_eigenvectors.T) * -pencil.M))
 
     def solve_sparse(self, pencil, N, target, rebuild_coeffs=False, left=False, normalize_left=True, raise_on_mismatch=True, **kw):