def S_c(c): from Simulator import * import numpy as np s = Simulator(50000, 5.0e-3, c) s.trajectory() #s.draw_trajectory() #s.draw_energy() #m = s.force_dipole() #print m #c = s.correlation(10000) #print c t,x = s.correlation_spectrum(110.0,100) spec = s.correlation_integral(100.0,t,x) return spec
cluster.close() res = np.sum(res) / float(n) print conc, res f = open('/home/dima/Dropbox/Results/Force_dipoles/S_c','a') f.write(str(conc)+' '+str(res)+'\n') f.close() c_max, n, k = 50.0, 30, 5 c = np.arange(0,c_max,c_max/(float(k))) for i in range(k): print 'running '+str(i)+'th cycle.' #S(c[i],n) #Simulator.plot_graph(5) s = Simulator(10000, 5.0e-3,0.0) s.trajectory() s.draw_trajectory() #s.draw_energy() #m = s.force_dipole() #print m #c = s.correlation(10000) #print c t,x = s.correlation_spectrum(200.0,100) spec = s.correlation_integral(5.0,t,x) print spec s.draw_correlation(t,x)