phynexis-deps/mfem
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Make style

Browse Source
This commit is contained in:
Justin Dong
2026-04-02 05:26:08 -07:00
parent e8f39487b9
commit a3a368dfd9
2 changed files with 21 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/fespace.hpp
+13 -13
View File
@@ -1502,19 +1502,19 @@ public:
return dynamic_cast<const L2_FECollection*>(fec) != NULL;
}
/// @brief Return true if the mesh contains only one topology, the elements are
/// all triangles or tetrahedrons, and the elements are ragged tensor elements
/// i.e. Bernstein/positive basis.
bool UsesRaggedTensorBasis() const
{
bool mixed = this->GetMesh()->IsMixedMesh();
bool simplex = (this->GetTypicalFE()->GetGeomType() == Geometry::TRIANGLE) ||
(this->GetTypicalFE()->GetGeomType() == Geometry::TETRAHEDRON);
bool positive =
dynamic_cast<const mfem::H1Pos_TriangleElement *>(this->GetTypicalFE()) ||
dynamic_cast<const mfem::H1Pos_TetrahedronElement *>(this->GetTypicalFE());
return !mixed && simplex && positive;
}
/// @brief Return true if the mesh contains only one topology, the elements are
/// all triangles or tetrahedrons, and the elements are ragged tensor elements
/// i.e. Bernstein/positive basis.
bool UsesRaggedTensorBasis() const
{
bool mixed = this->GetMesh()->IsMixedMesh();
bool simplex = (this->GetTypicalFE()->GetGeomType() == Geometry::TRIANGLE) ||
(this->GetTypicalFE()->GetGeomType() == Geometry::TETRAHEDRON);
bool positive =
dynamic_cast<const mfem::H1Pos_TriangleElement *>(this->GetTypicalFE()) ||
dynamic_cast<const mfem::H1Pos_TetrahedronElement *>(this->GetTypicalFE());
return !mixed && simplex && positive;
}
/** In variable-order spaces on nonconforming (NC) meshes, this function
controls whether strict conformity is enforced in cases where coarse
fem/integ/bilininteg_mass_kernels.hpp
+8 -4
View File
@@ -1793,8 +1793,10 @@ inline void PAMassApplyTriangle(const int NE,
const int d1d = 0,
const int q1d = 0)
{
MFEM_VERIFY(T_D1D ? T_D1D : d1d <= DeviceDofQuadLimits::Get().MAX_D1D_SIMPLEX, "");
MFEM_VERIFY(T_Q1D ? T_Q1D : q1d <= DeviceDofQuadLimits::Get().MAX_Q1D_SIMPLEX, "");
MFEM_VERIFY(T_D1D ? T_D1D : d1d <= DeviceDofQuadLimits::Get().MAX_D1D_SIMPLEX,
"");
MFEM_VERIFY(T_Q1D ? T_Q1D : q1d <= DeviceDofQuadLimits::Get().MAX_Q1D_SIMPLEX,
"");
const int BASIS_DIM = d1d * (d1d + 1) / 2;
const auto lex_map = lex_map_.Read();
@@ -1933,8 +1935,10 @@ inline void PAMassApplyTetrahedron(const int NE,
const int d1d = 0,
const int q1d = 0)
{
MFEM_VERIFY(T_D1D ? T_D1D : d1d <= DeviceDofQuadLimits::Get().MAX_D1D_SIMPLEX, "");
MFEM_VERIFY(T_Q1D ? T_Q1D : q1d <= DeviceDofQuadLimits::Get().MAX_Q1D_SIMPLEX, "");
MFEM_VERIFY(T_D1D ? T_D1D : d1d <= DeviceDofQuadLimits::Get().MAX_D1D_SIMPLEX,
"");
MFEM_VERIFY(T_Q1D ? T_Q1D : q1d <= DeviceDofQuadLimits::Get().MAX_Q1D_SIMPLEX,
"");
const int BASIS_DIM = d1d * (d1d + 1) * (d1d + 2) / 6;
const int BASIS_DIM2D = d1d * (d1d + 1) / 2;