#700 Faster (?) sparsity pattern calculation for matrix products
#700 Quicker sparsity pattern calculation if multipying diagonal matrices