Add advection matrix_sparse() and a test for it.