Nomenclature and figure prettynesss

ComputeMostHexagonalSampling.m

@@ -122,6 +122,7 @@ d_min = zeros(size(LArray));
 
 if doDebug
     figure('Name','Sampling grid and alps');
+    set(gcf, 'Position', [680   693   621   285])
     ah = subplot(1,2,1);
 end
 
@@ -140,7 +141,8 @@ for iL = 1:numel(LArray)
     if doDebug
         [kyArray, kzArray] = GetGridPoints(3*Lmax, deltak_inplane, sliceSpacing, LArray(iL));
         PlotGrid(kyArray, kzArray, sliceSpacing, ah);
-
+        title(sprintf('L = %.2f      %s = %.2f',LArray(iL), char(hex2dec('0394')), nSlicesSimultaneously/LArray(iL)));
+        
         subplot(1,2,2); hold on;
         plot(LArray(iL),packingDensityArray(iL),'.','Color','k');
 
@@ -151,8 +153,8 @@ for iL = 1:numel(LArray)
         end
         
         xlabel('L'); ylabel('Packing density [%]');
-        xlim([Lmin Lmax]); ylim([min(packingDensityArray) 90.7]);
-        sgtitle(sprintf('L = %.2f',LArray(iL)));
+        xlim([Lmin Lmax]); ylim([0 90.7]);
+        axis square;
 
         %% Make GIF
         if ~isempty(gifName)
@@ -213,7 +215,7 @@ if doPlot
     plot(L_CAIPI_B, packingDensityCAIPI_B, 'x', 'MarkerSize', 10, 'LineWidth', 3, ...
         'Color', GetSamplingPatternColour(2), 'DisplayName', 'Blipped-CAIPIRINHA - B');
     plot(L_optimum, packingDensity, 'x', 'MarkerSize', 10, 'LineWidth', 3, ...
-        'Color', GetSamplingPatternColour(0), 'DisplayName', 'Optimum sampling');
+        'Color', GetSamplingPatternColour(0), 'DisplayName', 'Proposed method');
     legend;
     xlabel('L');
     ylabel('Packing density [%]');
@@ -240,6 +242,9 @@ scatter(ky, kz, 'MarkerEdgeColor', [0.5 0.5 0.5], 'MarkerFaceColor', [0.5 0.5 0.
 
 hold on;
 
+ind0 = find(ky==0);
+scatter(ky(ind0:ind0+1), kz(ind0:ind0+1), 'MarkerEdgeColor', [0 0 0], 'MarkerFaceColor', [0 0 0], 'LineWidth', 5);
+
 slabBorderColour = [100 149 237]/255; % light blue
 yline(0,'r','LineWidth',3, 'Color', slabBorderColour); 
 yline(2*pi/sliceSpacing, 'r', 'LineWidth', 3, 'Color', slabBorderColour);

MakeGridAndAlpGif.m

@@ -5,6 +5,8 @@ baseFolder = pwd;
 doSave = true;
 if doSave
     gifName = fullfile(baseFolder,'GridAndAlpsGif.gif');
+    dfws = get(0,'DefaultFigureWindowStyle');
+    set(0,'DefaultFigureWindowStyle','normal');
 else
     gifName = [];
 end
@@ -18,5 +20,6 @@ sliceSpacing = FOVinplane/2/Nsim;
 [packingDensity, L_optimum, isDegenerate, packingDensityCAIPItmp, L_CAIPItmp] = ...
     ComputeMostHexagonalSampling(R, FOVinplane, sliceSpacing, Nsim, [], [], [], [], true, gifName, 0.02);
 
-[packingDensityCAIPI, indCAIPI] = max(packingDensityCAIPItmp);
-L_CAIPI = L_CAIPItmp(indCAIPI);
+if doSave
+    set(0,'DefaultFigureWindowStyle',dfws);
+end

MakeOptDvsRPlot.m

@@ -52,7 +52,7 @@ for iNSim = 1:numel(NSimArray)
         end
         nexttile;  hold on;
         plot(RArray(1:numel(packingDensity))*Nsim, packingDensity, 'LineWidth', 2,...
-            'DisplayName', ['As hexagonal as possible' name], 'Color', GetSamplingPatternColour(0));
+            'DisplayName', ['Proposed method' name], 'Color', GetSamplingPatternColour(0));
         plot(RArray(1:numel(packingDensity))*Nsim, packingDensityCAIPI_A, 'LineWidth', 2,...
             'DisplayName', ['Blipped-CAIPIRINHA - A' name], 'Color', GetSamplingPatternColour(1));
         plot(RArray(1:min(numel(packingDensity),numel(packingDensityCAIPI_B)))*Nsim, packingDensityCAIPI_B, 'LineWidth', 2,...