%% 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('t_1/t_N')
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('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');