Pablo D. Brubeck

Pablo D. Brubeck, Yizhou Liang, Charles Parker

We show that the cohomology of the finite element Stokes complex consisting of piecewise polynomials spaces on an Alfeld split mesh from Fu, Guzmán, & Neilan (2020, Math. Comp., 89, 1059--1091) is isomorphic to the cohomologies of the continuous Stokes and de Rham complexes. We also construct novel "minimal" conforming finite element complexes where the $H^1$-conforming space is the lowest-order space from Guzmán & Neilan (2018, SIAM J. Numer. Anal., 56, 2826--2844) and the $L^2$-conforming space is piecewise constants. These minimal complexes also have cohomologies isomorphic to the continuous Stokes and de Rham complexes. We further construct local, bounded, cochain projections for the minimal complexes. All the results hold for strongly Lipschitz domains with nontrivial topologies and in the presence of mixed boundary conditions.

Pablo D. Brubeck, Patrick E. Farrell, Robert C. Kirby et al.

We present new finite elements for solving the Riesz maps of the de Rham complex on triangular and tetrahedral meshes at high order. The finite elements discretize the same spaces as usual, but with different basis functions, so that the resulting matrices have desirable properties. These properties mean that we can solve the Riesz maps to a given accuracy in a $p$-robust number of iterations with $\mathcal{O}(p^6)$ flops in three dimensions, rather than the naïve $\mathcal{O}(p^9)$ flops. The degrees of freedom build upon an idea of Demkowicz et al., and consist of integral moments on an equilateral reference simplex with respect to a numerically computed polynomial basis that is orthogonal in two different inner products. As a result, the interior-interface and interior-interior couplings are provably weak, and we devise a preconditioning strategy by neglecting them. The combination of this approach with a space decomposition method on vertex and edge star patches allows us to efficiently solve the canonical Riesz maps at high order. We apply this to solving the Hodge Laplacians of the de Rham complex with novel augmented Lagrangian preconditioners.