    % ticGuide = tic;
    init;                   % Initialization (addpath...)

    seti.dim = 3;
    seti.rCD = 2.0;
    seti.k = 10;
    seti.contrast = 'twoTripods3D';
    seti.incNb = 35;
    seti.measNb = 35;
    seti.radSrc = 5;
    seti.radMeas = 5;
    seti.tau = 1.25;

    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)
    % tocGuide = toc(ticGuide)

    figure(1); contourPlotROI(seti.qROIexact, seti, 'real');
    % set(gca,'FontSize',20); print(1,'-depsc','guideRecon3D1.eps');

    figure(2); contourPlotROI(seti.qROIcomp, seti, 'real');
    % set(gca,'FontSize',20); print(2,'-depsc','guideRecon3D2.eps');
    % set(gca,'FontSize',20); print(2,'-dpng','guideRecon3D2dpiDefault.png');
    % set(gca,'FontSize',20); print(2,'-dpng','-r600','guideRecon3D2dpi600.png');

    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));

    % -- Result --
    % >> init; guideRecon3D
    %
    % tocGuide =
    %
    %    1.4730e+04 % 4.0916 h
    %
    % Reconstruction stopped after 9 outer iterations...
    % ... with relative discrepancy 0.0125...
    % ... and relative error 0.6590.
    % >>