Revise option to fully remove FiniteElementSpace associated with high…

…-order nodes and rely on mesh vertices instead
awslabs · Aug 23, 2023 · 672ccd3 · 672ccd3
1 parent fcef953
commit 672ccd3
Showing 1 changed file with 12 additions and 3 deletions.
diff --git a/palace/utils/geodata.cpp b/palace/utils/geodata.cpp
@@ -1110,10 +1110,19 @@ std::unique_ptr<mfem::Mesh> LoadMesh(const std::string &path, bool remove_curvat
     }
     mesh = std::make_unique<mfem::Mesh>(fi, 1, 1, true);
   }
-  mesh->EnsureNodes();
-  if (remove_curvature && mesh->GetNodes()->FESpace()->GetMaxElementOrder() > 1)
+  if (remove_curvature)
   {
-    mesh->SetCurvature(1, false);
+    mfem::GridFunction *nodes = nullptr;
+    int own_nodes = true;
+    mesh->SwapNodes(nodes, own_nodes);
+    if (own_nodes)
+    {
+      delete nodes;
+    }
+  }
+  else
+  {
+    mesh->EnsureNodes();
   }
   return mesh;
 }