On the parallel solution of parabolic equations

We propose new parallel algorithms for the solution of linear parabolic problems. The first of these methods is based on using polynomial approximation to the exponential. It does not require solving any linear systems and is highly parallelizable. The two other methods proposed are based on Padé and Chebyshev approximations to the matrix exponential. The parallelization of these methods is achieved by using partial fraction decomposition techniques to solve the resulting systems and thus offers the potential for increased time parallelism in time dependent problems. We also present experimental results from the Alliant FX/8 and the Cray Y-MP/832 vector multiprocessors.

[1]  H. Tal-Ezer.  Spectral methods in time for parabolic problems , 1989 .

[2]  N. Pomphrey,et al.  DIAGONAL PAD ] APPROXIMATIONS FOR INITIAL VALUE PROBLEMS * , 1988 .

[3]  R. Varga,et al.  Extended numerical computations on the “1/9” conjecture in rational approximation theory , 1984 .

[4]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[5]  R. Varga,et al.  Chebyshev rational approximations to e−x in [0, +∞) and applications to heat-conduction problems , 1969 .

[6]  Richard S. Varga,et al.  SOME RESULTS IN APPROXIMATION THEORY WITH APPLICATIONS TO NUMERICAL ANALYSIS , 1971 .

[7]  C. Loan,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix , 1978 .

[8]  G. Fairweather.  A note on the efficient implementation of certain Padé Methods for linear parabolic problems , 1978 .

[9]  S. P. N. Rsett.  Restricted Pade approximations to the exponential function , 1978 .

[10]  H. T. Kung.  New Algorithms and Lower Bounds for the Parallel Evaluation of Certain Rational Expressions and Recurrences , 1976, JACM.

[11]  D J Kuck,et al.  Parallel Supercomputing Today and the Cedar Approach , 1986, Science.

[12]  Henry C. Thacher,et al.  Applied and Computational Complex Analysis. , 1988 .

[13]  V. Zakian.  Properties of Imn and Jmn approximants and applications to numerical inversion of Laplace transforms and initial value problems , 1975 .

[14]  B. Nour-Omid.  Applications of the Lanczos method , 1989 .

[15]  Yousef Saad,et al.  A parallel block cyclic reduction algorithm for the fast solution of elliptic equations , 1987, Parallel Comput..

[16]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[17]  R. Varga.  On Higher Order Stable Implicit Methods for Solving Parabolic Partial Differential Equations , 1961 .

[18]  M. Schaefer.  A polynomial based iterative method for linear parabolic equations , 1990 .

[19]  S. P. Nørsett.  Restricted Pad Approximations to the Exponential Function , 1978 .

[20]  R. Freund.  On conjugate gradient type methods and polynomial preconditioners for a class of complex non-hermitian matrices , 1990 .

[21]  W. E. Culham,et al.  A comparison of Crank-Nicolson and Chebyshev rational methods for numerically solving linear parabolic equations , 1972 .

[22]  David A. Swayne,et al.  High-Order Near Best Uniform Approximations to the Solution of Heat Conduction Problems , 1980, IFIP Congress.

[23]  Neil Pomphrey,et al.  Diagonal Padé Approximations for Initial Value Problems , 1988 .