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).