Stable broken 𝐻(𝑐𝑢𝑟𝑙) polynomial extensions and 𝑝-robust a posteriori error estimates by broken patchwise equilibration for the curl–curl problem

We study extensions of piecewise polynomial data prescribed in a patch of tetrahedra sharing an edge. We show stability in the sense that the minimizers over piecewise polynomial spaces with prescribed tangential component jumps across faces and prescribed piecewise curl in elements are subordinate in the broken energy norm to the minimizers over the broken H ( curl ) \boldsymbol H(\boldsymbol {\operatorname {curl}}) space with the same prescriptions. Our proofs are constructive and yield constants independent of the polynomial degree. We then detail the application of this result to the a posteriori error analysis of the curl–curl problem discretized with Nédélec finite elements of arbitrary order. The resulting estimators are reliable, locally efficient, polynomial-degree-robust, and inexpensive. They are constructed by a broken patchwise equilibration which, in particular, does not produce a globally H ( curl ) \boldsymbol H(\boldsymbol {\operatorname {curl}}) -conforming flux. The equilibration is only related to edge patches and can be realized without solutions of patch problems by a sweep through tetrahedra around every mesh edge. The error estimates become guaranteed when the regularity pick-up constant is explicitly known. Numerical experiments illustrate the theoretical findings.

Spanning the isogeny class of a power of an elliptic curve

Let E E be an ordinary elliptic curve over a finite field and g g be a positive integer. Under some technical assumptions, we give an algorithm to span the isomorphism classes of principally polarized abelian varieties in the isogeny class of E g E^g . The varieties are first described as hermitian lattices over (not necessarily maximal) quadratic orders and then geometrically in terms of their algebraic theta null point. We also show how to algebraically compute Siegel modular forms of even weight given as polynomials in the theta constants by a careful choice of an affine lift of the theta null point. We then use these results to give an algebraic computation of Serre’s obstruction for principally polarized abelian threefolds isogenous to E 3 E^3 and of the Igusa modular form in dimension 4 4 . We illustrate our algorithms with examples of curves with many rational points over finite fields.