A Supernodal Approach to Incomplete LU Factorization with Partial Pivoting