Leveraging Python tensor contraction packages for evaluating finite element weak forms

Abstract

: A new implementation of finite element matrix evaluation functions in the finite element code SfePy is introduced, leveraging several Python tensor contraction packages that implement a general function for evaluating expressions given using the Einstein summation convention. An example of the new weak form implementation is shown, and then results of a numerical study are presented comparing performance of the new implementation, the original implementation, and FEniCS suite.