Fast evaluation of finite element weak forms using python tensor contraction packages

dc.contributor.authorCimrman, Robert
dc.date.accessioned2021-09-06T10:00:28Z
dc.date.available2021-09-06T10:00:28Z
dc.date.issued2021
dc.description.abstract-translatedIn finite element calculations, the integral forms are usually evaluated using nested loops over elements, and over quadrature points. Many such forms (e.g. linear or multi-linear) can be expressed in a compact way, without the explicit loops, using a single tensor contraction expression by employing the Einstein summation convention. To automate this process and leverage existing high performance codes, we first introduce a notation allowing trivial differentiation of multi-linear finite element forms. Based on that we propose and describe a new transpiler from Einstein summation based expressions, augmented to allow defining multi-linear finite element weak forms, to regular tensor contraction expressions. The resulting expressions are compatible with a number of Python scientific computing packages, that implement, optimize and in some cases parallelize the general tensor contractions. We assess the performance of those packages, as well as the influence of operand memory layouts and tensor contraction paths optimizations on the elapsed time and memory requirements of the finite element form evaluations. We also compare the efficiency of the transpiled weak form implementations to the C-based functions available in the finite element package SfePy.en
dc.format26 s.cs
dc.format.mimetypeapplication/pdf
dc.identifier.citationCIMRMAN, R. Fast evaluation of finite element weak forms using python tensor contraction packages. Advances in engineering software, 2021, roč. 159, č. 20. July 2021. ISSN 0965-9978.cs
dc.identifier.document-number676736200005
dc.identifier.doi10.1016/j.advengsoft.2021.103033
dc.identifier.issn0965-9978
dc.identifier.obd43933198
dc.identifier.uri2-s2.0-85110434316
dc.identifier.urihttp://hdl.handle.net/11025/45055
dc.language.isoenen
dc.project.IDEF17_048/0007280/Aplikace moderních technologií v medicíně a průmyslucs
dc.publisherElsevieren
dc.relation.ispartofseriesAdvances In Engineering Softwareen
dc.rights© Elsevieren
dc.rights.accessopenAccessen
dc.subject.translatedFinite element methoden
dc.subject.translatedTensor contractionsen
dc.subject.translatedWeak form transpileren
dc.subject.translatedNumerical performance studyen
dc.titleFast evaluation of finite element weak forms using python tensor contraction packagesen
dc.typečlánekcs
dc.typearticleen
dc.type.statusPeer-revieweden
dc.type.versionpublishedVersionen

Files

Collections