cross output matches input

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