PARALLEL L₀-STABLE METHODS FOR THE MULTI-DIMENSIONAL DIFFUSION EQUATION

On the numerical solution of the diffusion equation subject to the specification of mass

“…Following Lawson [7] and Swayne [9], the vector yj/ The development of numerical methods will be based on making appropriate approximations to the exponentials in this recurrence relation. Higher-order Pade approximants [8] are popularly employed (see, for instance, [4,5,11]) for such exponentials. Methods based on the use of these approximants are of high accuracy in time and, in the case of the subdiagonal Pades, have good stability properties.…”