From 078b9e2db2e25d663a1140cc71ee4d16c36cc161 Mon Sep 17 00:00:00 2001 From: "Jakob Jorgensen, algol@harwell" Date: Thu, 7 Sep 2017 15:09:50 +0100 Subject: bugsf --- main_func/FISTA_REC.m | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/main_func/FISTA_REC.m b/main_func/FISTA_REC.m index 8dd569f..00e59ab 100644 --- a/main_func/FISTA_REC.m +++ b/main_func/FISTA_REC.m @@ -104,7 +104,7 @@ else astra_mex_data3d('delete', sino_id); astra_mex_data3d('delete', id); end - clear proj_geomT vol_geomT + %clear proj_geomT vol_geomT else % divergen beam geometry fprintf('%s \n', 'Calculating Lipshitz constant for divergen beam geometry... will take some time!'); @@ -309,7 +309,7 @@ if (subsets == 0) if (strcmp(proj_geom.type,'parallel') || strcmp(proj_geom.type,'parallel3d')) sino_updt = zeros(size(sino),'single'); for kkk = 1:SlicesZ - [sino_id, sino_updt(:,kkk,:)] = astra_create_sino3d_cuda(X_t(:,:,kkk), proj_geomT, vol_geomT); + [sino_id, sino_updt(:,:,kkk)] = astra_create_sino3d_cuda(X_t(:,:,kkk), proj_geomT, vol_geomT); astra_mex_data3d('delete', sino_id); end else @@ -338,7 +338,7 @@ if (subsets == 0) if (strcmp(proj_geom.type,'parallel') || strcmp(proj_geom.type,'parallel3d')) x_temp = zeros(size(X),'single'); for kkk = 1:SlicesZ - [id, x_temp(:,:,kkk)] = astra_create_backprojection3d_cuda(squeeze(residual(:,kkk,:)), proj_geomT, vol_geomT); + [id, x_temp(:,:,kkk)] = astra_create_backprojection3d_cuda(squeeze(residual(:,:,kkk)), proj_geomT, vol_geomT); astra_mex_data3d('delete', id); end else -- cgit v1.2.3