There is an undocumented function SparseArray`SparseLUDecomposition which, i believe, does a LU decomposition on a sparse array, returning a sparse array (?) But when I apply this function on a test sparse array, it does nothing... Maybe it requires further arguments than the sparse array?
Does anyone have further information on this function?