Another way to add orientation information

[luoh15]

[ihnorton]

Note from talking to Huan: the return signature of _rectangular_distance_matrix has changed, as the function was edited in-place. We need to check all uses, and either move his new functionality to a separate function, or update all of the places where this function is called (if compatible with the changes).