doPlotRes = true;

%% set up the solver

% number of GMRES iterations
iterations = 3;

% As we want to control GMRES with number of iterations rather than
% the tolerance, we use a strict tolerance that we most probably would not
% reach in a few iterations 
tolerance = 1e-16;

%% Run the two cases
for casenr = 1:2
    
    switch casenr
        case 1 % coarse fractures constraints the coarse dual grid
            coarseGrid = cell(1,1);
            coarseGrid{1} = fset{1};
            coarseGrid{1}.merge 