From a24e6084132b9131b7a5179b98fd1b02835d894c Mon Sep 17 00:00:00 2001 From: Fabian Joswig Date: Mon, 25 Oct 2021 14:26:58 +0100 Subject: [PATCH] g5H operation of Npr_matrix optimized --- pyerrors/npr.py | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/pyerrors/npr.py b/pyerrors/npr.py index a8b5ab9a..acd21dff 100644 --- a/pyerrors/npr.py +++ b/pyerrors/npr.py @@ -26,10 +26,9 @@ class Npr_matrix(np.ndarray): if self.shape != (12, 12): raise Exception('g5H only works for 12x12 matrices.') extended_g5 = np.kron(np.eye(3, dtype=int), gamma5) - new_matrix = extended_g5 @ self.conj().T @ extended_g5 - new_matrix.mom_in = self.mom_out - new_matrix.mom_out = self.mom_in - return new_matrix + return Npr_matrix(matmul(extended_g5, self.conj().T, extended_g5), + mom_in=self.mom_out, + mom_out=self.mom_in) def _propagate_mom(self, other, name): s_mom = getattr(self, name, None)