WIP: Add allocation free matrix multiplaction and solving

Just to give some examples to start the discussion. In particular the allocation free sparse multiplication is useful for the iterative solvers that have been discussed on the mailing list lately.