An efficient direct solver for the boundary concentrated FEM in 2D
Boris N. Khoromskij and Jens Markus Melenk
Contact the author: Please use for correspondence this email.
Submission date: 05. Oct. 2001
published in: Computing, 69 (2002) 2, p. 91-117
DOI number (of the published article): 10.1007/s00607-002-1452-2
MSC-Numbers: 65N35, 65F10, 35D10
Keywords and phrases: hp-finite element methods, meshes refined towards boundary, direct solvers, schur complement
Download full preprint: PDF (1134 kB), PS ziped (531 kB)
The boundary concentrated FEM, a variant of the hp-version of the finite element method, is proposed for the numerical treatment of elliptic boundary value problems. It is particularly suited for equations with smooth coefficients and boundary conditions that may have low Sobolev regularity. In the two-dimensional case, it is shown that the Cholesky factorization of the resulting stiffness matrix requires O( N log4N) units of storage and can be computed with O( N log8N) work. Numerical results confirm theoretical estimates.