Fast Sparse Matrix Permutation for Mesh-Based Direct Solvers
This work accelerates direct solvers for mesh-based problems in graphics, offering a practical speedup for practitioners.
The authors present a fast sparse matrix permutation algorithm for triangle mesh linear systems that reduces permutation time and improves sparse Cholesky solve performance by up to 6.27x across graphics applications.
We present a fast sparse matrix permutation algorithm tailored to linear systems arising from triangle meshes. Our approach produces nested-dissection-style permutations while significantly reducing permutation runtime overhead. Rather than enforcing strict balance and separator optimality, the algorithm deliberately relaxes these design decisions to favor fast partitioning and efficient elimination-tree construction. Our method decomposes permutation into patch-level local orderings and a compact quotient-graph ordering of separators, preserving the essential structure required by sparse Cholesky factorization while avoiding its most expensive components. We integrate our algorithm into vendor-maintained sparse Cholesky solvers on both CPUs and GPUs. Across a range of graphics applications, including single factorizations and repeated factorizations, our method reduces permutation time and improves the sparse Cholesky solve performance by up to 6.27x. Our code is available at https://github.com/BehroozZare/fast-permute.