diff --git a/src/serac/numerics/functional/element_restriction.cpp b/src/serac/numerics/functional/element_restriction.cpp index f7f1a6db3..e5c613364 100644 --- a/src/serac/numerics/functional/element_restriction.cpp +++ b/src/serac/numerics/functional/element_restriction.cpp @@ -216,7 +216,8 @@ std::vector > geom_local_face_dofs(int p) return output; } -axom::Array GetElementRestriction(const serac::fes_t* fes, mfem::Geometry::Type geom) +axom::Array GetElementRestriction(const serac::fes_t* fes, + mfem::Geometry::Type geom) { std::vector elem_dofs{}; mfem::Mesh* mesh = fes->GetMesh(); @@ -271,7 +272,7 @@ axom::Array GetElementRestriction(cons if (n == 0) { return axom::Array{}; } else { - uint64_t dofs_per_elem = elem_dofs.size() / n; + uint64_t dofs_per_elem = elem_dofs.size() / n; axom::Array output(n, dofs_per_elem); std::memcpy(output.data(), elem_dofs.data(), sizeof(DoF) * n * dofs_per_elem); return output; @@ -279,7 +280,7 @@ axom::Array GetElementRestriction(cons } axom::Array GetElementDofs(const serac::fes_t* fes, mfem::Geometry::Type geom, - const std::vector& mfem_elem_ids) + const std::vector& mfem_elem_ids) { std::vector elem_dofs{}; @@ -337,15 +338,15 @@ axom::Array GetElementDofs(const serac if (n == 0) { return axom::Array{}; } else { - uint64_t dofs_per_elem = elem_dofs.size() / n; + uint64_t dofs_per_elem = elem_dofs.size() / n; axom::Array output(n, dofs_per_elem); std::memcpy(output.data(), elem_dofs.data(), sizeof(DoF) * n * dofs_per_elem); return output; } } -axom::Array GetFaceDofs(const serac::fes_t* fes, mfem::Geometry::Type face_geom, - FaceType type) +axom::Array GetFaceDofs(const serac::fes_t* fes, + mfem::Geometry::Type face_geom, FaceType type) { std::vector face_dofs; mfem::Mesh* mesh = fes->GetMesh(); @@ -452,15 +453,16 @@ axom::Array GetFaceDofs(const serac::f if (n == 0) { return axom::Array{}; } else { - uint64_t dofs_per_face = face_dofs.size() / n; + uint64_t dofs_per_face = face_dofs.size() / n; axom::Array output(n, dofs_per_face); std::memcpy(output.data(), face_dofs.data(), sizeof(DoF) * n * dofs_per_face); return output; } } -axom::Array GetFaceDofs(const serac::fes_t* fes, mfem::Geometry::Type face_geom, - const std::vector& mfem_face_ids) +axom::Array GetFaceDofs(const serac::fes_t* fes, + mfem::Geometry::Type face_geom, + const std::vector& mfem_face_ids) { std::vector face_dofs; mfem::Mesh* mesh = fes->GetMesh(); @@ -619,7 +621,7 @@ axom::Array GetFaceDofs(const serac::f if (n == 0) { return axom::Array{}; } else { - uint64_t dofs_per_face = face_dofs.size() / n; + uint64_t dofs_per_face = face_dofs.size() / n; axom::Array output(n, dofs_per_face); std::memcpy(output.data(), face_dofs.data(), sizeof(DoF) * n * dofs_per_face); return output; diff --git a/src/serac/numerics/functional/element_restriction.hpp b/src/serac/numerics/functional/element_restriction.hpp index d6e317645..9444c08a9 100644 --- a/src/serac/numerics/functional/element_restriction.hpp +++ b/src/serac/numerics/functional/element_restriction.hpp @@ -243,7 +243,8 @@ struct BlockElementRestriction { * @param fes the finite element space containing the dof information * @param geom the kind of element geometry */ -axom::Array GetElementDofs(const serac::fes_t* fes, mfem::Geometry::Type geom); +axom::Array GetElementDofs(const serac::fes_t* fes, + mfem::Geometry::Type geom); /** * @brief Get the list of dofs for each face element (of the specified geometry) from the fes_t @@ -252,9 +253,10 @@ axom::Array GetElementDofs(const serac * @param geom the kind of element geometry * @param type whether the face is of interior or boundary type */ -axom::Array GetFaceDofs(const serac::fes_t* fes, mfem::Geometry::Type face_geom, - FaceType type); +axom::Array GetFaceDofs(const serac::fes_t* fes, + mfem::Geometry::Type face_geom, FaceType type); /// @overload -axom::Array GetFaceDofs(const serac::fes_t* fes, mfem::Geometry::Type face_geom, - const std::vector& mfem_face_ids); +axom::Array GetFaceDofs(const serac::fes_t* fes, + mfem::Geometry::Type face_geom, + const std::vector& mfem_face_ids);