Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
A comparison of two multilevel Schur preconditioners for adaptive FEM
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing.
2014 (English)Independent thesis Advanced level (professional degree), 20 credits / 30 HE creditsStudent thesis
Abstract [en]

There are several algorithms for solving the linear system of equations that arise from the finite element method with linear or near-linear computational complexity. One way is to find an approximation of the stiffness matrix that is such that it can be used in a preconditioned conjugate residual method, that is, a preconditioner to the stiffness matrix. We have studied two preconditioners for the conjugate residual method, both based on writing the stiffness matrix in block form, factorising it and then approximating the Schur complement block to get a preconditioner. We have studied the stationary reaction-diffusion-advection equation in two dimensions. The mesh is refined adaptively, giving a hierarchy of meshes. In the first method the Schur complement is approximated by the stiffness matrix at one coarser level of the mesh, in the second method it is approximated as the assembly of local Schur complements corresponding to macro triangles.

For two levels the theoretical bound of the condition number is 1/(1-C²) for either method, where C is the Cauchy-Bunyakovsky-Schwarz constant. For multiple levels there is less theory. For the first method it is known that the condition number of the preconditioned stiffness matrix is O(l²), where l is the number of levels of the preconditioner, or, equivalently, the number mesh refinements. For the second method the asymptotic behaviour is not known theoretically. In neither case is the dependency of the condition number of C known. We have tested both methods on several problems and found the first method to always give a better condition number, except for very few levels. For all tested problems, using the first method it seems that the condition number is O(l), in fact it is typically not larger than Cl. For the second method the growth seems to be superlinear.

Place, publisher, year, edition, pages
2014. , 38 p.
Series
UPTEC F, ISSN 1401-5757 ; 14006
Keyword [en]
computational science, PDE, partial differential equations, FEM, finite element method, adaptive mesh refinement, Schur complement, multilevel method, iterative method, conjugate gradient method
Keyword [sv]
beräkningsvetenskap, PDE, partiella differentialekvationer, FEM, finita elementmetoden, adaptivt beräkningsnät, adaptiv mesh, Schurkomplement, multilevelmetod, iterativ metod, konjugerade gradientmetoden
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:uu:diva-219939OAI: oai:DiVA.org:uu-219939DiVA: diva2:703651
Educational program
Master Programme in Engineering Physics
Supervisors
Examiners
Available from: 2014-03-18 Created: 2014-03-07 Last updated: 2014-03-18Bibliographically approved

Open Access in DiVA

fulltext(737 kB)290 downloads
File information
File name FULLTEXT01.pdfFile size 737 kBChecksum SHA-512
b8eb5d55327a20d22972e52ade5701559c1c9cd3bb958abfa4f1f59afbf8ca8c09842cc19c429d61af8aef4b423bcb6d2ee62d1e4e819d4f0be2fec95e47611c
Type fulltextMimetype application/pdf

By organisation
Division of Scientific Computing
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 290 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

urn-nbn

Altmetric score

urn-nbn
Total: 594 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf