Comparison of the solvers of systems of linear equations: FETI-DP vs BDDC

Abstract

This bachelor thesis compares the dual–primal domain decomposition methods FETI-DP and BDDC for solving systems of linear equations arising from finite element discretizations of elliptic problems. The work is based on a model Poisson equation with mixed Dirichlet–Neumann boundary conditions on a two-dimensional unit square domain. After deriving the weak formula tion and performing the finite element discretization, the discrete problem is reduced by subdomain decomposition and Schur complements to a problem in the space corresponding to the interface between subdomains, which makes it possible to describe both methods in a common framework. Within this framework, sequential implementations of both methods are developed in Octave and compared according to the spectral properties of the preconditioned operators and the convergence of the PCG method. The study is complemented by a benchmark comparison of the official parallel implementations available in the PETSc library. The results show that both methods are numerically robust and behave very similarly on the tested problems.

Táto bakalárska práca porovnáva dual-primal metódy rozkladu oblasti FETI-DP a BDDC určené na riešenie sústav lineárnych rovníc vznikajúcich pri diskretizácii eliptických úloh metódou konečných prvkov. Práca vychádza z modelovej Poissonovej rovnice so zmiešanými Dirichletovými a Neumannovými okrajovými podmienkami na dvojrozmernej štvorcovej oblasti s jednotkovými rozmermi. Po odvodení slabej formulácie a diskretizácii metódou konečných prvkov je diskrétny problém pomocou rozkladu na subdomény a Schurových doplnkov redukovaný na problém v pries tore zodpovedajúcom rozhraniu podoblastí, čo umožňuje opísať obe metódy v spoločnom rámci. V tomto rámci sú vytvorené sekvenčné implementácie oboch metód v prostredí Octave a porovnané podľa spektrálnych vlastností predpodmienených operátorov a konvergencie metódy PCG. Štúdia je doplnená benchmarkovým porovnaním oficiálnych paralelných implementácií v knižnici PETSc. Výsledky ukazujú, že obe metódy sú na testovaných úlohách numericky robustné a veľmi podobné.