minor cosmetics

Author: pjmark

niftypet/nipet/sct/mmrsct.py

@@ -476,7 +476,7 @@ def vsm(
         muim = muo + muh
         emim = em
 
-    muim = ndi.interpolation.zoom(muim, Cnt['SCTSCLMU'], order=3) # (0.499, 0.5, 0.5)
+    muim = ndi.interpolation.zoom(muim, Cnt['SCTSCLMU'], order=3) # (0.499, 0.5, 0.5) <<<<<<<<<<<< TODO: deprecation
     emim = ndi.interpolation.zoom(emim, Cnt['SCTSCLEM'], order=3) # (0.34, 0.33, 0.33)
 
     # -smooth the mu-map for mask creation.
@@ -627,6 +627,7 @@ def vsm(
 
     # ======= SCALING ========
     # > scale scatter using non-TOF SSRB sinograms
+    #import pdb; pdb.set_trace()
 
     # > gap mask
     rmsk = (txLUT['msino'] > 0).T

niftypet/nipet/sct/src/sctaux.cu

@@ -239,7 +239,6 @@ iMSK get_imskEm(IMflt imvol, float thrshld, Cnst Cnt) {
   HANDLE_ERROR(cudaFreeHost(h_v2i));
 
 #else
-  // printf(">>>>> NVX:%d, THRESHOLD:%f\n", nvx, thrshld);
   HANDLE_ERROR(cudaMallocManaged(&d_i2v, nvx * sizeof(int)));
   HANDLE_ERROR(cudaMallocManaged(&d_v2i, SSE_IMX * SSE_IMY * SSE_IMZ * sizeof(int)));
 
@@ -258,7 +257,7 @@ iMSK get_imskEm(IMflt imvol, float thrshld, Cnst Cnt) {
 #endif
 
   if (Cnt.LOG <= LOGDEBUG)
-    printf("d> number of voxel values greater than %3.2f is %d out of %d (ratio: %3.2f)\n",
+    printf("d> number of voxel values greater than %3.2f (THRSHL=%3.2f) is %d out of %d (ratio: %3.2f)\n",
            thrshld, nvx, SSE_IMX * SSE_IMY * SSE_IMZ, nvx / (float)(SSE_IMX * SSE_IMY * SSE_IMZ));
   msk.nvx = nvx;
   msk.i2v = d_i2v;