Hi Gabe, SuiteSparse which comes together with Julia includes a library KLU which does sparse LU factorization. It has the option to return the fill-in reducing order it manages to find for a given sparsity pattern. With this returned structure it can subsequently factorize many matrices (and solve equations). After some tinkering, I have got it to work, so it is possible.
If this is still an issue and using KLU fits the bill, I'll provide more details of my tinkering. DG On Friday, July 8, 2016 at 4:40:45 AM UTC-4, Gabriel Goh wrote: > > Hey, > > I have a sequence of sparse matrix factorizations I need to do, each one a > different matrix but with the same sparsity structure. Is there a way I can > save the AMD (or any other) ordering that sparsesuite returns, it does not > need to be recomputed each time? > > Thanks a lot for the help! > > Gabe >
