Avoid constructing discrete gradient matrix on coarsest mesh if not n…

…eeded
awslabs · Feb 2, 2024 · 63ad43c · 63ad43c
1 parent 1ba58e4
commit 63ad43c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/palace/fem/fespace.hpp b/palace/fem/fespace.hpp
@@ -300,7 +300,8 @@ class AuxiliaryFiniteElementSpaceHierarchy
   std::vector<const Operator *> GetDiscreteInterpolators() const
   {
     std::vector<const Operator *> G_(GetNumLevels());
-    for (std::size_t l = 0; l < G_.size(); l++)
+    G_[0] = nullptr;  // No discrete interpolator for coarsest level
+    for (std::size_t l = 1; l < G_.size(); l++)
     {
       G_[l] = &GetDiscreteInterpolatorAtLevel(l);
     }

diff --git a/palace/fem/libceed/basis.cpp b/palace/fem/libceed/basis.cpp
@@ -192,7 +192,7 @@ void InitInterpolatorBasis(const mfem::FiniteElement &trial_fe,
   if constexpr (false)
   {
     std::cout << "New interpolator basis (" << ceed << ", " << &trial_fe << ", " << &test_fe
-              << ")\n";
+              << ", " << (trial_fe.GetMapType() == test_fe.GetMapType()) << ")\n";
   }
   if constexpr (false)
   {