change mlx to m file

dev/MEEM/compare_MEEM_figs.m

@@ -0,0 +1,79 @@
+%% Comparing using small and large
+
+folder = ['dev' filesep 'MEEM_figs'];
+
+phis = {'small-small','large-large'};
+vels = {'small-small','large-small','large-large'};
+figure 
+tcl = tiledlayout(length(vels),length(phis)*2);
+for i_phi = 1:length(phis)
+    for i_vel = 1:length(vels)
+        file = sprintf('*%s*.fig',['phi = ' phis{i_phi} ', vel = ' vels{i_vel}]);
+        all_fnames = dir(fullfile(folder,file));
+        fname = all_fnames(end).name; % latest
+        fig = openfig(fname);
+
+        ax = fig.Children(4);
+        ax.Parent = tcl;
+        ax.Layout.Tile = (i_vel-1)*2*length(phis) + 2*i_phi - 1;
+        if i_vel==length(vels)
+            xl = xlim;
+            yl = ylim;
+            text(ax, 1.25*xl(2), 1.25*yl(1), phis{i_phi})
+        end
+        if i_phi == 1
+            ylabel(ax, {vels{i_vel}, 'Z'})
+        end
+        title(ax,'Real')
+
+        ax2 = fig.Children(2);
+        ax2.Parent = tcl;
+        ax2.Layout.Tile = (i_vel-1)*2*length(phis) + 2*i_phi;
+        title(ax2,'Imaginary')
+    end
+end
+title(tcl,'Total Potential for Different Matching Strategies (N=6)')
+subtitle(tcl,'Integration Region Size - Eigenfunction Size')
+xlabel(tcl,'Potential Matching Strategy')
+ylabel(tcl,'Velocity Matching Strategy')
+
+%% Generate binary combinations
+clear all
+for num=0:31
+    disp(['running num=' num2str(num)])
+    try
+        Chau_Yeung_live()
+    catch
+        disp(['num=' num2str(num) ' failed.'])
+    end
+end
+
+%% Comparing using binary combinations
+folder = ['dev' filesep 'MEEM_figs'];
+
+figure 
+
+file = sprintf('*%s*%s*%s*.fig','04-30',', BC = ','VelMatch');
+all_fnames = dir(fullfile(folder,file));
+tile_length = ceil(sqrt(length(all_fnames)));
+tcl = tiledlayout(tile_length,tile_length);
+
+for i = 1:length(all_fnames)   
+    fname = all_fnames(i).name;
+    BC_num = regexp(fname, 'BC = (\d+)', 'tokens', 'once');
+    fig = openfig(fname);
+    if i~=length(all_fnames)
+        legend off
+        ax = fig.Children;
+    else
+        ax = fig.Children(2);
+    end
+    ax.Parent = tcl;
+    ax.Layout.Tile = i;
+    title(ax,['BC = ', BC_num{:}])
+
+
+end
+
+title(tcl,'Radial Velocity for Different Matching Strategies (N=2)')
+