pass
if DNA:
throw_in_particles(s.world, Ad, N_PARTICLES_FACTOR * 1)
throw_in_particles(s.world, Bd, N_PARTICLES_FACTOR * 1)
throw_in_particles(s.world, Cd, N_PARTICLES_FACTOR * 4)
'''Define loggers.
'''
if VTK_LOGGER == True:
from vtklogger import VTKLogger
# Write vtk files. See vtklogger.py. Use VTK or Paraview to
# visualize.
vtklogger = VTKLogger(s, 'data/run', extra_particle_step=True,
buffer_size=20)
if LOGGER == True:
from logger import Logger
# Write log files. See logger.py. Use visualizer.py to visualize.
l = Logger('example')
l.start(s)
#l.set_particle_out_interval(1e-4) #1e-3)
'''Run simulation.
'''
def run():
print s.dump_reaction_rules()
