diff --git a/MATLAB/Algorithms/PICCS.m b/MATLAB/Algorithms/PICCS.m index 792de229..7b6e04f5 100644 --- 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('');