  One conforming and one non-conforming virtual element Hessian complexes on
tetrahedral grids are constructed based on decompositions of polynomial tensor
space. They are applied to discretize the linearized time-independent
Einstein-Bianchi system.
