Massively Parallel Quadratic Programming Solvers with Applications in Mechanics

Abstract

This thesis focuses on practical solution of large-scale contact problems of structure mechanics by means of a derived quadratic programming (QP) formulation. An approach proposed by professor Dostál, combining a FETI-type non-overlapping domain decomposition method, the SMALBE algorithm based on augmented Lagrangians, and the MPRGP algorithm belonging to active set methods, has been adopted. This approach enjoys theoretically supported numerical scalability and a favourable potential for parallel scalability. The thesis consists of two parts: Background and Implementation.

Background is devoted to rather theoretical aspects of QP and FETI, although tightly connected to practical implementation. Original topics include QP transforms, implicit orthonormalization of equality constraints, and a minor modification of SMALBE shortening its termination phase considerably.

Second part, Implementation, deals with the massively parallel implementation of the aforementioned approach within PERMON, a new set of software libraries established by the author. The most important part, PERMON Solver Core, is formed mainly by the general-purpose QP solver PermonQP, and its extension PermonFLLOP providing support for domain decomposition. These libraries make use of and extend PETSc, an open source software framework for numerical computing. Performance of PERMON is demonstrated on several numerical experiments.

Tato dizertační práce se zaměřuje na praktické řešení rozsáhlých kontaktních úloh strukturální mechaniky přes odvozenou úlohu kvadratického programování (QP). Přebírá přístup navržený prof. Dostálem, kombinující nepřekrývající metodu rozložení oblasti typu FETI, algoritmus SMALBE založený na rozšířených Lagrangiánech a algoritmus MPRGP náležící k metodám aktivních množin. Tento přístup se těší teoreticky podložené numerické škálovatelnosti a příznivému potenciálu pro paralelní škálovatelnost. Tato práce se skládá ze dvou částí: Teoretické pozadí a Implementace.

Teoretické pozadí se věnuje spíše teoretickým stránkám QP a FETI, i když těsně spjatými s praktickou implementací. K originálním tématům patří QP transformace, implicitní ortonormalizace rovnostních omezení a drobná modifikace SMALBE podstatně zkracující jeho ukončovací fázi.

Druhá část, Implementace, je věnována masivně paralelní implementaci výše uvedeného přístupu v rámci PERMONu, nové sady softwarových knihoven založené autorem. Nejdůležitější částí PERMONu, tzv. Řešičové jádro (PERMON Solver Core), je tvořeno zejména QP řešičem pro obecné použití PermonQP a jeho extenzí PermonFLLOP, poskytující podporu pro rozložení oblasti. Tyto knihovny využívají a rozšiřují PETSc, softwarový rámec pro numerické výpočty s otevřeným kódem. Výkonnost PERMONu je demonstrována na několika numerických experimentech.