def create_all_H_and_eqcs(n): codomain = determine_codomain(n) eqc_array = [] for H in codomain: eqc_for_H = determine_equivalence_class_iterative_precompute(H) eqc_array.append(eqc_for_H) zkl.save(codomain,"H.zkl") zkl.save(eqc_array,"eqcs_for_H.zkl")
def make_unique_urg_and_pics(num_nodes):
uurgs = determine_unique_unreachable_graphs('rgs_%s.zkl' % num_nodes,num_nodes)
zkl.save(uurgs,'unique_urgs_%s.zkl' % num_nodes)
print "Finished importing. Now drawing:"
d = zkl.load('unique_urgs_%s.zkl' % num_nodes)
i = 1
for item in d:
draw_codomain_graph(item,'pic_unique_urgs_%s_graph_%d.png' % (num_nodes,i))
print '\n'
print item
i = i + 1
def create_H_and_eqcs(numH): H_array = [] eqc_array = [] for i in range(0,numH): num_nodes = np.random.randint(low = 3, high = 7) H = generate_H(num_nodes) H_array.append(H) eqc_for_H = determine_equivalence_class_iterative_precompute(H) eqc_array.append(eqc_for_H) zkl.save(H_array,"random_H.zkl") zkl.save(eqc_array,"eqcs_for_H")
},
'4': {
'1': set([(0, 1)])
},
}
#d2l.matrix_unfold(l[0],2,1,R=5, w_gap=1, h_gap=2, mname='TT1')
n = 20
dens = 0.07
for i in range(10):
print i
g = bfu.ringmore(n, bfu.dens2edgenum(dens, n))
gl = bfu.all_undersamples(g)
zkl.save(gl, 'list1.zkl')
# output file
foo = open('figures/shipfig_figure.tex', 'wb')
sys.stdout = foo
listplot('list1.zkl', width=17, R=6)
sys.stdout = sys.__stdout__ # remember to reset sys.stdout!
foo.flush()
foo.close()
PPP = os.getcwd()
os.chdir('figures')
os.system('pdflatex --shell-escape shipfig.tex 2>&1 > /dev/null')
os.system('mv shipfig.pdf /tmp/all_graphs_' + str(n) + 'nodes_density_' +
str(dens) + '_' + str(i) + '.pdf')
os.chdir(PPP)
def savegraphs(l, fname):
zkl.save(l, fname)
for dens in densities[nodes]:
print "{:2}: {:8} : {:10} {:10}".format('id', 'density', 'OCE',
'time')
if PARALLEL:
errors = pool.map(functools.partial(wrp, n=nodes, dens=dens),
range(REPEATS))
print 'done'
else:
errors = []
for i in range(REPEATS):
errors.append(wrp(i, n=nodes, dens=dens))
print 'computed'
z[dens] = errors
zkl.save(
z[dens],
socket.gethostname().split('.')[0] + '_nodes_' + str(nodes) +
'_samples_' + str(SAMPLESIZE) + '_density_' + str(dens) +
'_noise_' + NOISE_STD + '_OCE_b_' + EST + '_' + DIST + POSTFIX +
'.zkl')
print ''
print '----'
print ''
pool.close()
pool.join()
zkl.save(
z,
socket.gethostname().split('.')[0] + '_nodes_' + str(nodes) +
'_samples_' + str(SAMPLESIZE) + '_noise_' + NOISE_STD + '_OCE_b_' +
EST + '_' + DIST + POSTFIX + '.zkl')
def savegraphs(l, fname):
zkl.save(l, fname)
densities = {6: [0.2, 0.25, 0.3, 0.35],
8: [0.3],
10:[0.1],# 0.15, 0.2, 0.25, 0.3],
15:[0.25, 0.3],
20:[0.1],# 0.15, 0.2, 0.25, 0.3],
25:[0.1],
30:[0.1],
35:[0.1],
40:[0.1],
50:[0.05, 0.1],
60:[0.05, 0.1]}
for nodes in [15]:
z = {}
pool=Pool(processes=PNUM)
for dens in densities[nodes]:
print "{:2}: {:8} : {:10} {:10}".format('id', 'density', 'eq class', 'time')
e = bfutils.dens2edgenum(dens, n=nodes)
eqclasses = pool.map(functools.partial(fan_wrapper, n=nodes, k=e),
range(REPEATS))
z[dens] = eqclasses
zkl.save(z[dens],
socket.gethostname().split('.')[0]+\
'_nodes_'+str(nodes)+'_density_'+str(dens)+'_newp_.zkl')
print ''
print '----'
print ''
pool.close()
pool.join()
zkl.save(z,socket.gethostname().split('.')[0]+'_nodes_'+str(nodes)+'_newp_.zkl')
25:[0.1],
30:[0.1],
35:[0.1],
40:[0.1],
50:[0.05, 0.1],
60:[0.05, 0.1]}
for nodes in [5]:
z = {}
#pool=Pool(processes=PNUM)
for dens in densities[nodes]:
print "{:2}: {:8} : {:10} : {:10} {:10}".format('id', 'densityi(G)', 'density(H)', 'eq class', 'time')
e = bfutils.dens2edgenum(dens, n=nodes)
eqclasses = []
for x in pool.imap(functools.partial(ra_wrapper, n=nodes, k=e),
range(REPEATS)):
eqclasses.append(x)
z[dens] = eqclasses
zkl.save(eqclasses,
socket.gethostname().split('.')[0]+\
'_nodes_'+str(nodes)+'_density_'+\
str(dens)+'_'+KEY+'_.zkl')
print ''
print '----'
print ''
#pool.close()
#pool.join()
#zkl.save(z,socket.gethostname().split('.')[0]+'_nodes_'+str(nodes)+'_'+KEY+'_.zkl')
30:[0.1],
35:[0.1]}
wrp = wrapper
for nodes in [6]:
z = {}
pool=Pool(processes=PNUM)
for dens in densities[nodes]:
print "{:2}: {:8} : {:10} {:10}".format('id', 'density', 'OCE', 'time')
if PARALLEL:
errors = pool.map(functools.partial(wrp, n=nodes,
dens=dens),
range(REPEATS))
print 'done'
else:
errors = []
for i in range(REPEATS):
errors.append(wrp(i,n=nodes,dens=dens))
print 'computed'
z[dens] = errors
zkl.save(z[dens],
socket.gethostname().split('.')[0]+'_nodes_'+str(nodes)+'_samples_'+str(SAMPLESIZE)+'_density_'+str(dens)+'_noise_'+NOISE_STD+'_OCE_b_'+EST+'_'+DIST+POSTFIX+'.zkl')
print ''
print '----'
print ''
pool.close()
pool.join()
zkl.save(z,socket.gethostname().split('.')[0]+'_nodes_'+str(nodes)+'_samples_'+str(SAMPLESIZE)+'_noise_'+NOISE_STD+'_OCE_b_'+EST+'_'+DIST+POSTFIX+'.zkl')
25: [0.1],
30: [0.1],
35: [0.1],
40: [0.1],
50: [0.05, 0.1],
60: [0.05, 0.1]
}
for nodes in [15]:
z = {}
pool = Pool(processes=PNUM)
for dens in densities[nodes]:
print "{:2}: {:8} : {:10} {:10}".format('id', 'density', 'eq class',
'time')
e = bfutils.dens2edgenum(dens, n=nodes)
eqclasses = pool.map(functools.partial(fan_wrapper, n=nodes, k=e),
range(REPEATS))
z[dens] = eqclasses
zkl.save(z[dens],
socket.gethostname().split('.')[0]+\
'_nodes_'+str(nodes)+'_density_'+str(dens)+'_newp_.zkl')
print ''
print '----'
print ''
pool.close()
pool.join()
zkl.save(
z,
socket.gethostname().split('.')[0] + '_nodes_' + str(nodes) +
'_newp_.zkl')
