Fix FDK option usegpufft option not working

CERN · May 29, 2023 · 7c97471 · 7c97471
1 parent aed4379
commit 7c97471
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 6 deletions.
diff --git a/MATLAB/Algorithms/FDK.m b/MATLAB/Algorithms/FDK.m
@@ -71,7 +71,7 @@
     proj(:,:,ii) = w.*proj(:,:,ii);
 end
 %% Fourier transform based filtering
-proj = filtering(proj,geo,angles,parker); % Not sure if offsets are good in here
+proj = filtering(proj,geo,angles,parker, 'usegpufft', usegpufft); % Not sure if offsets are good in here
 
 % RMFIELD Remove fields from a structure array.
 geo=rmfield(geo,'filter');

diff --git a/MATLAB/Utilities/filtering.m b/MATLAB/Utilities/filtering.m
@@ -68,11 +68,13 @@
         fproj = (zeros(filt_len,geo.nDetector(2),'single'));
 
         fproj(round(filt_len/2-geo.nDetector(1)/2+1):round(filt_len/2+geo.nDetector(1)/2),:) = proj(:,:,ii);
-
-        fproj = fft(fproj);   
-        fproj = fproj.*filt;
-        fproj = (real(ifft(fproj)));
-
+        if usegpufft==1
+            fproj = ApplyFbpFiltration(fproj, filt, gpuids);
+        else
+            fproj = fft(fproj);   
+            fproj = fproj.*filt;
+            fproj = (real(ifft(fproj)));
+        end
         if parker
             proj(:,:,ii) = fproj(round(end/2-geo.nDetector(1)/2+1):round(end/2+geo.nDetector(1)/2),:)/2/geo.dDetector(1)*(2*pi/  (pi/angle_step)   )/2*(geo.DSD(ii)/geo.DSO(ii));
         else