code: def quaternion_transform(r):
      
      Wt_r = makeW(*r).T
      Q_r = makeQ(*r)
      rot = Wt_r.dot(Q_r)[:3, :3]
      return rot
docs: Get optimal rotation
      note: translation will be zero when the centroids of each molecule are the
      same