minor bug in PICCS

CERN · Jan 15, 2025 · 872134c · 872134c
1 parent f566283
commit 872134c
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/MATLAB/Algorithms/PICCS.m b/MATLAB/Algorithms/PICCS.m
@@ -135,7 +135,7 @@
         %         weigth_backprj=bsxfun(@times,1./V(:,:,jj),backprj); %                 V * At * W^-1 * (b-Ax)
         %         f=f+beta*weigth_backprj;                          % x= x + lambda * V * At * W^-1 * (b-Ax)
         % Enforce positivity
-        f=f+beta* bsxfun(@times,1./V(:,:,index_angles(:,jj)),Atb(W(:,:,index_angles(:,jj)).*(proj(:,:,index_angles(:,jj))-Ax(f,geo,angles_reorder(:,jj),'gpuids',gpuids)),geo,angles_reorder(:,jj),'gpuids',gpuids));
+        f= f+beta* bsxfun(@times,1./V(:,:,jj),Atb(W(:,:,index_angles(:,jj)).*(proj(:,:,index_angles(:,jj))-Ax(f,geo,angles_reorder(:,jj),'gpuids',gpuids)),geo,angles_reorder(:,jj),'gpuids',gpuids));
         % non-negativity constrain
         if nonneg
             f=max(f,0);
@@ -199,7 +199,7 @@
     end
     if (iter==1 && verbose==1)
         expected_time=toc*maxiter;
-        disp('ADS-POCS');
+        disp('PICCS');
         disp(['Expected duration  :    ',secs2hms(expected_time)]);
         disp(['Exected finish time:    ',datestr(datetime('now')+seconds(expected_time))]);
         disp('');