I assumed that we are not allowed to modify the matrices received by parameter (it is very common to pass large types like a matrix by reference to const) but you could use one of them to store the result and return that instead of allocating a new one.