Matrix: Added support for the operations A = trans(B), A += trans(B), A -= trans(B)

Also works in the case that A and B are the same matrix.
Currently this operation is fairly slow, because the transposition is performed
in main memory. On the other hand, it's certainly better to have a slow
implementation rather than no implementation at all...