diff --git a/fem/fe.cpp b/fem/fe.cpp
index 568ee18e68..d0d3825e13 100644
--- a/fem/fe.cpp
+++ b/fem/fe.cpp
@@ -1395,6 +1395,8 @@ void VectorFiniteElement::LocalRestriction_ND(
 PointFiniteElement::PointFiniteElement()
    : NodalFiniteElement(0, Geometry::POINT, 1, 0)
 {
+   lex_ordering.SetSize(1);
+   lex_ordering[0] = 0;
    Nodes.IntPoint(0).x = 0.0;
 }