diff --git a/palace/models/curlcurloperator.cpp b/palace/models/curlcurloperator.cpp
index e1385417c..ba00d3a63 100644
--- a/palace/models/curlcurloperator.cpp
+++ b/palace/models/curlcurloperator.cpp
@@ -128,7 +128,8 @@ std::unique_ptr<Operator> CurlCurlOperator::GetStiffnessMatrix()
     auto &nd_fespace_l = nd_fespaces.GetFESpaceAtLevel(l);
     if (print_hdr)
     {
-      Mpi::Print(" Level {:d}: {:d} unknowns", l, nd_fespace_l.GlobalTrueVSize());
+      Mpi::Print(" Level {:d} (p = {:d}): {:d} unknowns", l,
+                 nd_fespace_l.GetMaxElementOrder(), nd_fespace_l.GlobalTrueVSize());
     }
     constexpr auto MatType = MaterialPropertyType::INV_PERMEABILITY;
     MaterialPropertyCoefficient<MatType> muinv_func(mat_op);
diff --git a/palace/models/laplaceoperator.cpp b/palace/models/laplaceoperator.cpp
index 799c42b64..146f5cd94 100644
--- a/palace/models/laplaceoperator.cpp
+++ b/palace/models/laplaceoperator.cpp
@@ -150,7 +150,8 @@ std::unique_ptr<Operator> LaplaceOperator::GetStiffnessMatrix()
     auto &h1_fespace_l = h1_fespaces.GetFESpaceAtLevel(l);
     if (print_hdr)
     {
-      Mpi::Print(" Level {:d}: {:d} unknowns", l, h1_fespace_l.GlobalTrueVSize());
+      Mpi::Print(" Level {:d} (p = {:d}): {:d} unknowns", l,
+                 h1_fespace_l.GetMaxElementOrder(), h1_fespace_l.GlobalTrueVSize());
     }
     constexpr auto MatType = MaterialPropertyType::PERMITTIVITY_REAL;
     MaterialPropertyCoefficient<MatType> epsilon_func(mat_op);
diff --git a/palace/models/spaceoperator.cpp b/palace/models/spaceoperator.cpp
index ef3c9ab30..8022247f6 100644
--- a/palace/models/spaceoperator.cpp
+++ b/palace/models/spaceoperator.cpp
@@ -688,7 +688,8 @@ std::unique_ptr<OperType> SpaceOperator::GetPreconditionerMatrix(double a0, doub
       auto &fespace_l = fespaces.GetFESpaceAtLevel(l);
       if (print_prec_hdr)
       {
-        Mpi::Print(" Level {:d}{}: {:d} unknowns", l, (s == 0) ? "" : " (auxiliary)",
+        Mpi::Print(" Level {:d}{} (p = {:d}): {:d} unknowns", l,
+                   (s == 0) ? "" : " (auxiliary)", fespace_l.GetMaxElementOrder(),
                    fespace_l.GlobalTrueVSize());
       }
       const int sdim = GetNDSpace().GetParMesh()->SpaceDimension();