diff --git a/dedalus/core/solvers.py b/dedalus/core/solvers.py index c9eb4c94..b333b794 100644 --- 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