sscalling.m
a3531b44407cgeom_test

sscalling.m
View Options

			%% Strong scaling for a collisional system with 10M electrons using FENNECS
			% the geometry is a coaxial configuration of constant radius

			Ncores = [ 1 2 4 8 16];% 36 72];
			time = [ (10890+10150)/2 6402 3799 2366 2128];% 1665 1155];

			f=figure;
			plot(Ncores, time(1)*Ncores(1)./time,'r-','displayname','Measured','linewidth',1.5)
			hold on
			plot(Ncores,Ncores,'k--','displayname','Ideal','linewidth',1.5)
			ylabel('t_1/t_N')
			yyaxis right
			plot(Ncores,time(1).Ncores(1)./(time.Ncores),'bx:','linewidth',1.5,'markersize',12,'displayname','Efficiency')
			xlabel('N_{cores}')
			ylabel('Efficiency')
			set(gca,'fontsize',12)
			grid on
			legend('location','southeast')
			ax = gca;
			ax.YAxis(1).Color = 'k';
			ax.YAxis(2).Color = 'b';
			ylim([0 1])
			M.savegraph(f,'sscalingOMP');

			%% MPI with 4 cores per task

			Ntasks = [ 1 2 4 8];
			time = [3799 1533 976 654];

			f=figure;
			plot(Ntasks, time(1)./time,'r-','displayname','Measured','linewidth',1.5)
			hold on
			ylabel('Speedup [-]')
			plot(Ntasks,Ntasks,'k--','displayname','Ideal','linewidth',1.5)
			yyaxis right
			plot(Ntasks,time(1).Ntasks(1)./(time.Ntasks),'bx:','linewidth',1.5,'markersize',12,'displayname','Efficiency')
			xlabel('N_{tasks} [-]')
			ylabel('Parallel efficiency [-]')
			set(gca,'fontsize',12)
			grid on
			legend('location','southeast')
			ylim([0 2])
			xlim([0 9])
			ax = gca;
			ax.YAxis(1).Color = 'k';
			ax.YAxis(2).Color = 'b';
			M.savegraph(f,'sscalingMPI');