Add A_mul_B!(output, A, b) etc for sparse A