Update polynomial bases

gridap · Jan 15, 2025 · 4e56006 · 4e56006
1 parent 4f5b12f
commit 4e56006
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/src/ReferenceFEs/AWRefFEs.jl b/src/ReferenceFEs/AWRefFEs.jl
@@ -22,7 +22,7 @@ function ArnoldWintherRefFE(::Type{T},p::Polytope,order::Integer) where T
 
   VT = SymTensorValue{2,T}
   prebasis = MonomialBasis(Val(2),VT,3,Polynomials._p_filter)
-  fb = MonomialBasis(Val(D-1),T,0,Polynomials._p_filter)
+  fb = MonomialBasis(Val(1),T,0,Polynomials._p_filter)
   cb = map(constant_field,component_basis(VT))
 
   function cmom(φ,μ,ds) # Cell and Node moment function: σ_K(φ,μ) = ∫(φ:μ)dK

diff --git a/src/ReferenceFEs/HHJRefFEs.jl b/src/ReferenceFEs/HHJRefFEs.jl
@@ -20,9 +20,9 @@ function HellanHerrmannJhonsonRefFE(::Type{T},p::Polytope,order::Integer) where
   @assert p == TRI "HellanHerrmannJhonson Reference FE only defined for TRIangles"
 
   VT = SymTensorValue{2,T}
-  prebasis = MonomialBasis{2}(VT,order,Polynomials._p_filter)
-  fb = MonomialBasis{1}(T,order,Polynomials._p_filter)
-  cb = MonomialBasis{2}(VT,order-1,Polynomials._p_filter)
+  prebasis = MonomialBasis(Val(2),VT,order,Polynomials._p_filter)
+  fb = MonomialBasis(Val(1),T,order,Polynomials._p_filter)
+  cb = MonomialBasis(Val(2),VT,order-1,Polynomials._p_filter)
 
   function cmom(φ,μ,ds) # Cell and Node moment function: σ_K(φ,μ) = ∫(φ:μ)dK
     Broadcasting(Operation(⊙))(φ,μ)