def solverParticlePos(position, velocity, rotation, omega): position_new = utils.pbc(position + velocity * dt, sys.grid.length) rotation_new = rotation + dt * sys.sloverRotation(omega, rotation) return position_new, sys.normalize(rotation_new)
def constantRotation(position, velocity, rotation, omega): return utils.pbc(position + velocity * dt, sys.grid.length), rotation