Implementation of Block Preconditioners for Piezoelasticity in SfePy

Motivated by the need to model 3D printed electroactive metamaterials, an approach is presented to improve the computational efficiency of the SfePy package for the simulation of structurally complex piezoelectric lattice-like structures.Enhancements to the iterative solver interface allow the use of state-of-the-art iterative solvers and Schur complement-based preconditioners.Their performance is illustrated by a numerical study where the convergence of the iterative solvers is observed for several levels of uniform mesh refinement.