comp = [
%8x8       16x16       32x32      64x64      128x128      256x256
3.89531   15.038790   67.609310  293.333267 1316.030000  6260.702702 % 1proc
NaN       6.169225    25.938822  104.177487 445.498488   2145.335218 % 4proc
NaN       NaN         14.092557  52.729437  207.071409   929.346426  % 16proc
NaN       NaN         NaN        25.491762  93.926878    388.967145  % 64proc
];


% the array of T1(alpha x p)
t1_alphap = [
% Alpha (particles/proc)
% 64          256          1024         4096        16K          64K
3.89531       15.038790    67.609310    293.333267  1316.030000  6260.702702
15.038790     67.609310    293.333267 	1316.030000 6260.702702  NaN
67.609310     293.333267   1316.030000	6260.702702 NaN          NaN
293.333267    1316.030000  6260.702702  NaN         NaN          NaN
];

% The array of Tp(alpha x p)
tp_alphap = [
% Alpha (particles/proc)
% 64          256          1024         4096        16K          64K
3.89531    15.038790  67.609310  293.333267 1316.030000  6260.702702
6.169225   25.938822  104.177487 445.498488 2145.335218  NaN
14.092557  52.729437  207.071409 929.346426 NaN          NaN
25.491762  93.926878  388.967145 NaN        NaN          NaN
];

% scaled speedup
scspeedup = t1_alphap ./ tp_alphap;
% particles/proc
alpha = [64 256 1024];

% To plot scaled speedup:
figure(1);
meshz(procs,alpha,scspeedup)
xlabel('Num Procs');
ylabel('alpha (particles/proc)');
zlabel('Scaled Speedup');
title('Scaled Speedup');

n = [8 16 32 64 128 256];
procs = [1 4 16 64];

speedup = ones(4,6);
speedup(2,1:6) = comp(1,:)./comp(2,1:6);
speedup(3,1:6) = comp(1,1:6)./comp(3,1:6);
speedup(4,1:6) = comp(1,1:6)./comp(4,1:6);

figure(2);
meshz(procs,n,speedup)
xlabel('Num Procs');
ylabel('n (nxn world)');
zlabel('Scaled Speedup');
title('Scaled Speedup');