adding matrix vector product, working on linear systems