Functionalities not supported for matrices

PR #330 adds the first batch of implementations for operator-matrix products, but not everything is supported. It only contains barebones code for normal, transposed, adjoint and conjugated operators.
Opening this issue to keep track of the rest, and add more tests.