def getRand(b_):
"""Returns uniformly distributed vectors on the surface of a sphere with radius b_."""
x1 = x2 = 1.0
while (x1**2 + x2**2 >= 1.0):
x1 = 2.0 * random.random() - 1.0
x2 = 2.0 * random.random() - 1.0
x = 2.0 * x1 * (1 - x1**2 - x2**2)**0.5
y = 2.0 * x2 * (1 - x1**2 - x2**2)**0.5
z = 1.0 - 2.0 * (x1**2 + x2**2)
return (x * b_, y * b_, z * b_)
for file in blockfiles:
fmanager = pypolymer.read_blockfile(file)
timestep = fmanager[0]
total_particles = fmanager[1]
M_rings = fmanager[2]
M_linear = fmanager[3]
monomers = fmanager[4]
monomers_half = fmanager[5]
x = fmanager[6]
y = fmanager[7]
z = fmanager[8]
if (file == blockfiles[0]):
total_particles1st = total_particles
sys.stdout.write('System parameters:\n')
sys.stdout.write(' M_rings = ' + str(M_rings) + '\n')
sys.stdout.write(' M_linear = ' + str(M_linear) + '\n')
def getRand(b_):
"""Returns uniformly distributed vectors on the surface of a sphere with radius b_."""
x1 = x2 = 1.0
while x1 ** 2 + x2 ** 2 >= 1.0:
x1 = 2.0 * random.random() - 1.0
x2 = 2.0 * random.random() - 1.0
x = 2.0 * x1 * (1 - x1 ** 2 - x2 ** 2) ** 0.5
y = 2.0 * x2 * (1 - x1 ** 2 - x2 ** 2) ** 0.5
z = 1.0 - 2.0 * (x1 ** 2 + x2 ** 2)
return (x * b_, y * b_, z * b_)
for file in blockfiles:
fmanager = pypolymer.read_blockfile(file)
timestep = fmanager[0]
total_particles = fmanager[1]
M_rings = fmanager[2]
M_linear = fmanager[3]
monomers = fmanager[4]
monomers_half = fmanager[5]
x = fmanager[6]
y = fmanager[7]
z = fmanager[8]
if file == blockfiles[0]:
total_particles1st = total_particles
sys.stdout.write("System parameters:\n")
sys.stdout.write(" M_rings = " + str(M_rings) + "\n")
sys.stdout.write(" M_linear = " + str(M_linear) + "\n")
