It's not clear to me that this is the equivalent of the 2D function of the same name -- I think that actually finds a unit vector that is normal to the line between the two input vectors, as opposed to being normal to both of them.