Sparse matrices: No more temporaries for y += Ax; and y -= Ax;

Now using two sets of SpMV kernels:
 - One for y = Ax for performance reasons
 - One for y = alpha * Ax + beta * y; for inplace-operations.
Previous use of temporaries for inplace-operations was too slow and
resulted in unnecessary memory overhead.