Ejemplo n.º 1
0
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)
Ejemplo n.º 2
0
def constantRotation(position, velocity, rotation, omega):
    return utils.pbc(position + velocity * dt, sys.grid.length), rotation