예제 #1
0
        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()