Matrix decomposition algorithms for arbitrary order CO tensor product finite element systems

• Autores: Kui Du, Graeme Fairweather, Weiwei Sun
• Localización: Journal of computational and applied mathematics, ISSN 0377-0427, Vol. 275, Nº 1, 2015, págs. 162-182
• Idioma: inglés
• DOI: 10.1016/j.cam.2014.08.015
• Texto completo no disponible (Saber más ...)
• Resumen

  • Matrix decomposition algorithms (MDAs) are fast direct methods for the solution of systems of linear algebraic equations which arise in the approximation of Poisson�fs equation on the unit square using various techniques such as finite difference, spline collocation and spectral methods. The attraction of MDAs is that they employ fast Fourier transforms and require O(N2 log N) operations on an N �~ N uniform partition of the unit square. In this paper, MDAs are formulated for the solution of the finite element Galerkin equations arising when spaces of C0 piecewise polynomials of degree k . 3 are employed. Results of numerical experiments exhibit the expected optimal global convergence rates and superconvergence phenomena.