A Fast Solver Free of Fill-In for Finite Element Problems

A new algorithm for solving finite element problems is presented. It blends a preconditioned conjugate gradient iteration into a direct factorization method. The goal was to reduce fill to a negligible level and thus reduce storage requirements, but the algorithm turned out to be faster than its rivals for an important class of problems.