Sparse LU factorization and update
