X(matrix_norm2) should not destroy the matrix by default

Add an optional argument preserve_mat to preserve the matrix. At the moment, it is destroyed on exit