Abstract

A new multigrid or incomplete LU technique is developed in this paper for solving large sparse algebraic systems from discretizing partial differential equations. By exploring some deep connection between the hierarchical basis method and incomplete LU decomposition, the resulting algorithm can be effectively applied to problems discretized on completely unstructured grids. Numerical experiments demonstrating the efficiency of the method are also reported.