fileo="/Users/francovazza/Desktop/articles/proposals/weak.png"

fileo2="/Users/francovazza/Desktop/articles/proposals/strong.png"

np=[32.0,64.0,128.0,256.0,512.0,1024.0]

grid1=[322,406,512,644,812,1024]

time=[99.29,103.,107.15,134.88,318.,1846.]

tperf=[99.29,99.29,99.29,99.29,99.29]

np2=[5.0,40.,320.]

time2=[50.,38.,67.]

grid2=[864.*0.5,864.,864.*2.]

using PyPlot,Colors

clf()

rc("font", family="arial")

loglog(np[1:5],time[1:5],color="#0f87bf",linestyle="-",marker="o")

plot(np[1:5],tperf[1:5],color="red",linestyle=":")

axis("tight")

ylabel("Wall Clock Time [s]",size=15)

xlabel("# cores",size=15)

title("weak scaling on Jureca",size=18)

axis([20,600,50,1000])

savefig(fileo)

clf()

np=[7.0,14.0,28.0,56.0,112.0]

time=5.*[2.38e1,1.411546e+01,7.623709e+00,4.500006e+00,2.90]/(1.75)

tperfect=[time[1],time[1],time[1],time[1],time[1]]

for i in 1:5

tperfect[i]=tperfect[i]/(np[i]/7)

rc("font", family="arial")

loglog(np,time,color="#0f87bf",linestyle="-",marker="o")

plot(np,tperfect,color="red",linestyle=":")

ylabel("Wall Clock Time [s]",size=15)

xlabel("# cores",size=15)

title("strong scaling on Jureca",size=18)

axis([2,200,4,1e2])

savefig(fileo2)