def rot(point, axis, angle): q = Quaternion(axis, angle) P = point.copy() P.rotate(q) #print(point, P) return P