Normalise left eigenvectors in solve sparse with the conjugate of norms.

DedalusProject · Aug 2, 2023 · 21471db · 21471db
1 parent 08d20ef
commit 21471db
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/dedalus/core/solvers.py b/dedalus/core/solvers.py
@@ -183,7 +183,7 @@ def solve_sparse(self, pencil, N, target, rebuild_coeffs=False, left=False, norm
             if normalize_left:
                 unnormalized_modified_left_eigenvectors = np.conjugate(np.transpose(np.conjugate(self.left_eigenvectors.T) * -pencil.M))
                 norms = np.diag(unnormalized_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))
         self.eigenvalue_pencil = pencil