diff --git a/palace/fem/fespace.hpp b/palace/fem/fespace.hpp index f7fbd3033..db3e4e490 100644 --- a/palace/fem/fespace.hpp +++ b/palace/fem/fespace.hpp @@ -300,7 +300,8 @@ class AuxiliaryFiniteElementSpaceHierarchy std::vector GetDiscreteInterpolators() const { std::vector 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 index 3623f91ec..672f3e97b 100644 --- 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) {