Sparse matrix–vector multiplication
computation routine