UMFPACK: sparse multifrontal LU factorization