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')
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')