fix norm in CWW

passscoed · passscoed · commit fec07428e4d2 · 2025-11-06T18:41:52.000+01:00
diff --git a/src/pyGroupedTransforms/GroupedCoefficients.py b/src/pyGroupedTransforms/GroupedCoefficients.py
@@ -193,7 +193,7 @@ def norms(self, Dict=False, other=None, m=None):
                                             Psi, scipy.linalg.circulant(variances(j[dd], m))
                                         )
 
-                                n[i] += a @ Psi @ a
+                                n[i] += a.T @ Psi @ a
                                 ac_in = ac_in + 2 ** np.sum(j)
                         return np.sqrt(n)
                     else:
@@ -221,7 +221,7 @@ def norms(self, Dict=False, other=None, m=None):
                                             Psi, scipy.linalg.circulant(variances(j[dd], m))
                                         )
 
-                                n[i] += a @ Psi @ a
+                                n[i] += a.T @ Psi @ a
                                 ac_in = ac_in + 2 ** np.sum(j)
 
                             if len(self.settings[i].u) != 0:

Original file line number	Diff line number	Diff line change
`@@ -193,7 +193,7 @@ def norms(self, Dict=False, other=None, m=None):`
`193`	`193`	`Psi, scipy.linalg.circulant(variances(j[dd], m))`
`194`	`194`	`)`
`195`	`195`
`196`		`- n[i] += a @ Psi @ a`
	`196`	`+ n[i] += a.T @ Psi @ a`
`197`	`197`	`ac_in = ac_in + 2 ** np.sum(j)`
`198`	`198`	`return np.sqrt(n)`
`199`	`199`	`else:`
`@@ -221,7 +221,7 @@ def norms(self, Dict=False, other=None, m=None):`
`221`	`221`	`Psi, scipy.linalg.circulant(variances(j[dd], m))`
`222`	`222`	`)`
`223`	`223`
`224`		`- n[i] += a @ Psi @ a`
	`224`	`+ n[i] += a.T @ Psi @ a`
`225`	`225`	`ac_in = ac_in + 2 ** np.sum(j)`
`226`	`226`
`227`	`227`	`if len(self.settings[i].u) != 0:`