A lot of big sparse matrices

[stevengj]

@Avec

The matrix in my original problem really is sparse, but its construction is too long to put it here.

In that case, you lose the sparse data structure when you multiply by Sm, unless you use Sm = sparse(rand(m,m))… And pom is also not a sparse matrix.