init; % Initialization (addpath...) seti = struct; seti = setData(seti); % Set data: experimental set-up, contrast, simulate data seti = setRecon(seti); % Settings for the variational reconstruction seti = recon(seti); % Variational reconstruction (process) figure(1); imagesc(real(seti.G(seti.qROIexact))); axis xy; colorbar; %set(gca,'FontSize',20); axis square; caxis([-0.15 +1.42]); print(1,'-depsc','guideRecon1.eps'); figure(2); imagesc(real(seti.G(seti.qROIcomp))); axis xy; colorbar; %set(gca,'FontSize',20); axis square; caxis([-0.15 +1.42]); print(2,'-depsc','guideRecon2.eps'); fprintf('Reconstruction stopped after %i outer iterations...\n', seti.iOutStop); fprintf('... with relative discrepancy %.4f...\n', seti.dis(seti.iOutStop)); fprintf('... and relative error %.4f.\n', seti.err(seti.iOutStop));