Abstract
Contemporary applications in computational science and engineering often require the solution of linear systems which may be of different sizes, shapes, and structures. The goal of this paper is to explain how two libraries, PETSc and HPDDM, have been interfaced in order to offer end-users robust overlapping Schwarz preconditioners and advanced Krylov methods featuring recycling and the ability to deal with multiple right-hand sides. The flexibility of the implementation is showcased and explained with minimalist, easy-to-run, and reproducible examples, to ease the integration of these algorithms into more advanced frameworks. The examples provided cover applications from eigenanalysis, elasticity, combustion, and electromagnetism.
| Original language | English (US) |
|---|---|
| Pages (from-to) | 277-295 |
| Number of pages | 19 |
| Journal | Computers and Mathematics with Applications |
| Volume | 84 |
| DOIs | |
| State | Published - Jan 22 2021 |
Bibliographical note
KAUST Repository Item: Exported on 2021-02-21Acknowledgements: The authors would like to thank S. Balay, J. Brown, V. Hapla, M. Knepley, and B. Smith for reviewing the successive merge requests in PETSc repository and for their feedback on this manuscript. This work was granted access to the GENCI-sponsored HPC resources of:, • TGCC@CEA under allocation A0070607519;, • IDRIS@CNRS under allocation AP010611780. Jose E. Roman was supported by the Spanish Agencia Estatal de Investigación (AEI) under project SLEPc-DA (PID2019-107379RB-I00).