Add allocation free matrix multiplaction and solving