SparseArrays.jl

SparseArrays.jl copied to clipboard

Hi,

it would be nice to have dmperm implemented as well.

I coudn't find a julia implementation and am using a python workaround for now.

We don't ship CSparse unfortunately, which would make this easy. A contribution would be welcome. The code can probably be translated - and being LGPL, there is no licensing concern.

Do any of the other suitesparse libraries have a dmperm in them?

No not as far as I can tell. I can translate in another package? It's probably also in his textbook, in which case I may be able to do one from scratch (I have several others I need to do).

I'll also wrap CSparse and CXSparse in another package just to provide that functionality in the meantime.

doesn't CSparse use zero based indexing? so any call to CSparse needs to be converted to a 0 based indexing?

All of our C libs are like this.

I have it wrapped, I'll release it Monday or Tuesday.

Whenever we have the wrapped packages, we should note them in the README for this package.

We already ship BTF from SuiteSparse. So the solution outlined in https://discourse.julialang.org/t/block-triangular-form-permutation/818 should be usable here.

PR anyone?

With respect to BTF I have a BlockTriangularForm.jl about to be released. It's a direct translation of the SuiteSparse version, just writing a few more tests.

For dmperm I was hitting some weird issues with CSparse, I will release the wrapper this weekend