A domain-specific language for the hybridization and static condensation of finite element methods

File Description SizeFormat 
1802.00303v1.pdfWorking paper1.29 MBAdobe PDFDownload
Title: A domain-specific language for the hybridization and static condensation of finite element methods
Author(s): Gibson, TH
Mitchell, L
Ham, DA
Cotter, CJ
Item Type: Working Paper
Abstract: In this paper, we introduce a domain-specific language (DSL) for concisely expressing localized linear algebra on finite element tensors and its integration within a code-generation framework. This DSL is general enough to facilitate the automatic generation of element-based dense linear algebra kernels necessary for the implementation of static condensation methods and local solvers for a variety of problems. We demonstrate its use for the static condensation of continuous Galerkin problems, and systems arising from hybridizing finite element discretizations. We also describe how this DSL can be used to execute local post-processing procedures to construct superconvergent approximations to mixed problems. This work also features high-level implementations of static condensation and hybrid-mixed methods as preconditioners conforming to PETSc's interface for solving linear systems. These preconditioning interfaces provide reduced operators, which are obtained from locally assembled expressions, with the necessary context to specify full solver configurations on the resulting linear systems. We validate our implementations within the context of model second-order elliptic problems. Finally, we conclude with a performance comparison highlighting the use of a hybrid-mixed method as an efficient preconditioner in a semi-implicit method for the nonlinear rotating shallow water equations.
URI: http://hdl.handle.net/10044/1/57366
Copyright Statement: © 2018 The Authors
Sponsor/Funder: Natural Environment Research Council (NERC)
Engineering & Physical Science Research Council (EPSRC)
Funder's Grant Number: NE/K006789/1
Keywords: cs.MS
G.4; G.1.8; I.1
Appears in Collections:Faculty of Engineering
Applied Mathematics and Mathematical Physics
Faculty of Natural Sciences

Items in Spiral are protected by copyright, with all rights reserved, unless otherwise indicated.

Creative Commons