put debug prints in calc_cartesian_dcw

SyneRBI · Jun 19, 2024 · df80044 · df80044
1 parent 5849c07
commit df80044
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/xGadgetron/pGadgetron/Gadgetron.py b/src/xGadgetron/pGadgetron/Gadgetron.py
@@ -1855,12 +1855,18 @@ def calc_cartesian_dcw(ad):
     traj, inverse, counts = numpy.unique(traj, return_inverse=True, return_counts=True, axis=1)
 
     density_weight = (1.0 / counts)[inverse]
+    print(density_weight.shape)
 
     density_weight = numpy.expand_dims(density_weight, axis=1)
+    print(density_weight.shape)
     density_weight = numpy.expand_dims(density_weight, axis=2)
+    print(density_weight.shape)
+    print(ad.shape)
     density_weight = numpy.tile(density_weight, (1, ad.shape[1], ad.shape[2]))
+    print(density_weight.shape)
 
     dcw = ad.copy()
+    print(dcw.shape)
     dcw.fill(density_weight)
 
     return dcw