diff --git a/src/L2_projection.jl b/src/L2_projection.jl index 69717fc699..31ed414a08 100644 --- a/src/L2_projection.jl +++ b/src/L2_projection.jl @@ -326,6 +326,13 @@ function assemble_proj_rhs!(f::Matrix, cellvalues::CellValues, sdh::SubDofHandle end end +""" + evaluate_at_grid_nodes(proj::L2Projector, vals::AbstractVector) +Return a vector of length `getnnodes(grid)` where the order corresponds with the node order +in the `grid` used to create `proj`. + +`vals` must have length `getnnodes(grid)`. +""" evaluate_at_grid_nodes(proj::L2Projector, vals::AbstractVector) = _evaluate_at_grid_nodes(proj, vals, Val(false)) diff --git a/src/Quadrature/quadrature.jl b/src/Quadrature/quadrature.jl index 27379bddfa..af6c708b43 100644 --- a/src/Quadrature/quadrature.jl +++ b/src/Quadrature/quadrature.jl @@ -15,8 +15,8 @@ using Base.Cartesian: @nloops, @ntuple, @nexprs QuadratureRule{shape}(weights::AbstractVector{T}, points::AbstractVector{Vec{rdim, T}}) Create a `QuadratureRule` used for integration on the refshape `shape` (of type [`AbstractRefShape`](@ref)). -`order` is the order of the quadrature rule. -`quad_rule_type` is an optional argument determining the type of quadrature rule, +`order` is related to the number of quadrature points per direction (i.e. a `QuadratureRule{RefQuadrilateral}(2)` +uses 2 quadrature points per direction). `quad_rule_type` is an optional argument determining the type of quadrature rule, currently the `:legendre` and `:lobatto` rules are implemented for hypercubes. For triangles up to order 8 the default rule is the one by `:dunavant` (see [Dun:1985:hde](@cite)) and for tetrahedra the default rule is `keast_minimal` (see [Keast:1986:mtq](@cite)). Wedges and pyramids default