Parallel Rapid Operator for Iterative Finite Element Solvers on a Shared Memory Machine

Abstract The p -version of the finite element method, which uses high order hierarchic basis functions within an element, allows for ease of adaptive computation, and high accuracy with relatively few degrees of freedom. However, the solution time and storage requirements for large three-dimensional problems can be very high when a global matrix is assembled. The use of rapid operator application in iterative methods leads to minimal storage requirements, because neither elemental nor global matrices are formed. We show in this paper that this algorithm is amenable to parallel and vector computation, and is readily implemented on shared memory machines. Elements are assigned to processors as they become available, and the computations are completely dominated by matrix-matrix multiplications, which can be tuned for maximum speed. Results showing the high efficiency and scalability of such a scheme on an IBM 3090/600 with Vector Facility are given.

[1]  Richard J. Sahulka,et al.  IBM Parallel FORTRAN , 1988, IBM Syst. J..

[2]  Yvon Maday,et al.  Optimal error analysis of spectral methods with emphasis on non-constant coefficients and deformed geometries , 1990 .

[3]  S. G. Tucker,et al.  The IBM 3090 System: An Overview , 1986, IBM Syst. J..

[4]  A. Robin Forrest,et al.  Curves and surfaces for computer-aided design , 1968 .

[5]  Werner Buchholz.  The IBM System/370 Vector Architecture , 1986, IBM Syst. J..

[6]  S. A. Coons.  SURFACES FOR COMPUTER-AIDED DESIGN OF SPACE FORMS , 1967 .

[7]  Vijay Sonnad,et al.  Multilevel solution method for the p-version of finite elements , 1989 .

[8]  Thomas J. R. Hughes,et al.  Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies , 1985 .

[9]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[10]  Ivo Babuška,et al.  The p-Version of the Finite Element Method for Parabolic Equations. Part 1 , 1981 .

[11]  Vijay Sonnad,et al.  Minimal Storage Finite Element Solution of Large-Scale Three-Dimensional Elastodynamic Problems , 1992 .

[12]  Vijay Sonnad,et al.  A comparison of preconditioned iterative techniques using rapid operator application against direct solution methods , 1991 .

[13]  S. Orszag.  Spectral methods for problems in complex geometries , 1980 .

[14]  Ivo Babuška,et al.  The p - and h-p version of the finite element method, an overview , 1990 .

[15]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[16]  I. Babuska,et al.  Finite Element Analysis , 2021 .

[17]  T. Hughes,et al.  Element-by-Element Implicit Algorithms for Heat Conduction , 1983 .