Even better Williamson

XanaduAI · Jul 17, 2023 · 9544ff9 · 9544ff9
1 parent 1579a7d
commit 9544ff9
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/thewalrus/decompositions.py b/thewalrus/decompositions.py
@@ -107,8 +107,8 @@ def symplectic_eigenvals(cov):
         (array): symplectic eigenvalues
     """
     M = int(len(cov) / 2)
-    D, _ = williamson(cov)
-    return np.diag(D)[:M]
+    Omega = sympmat(M)
+    return np.real_if_close(-1j * np.linalg.eigvals(Omega @ cov))[::2]
 
 
 def blochmessiah(S):