Details of Research Outputs

Explicit expressions for the eigensystems of one-dimensional finite element Galerkin (FEG) matrices based on C0 piecewise quadratic polynomials are determined. These eigensystems are then used in the formulation of fast direct methods, matrix decomposition algorithms (MDAs), for the solution of the FEG equations arising from the discretization of Poisson's equation on the unit square subject to several standard boundary conditions. The MDAs employ fast Fourier transforms and require O(N2log N) operations on an N×N uniform partition. Numerical results are presented to demonstrate the efficacy of these algorithms. © Springer Science + Business Media B.V. 2009.