On finding approximate supernodes for an efficient ILU(k) factorization