Add order print statements during multigrid level construction

awslabs · Aug 23, 2023 · e09350c · e09350c
1 parent 5d04500
commit e09350c
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 3 deletions.
diff --git 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
@@ -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
@@ -675,7 +675,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();