def hexPhaseSpace3(numRepeats=1): ### nuTrans sets both non parallel axes x = np.linspace(0.4, 0.6, 5, endpoint=False) y = np.linspace(0, 1, 20) results = np.zeros((len(x), len(y))) ### List of all results subRes = np.zeros( numRepeats) ### 50 repeats of same nu values (different atria) np.random.seed(1259884) seeds = np.random.randint( 0, 1000000, (len(x), len(y), numRepeats)) #### Generates a bunch of random seeds for use for i in range(len(x)): print x[i] for j in range(len(y)): for k in range(numRepeats): ### Need some way of defining seeds A = Atrium( tTot=10**5, nu_downLeft=y[j], nu_downRight=y[j], nu_parallel=x[i], hexagonal=True, seed_prop=seeds[i][j][k], seed_dysfunc=3 * seeds[i][j][k], seed_connections=7 * seeds[i][j][k] ) ### Do we modify the times and stuff for scaling??? A.cmpFull() subRes[k] = float(A.totalTimeAF) / A.tTot results[i][j] = np.average( subRes ) ### Each nu pair has a len(numRepeats) array of different with open(os.getcwd() + '/PickleData/hexPhaseSpace%s.pkl' % x[0], 'wb') as writeData: pickle.dump(results, writeData)