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;