Fix _normalize_left_eigenvectors in EigenvalueSolver by dividing by t…

…he conjugate of the norm.
DedalusProject · May 17, 2023 · 19a8d7c · 19a8d7c
1 parent 29c96fe
commit 19a8d7c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/dedalus/core/solvers.py b/dedalus/core/solvers.py
@@ -160,7 +160,7 @@ def _build_modified_left_eigenvectors(self):
     def _normalize_left_eigenvectors(self):
         modified_left_eigenvectors = self._build_modified_left_eigenvectors()
         norms = np.diag(modified_left_eigenvectors.T.conj() @ self.eigenvectors)
-        self.left_eigenvectors /= norms
+        self.left_eigenvectors /= np.conj(norms)
 
     def solve_dense(self, subproblem, rebuild_matrices=False, left=False, normalize_left=True, **kw):
         """