Homec4science
Diffusion Project_multiproc (ancarola) Browse

multiproc/A1/graphs/gengraphs.py
f679bce86ec6ancarola

Last Change
Raw File

gengraphs.py
View Options

#!/usr/bin/python
import matplotlib.pyplot as plt
import numpy as np
from numpy import array
import tikzplotlib
N = 100000000
beta = 50000
plt.style.use("ggplot")
# 20 -> 0.95
# 30 -> 0.966666
# 50 -> 0.98
def speedup(n):
p = 1.0 / (1 + (beta * (n-1)) / float(N))
return 1.0 / (1 - p + p / n)
x = array([1, 2, 4, 8, 16, 32, 48, 64])
y = speedup(x)
def speedup_T(Tref, T):
return Tref / T
# load data
def loadvalue(line):
values = line.split(' ') # split into string set
values = [i for i in values if i != ''] # fliter null strings
pyarr = []
for v in values:
pyarr.append(float(v))
return (np.mean(pyarr), np.var(pyarr))
Tpi = []
Tpiv = []
timespi = open("timespi")
for line in timespi:
(t,v) = loadvalue(line)
Tpi.append(t)
Tpiv.append(v)
timespi.close()
Tint = []
Tintv = []
timesint = open("timesint")
for line in timesint:
(t,v) = loadvalue(line)
Tint.append(t)
Tintv.append(v)
timesint.close()
# good data
#Tpi = array([2.131, 1.067, 0.5339, 0.2673, 0.1362, 0.1051, 0.1083, 0.1089])
#Tint = array([1.064, 0.5328, 0.2672, 0.1333, 0.06728, 0.06198, 0.05691, 0.05873])
Tpi = array(Tpi)
Tint = array(Tint)
Trefpi = Tpi[0]
Trefint = Tint[0]
Spi = speedup_T(Trefpi, Tpi)
Sint = speedup_T(Trefint, Tint)
fout = open("speedup.dat", "w")
print("N Sth Spi Sint", file=fout)
for (n, theor, spi, sint) in zip(x, y, Spi, Sint):
print("%d %g %g %g" % (n, theor, spi, sint), file=fout)
fout.close()
print("Times")
for (n, tpi, tpiv, tint, tintv) in zip(x, Tpi, Tpiv, Tint, Tintv):
print("%d $%.3g \pm %.3g$ $%.3g \pm %.3g$ " % (n, tpi, tpiv**0.5, tint, tintv**0.5))
plt.plot(x, y, 'o-', lw = 4.0, label='theor')
plt.plot(x, Spi, 'o-', lw = 4.0, label='pi')
plt.plot(x, Sint, 'o-', lw = 4.0, label='integral')
plt.xlabel('Number of threads $n_{tr}$')
plt.ylabel('Speed-up $S$')
plt.grid(True)
#plt.savefig('exectime.pgf', format='pgf')
#tikzplotlib.save("theor_speedup.tex")
plt.show()
c4science · Help