On sparse lu factorization procedures for the solution of parabolic differential equations in three space dimensions

The numerical solution of partial differential equations in 3 dimensions by finite difference methods leads to the problem of solving large order sparse structured linear systems. In this paper, a factorization procedure in algorithmic form is derived yielding direct and iterative methods of solution of some interesting boundary value problems in physics and engineering.