doflocs for p0 elements should be in the center of gravity

kinnala · Apr 11, 2024 · a477591 · a477591
1 parent a0d679e
commit a477591
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/skfem/element/element_tet/element_tet_p0.py b/skfem/element/element_tet/element_tet_p0.py
@@ -10,7 +10,7 @@ class ElementTetP0(ElementH1):
     interior_dofs = 1
     maxdeg = 0
     dofnames = ['u']
-    doflocs = np.array([[.5, .5, .5]])
+    doflocs = np.array([[.25, .25, .25]])
     refdom = RefTet
 
     def lbasis(self, X, i):

diff --git a/skfem/element/element_tri/element_tri_p0.py b/skfem/element/element_tri/element_tri_p0.py
@@ -10,7 +10,7 @@ class ElementTriP0(ElementH1):
     interior_dofs = 1
     maxdeg = 0
     dofnames = ['u']
-    doflocs = np.array([[.5, .5]])
+    doflocs = np.array([[1/3, 1/3]])
     refdom = RefTri
 
     def lbasis(self, X, i):