def reaction_function(topology): recipe = top.Recipe(topology) if topology.get_n_particles() == 1: recipe.change_particle_type(0, "B") return recipe
def reaction_function(topology): recipe = top.Recipe(topology) if topology.get_n_particles() > 1: edge = np.random.randint(0, topology.get_n_particles() - 1) recipe.remove_edge(edge, edge + 1) return recipe