cross output matches input

This is just an update of PR 715 with an added docstring and a requirement that the dimensions in crossed matrices be the same.