I need to preform some sparse matrix-matrix calculations, but I only need a few of the non-zeros. Specifically, for `A*B=C`

I only need the resulting `C`

nonzeros where `A`

has nonzeros. Forming all of `C`

and then dropping the values I don’t need is prohibitively memory expensive.

It is my understanding that Julia makes use of `cholmod_ssmult`

for this computation. What I’m asking is not implemented in `cholmod_ssmult`

as far as I can tell. Are there any other libraries that can accomplish this? Is there a Julia way to accomplish this?