Subcontract B635525 (Final Report)
- Univ. of California, San Diego, CA (United States)
During my visits to LLNL July 8–12 and July 22–26, 2019, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $$\textit{hp}$$ adaptive finite element calculations, and is implemented in the $$\textit{PLTMG}$$ software package, version 13. This preconditioner typically requires 3–20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $$\textit{p}$$ as $$\textit{p} → ∞$$, but remains strong for $$\textit{p}$$ ≤ 9, which is the highest polynomial degree allowed in $$\textit{PLTMG}$$, due to limitations of the numerical quadrature rules implemented in the soft- ware package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- DOE Contract Number:
- AC52-07NA27344
- OSTI ID:
- 1581900
- Report Number(s):
- LLNL-SR-800724; 1003444
- Country of Publication:
- United States
- Language:
- English
Similar Records
Final Report - Subcontract B629114
Final Report - Subcontract B6183769