init; % Initialization (addpath...) seti.useTolOut = 1; % Use outer tolerance principle to stop the inner iteration 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) imagesc(real(seti.G(seti.qROIcomp))); axis xy; colorbar;