Faster sparse matrix multiplication (#8316)
This uses sortSparseMatrixCSC!(S) instead of double transpose to sort the result.
(Viral committing an improved version of Glen Chisholm's original PR).