results = [] P = [] for size in sizes: p = 0.01 results.append([]) print 'Running for size =', size while p < 1: stdout.write("-") stdout.flush() perc_count = 0 for k in range(trial_count): g = random_grid(size, p, 2) flow, perc = percolation_hk(g, True) if perc: perc_count += 1 #endif #endfor prob = float(perc_count) / trial_count #print 'percolation q=',prob results[n].append(prob) P.append(p) p += step #endwhile stdout.write(">|" + "\n") n += 1 #endfor
from percolation_recursive import * from percolation_hk import * from sys import * step = 0.01 trial_count = 20 size = 20 p = 0 print 'Running for size =', size data = open("percolation.data", 'w') data.write("--------size = %d--------\n" % size) data.write("Occupation \t spanning\n") data.write("probability \t probability\n") while p < 1: stdout.write("-") stdout.flush() perc_count = 0 for k in range(trial_count): g = random_grid(size, p, 2) #flow,perc = percolation_wave(g) #flow,perc = percolation_recursive(g) flow, perc = percolation_hk(g) if perc: perc_count += 1 prob = float(perc_count) / trial_count data.write("%f\t %f\n" % (p, prob)) p += step stdout.write(">|" + "\n") data.close()
from percolation_recursive import * from percolation_hk import * from sys import * step = 0.01 trial_count = 20 size = 20 p = 0 print 'Running for size =', size data = open("percolation.data", 'w') data.write("--------size = %d--------\n" %size) data.write("Occupation \t spanning\n") data.write("probability \t probability\n") while p < 1: stdout.write("-") stdout.flush() perc_count = 0 for k in range(trial_count): g = random_grid(size, p, 2) #flow,perc = percolation_wave(g) #flow,perc = percolation_recursive(g) flow,perc = percolation_hk(g) if perc: perc_count += 1 prob = float(perc_count)/trial_count data.write("%f\t %f\n" % (p, prob)) p += step stdout.write(">|" + "\n") data.close()
results=[] P = [] for size in sizes: p = 0.01 results.append([]) print 'Running for size =', size while p < 1: stdout.write("-") stdout.flush() perc_count = 0 for k in range(trial_count): g = random_grid(size, p, 2) flow,perc = percolation_hk(g,True) if perc: perc_count += 1 #endif #endfor prob = float(perc_count)/trial_count #print 'percolation q=',prob results[n].append(prob) P.append(p) p += step #endwhile stdout.write(">|" + "\n") n += 1 #endfor