All inputs are 2D arrays
dot-product along last axis
normalise a (3,) vector or a (n,3) array of vectors