changing default quadrature for tensor product elements to N+1 pt Gauss

jlchan · Aug 3, 2024 · f086cac · f086cac
1 parent 6a47a47
commit f086cac
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/RefElemData_polynomial.jl b/src/RefElemData_polynomial.jl
@@ -10,7 +10,7 @@ RefElemData(elem::Union{Line, Tri, Tet, Wedge, Pyr},
 
 # on quad and hex elements, default to a tensor product quadrature 
 RefElemData(elem::Union{Quad, Hex}, approximation_type::Polynomial{DefaultPolynomialType}, N; kwargs...) = 
-    RefElemData(elem, Polynomial(TensorProductQuadrature(gauss_quad(0, 0, N+1))), N; kwargs...)
+    RefElemData(elem, Polynomial(TensorProductQuadrature(gauss_quad(0, 0, N))), N; kwargs...)
 
 # special case: for lines, tensor product and multidimensional quadrature are the same
 RefElemData(elem::Line, approx_type::Polynomial{<:TensorProductQuadrature}, N; kwargs...) =