Repository navigation
Add edge and wedge quadrature with separate triangle and zeta orders - #118
Draft
abhiyanpaudel wants to merge 41 commits into
Draft
abhiyanpaudel wants to merge 41 commits into
abhiyanpaudel wants to merge 41 commits into
Conversation
Add Wedge to Mesh_Topology and numWedge to MeshInfo for wedge (triangular prism) elements. Wedge is appended after Pyramid so no existing enum value changes; Prism and numPrism are unchanged.
LinearWedgeShape is the linear wedge (triangular prism): the tensor product of the linear triangle in (xi0, xi1) and the linear edge in zeta in [-1, 1]. Nodes 0-2 are the triangle vertices at zeta = -1 and nodes 3-5 those at zeta = +1, following Zienkiewicz, Taylor and Zhu, Section 6.4, Fig. 6.16a. test/testWedgeShape.cpp checks the Kronecker delta property at the nodes and, at interior points, partition of unity, reproduction of the bilinear functions spanned by the basis and their gradients, and the analytic gradients against finite differences. test/testWedgeElement.cpp builds a FieldElement from a coordinate field, LinearWedgeShape and a small wedge-to-vertex map, and checks MeshField::evaluate, getJacobians and getJacobianDeterminants on a right prism and a skewed wedge against the geometric map written out by hand.
# Conflicts: # CMakeLists.txt
abhiyanpaudel
marked this pull request as draft
October 27, 2026 05:50
This branch has not been deployed
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Adds quadrature for
Wedgeelements.EdgeIntegration(N1–N6): Gauss-Legendre rules on [-1, 1] copied from SCOREC/coreapfIntegrate.cc; also available asgetIntegration<Edge>().getWedgeIntegrationPoints(triOrder, zetaOrder): triangle rule × edge rule; exact forgetIntegrationPoints<topo>(order, zetaOrder): optional zeta order for wedges (defaults toorder), so existing calls are unchanged.Integrator(triOrder, zetaOrder)constructor, used byprocess()for wedges;Integrator(order)is unchanged.WedgeIntegrationTest: edge and wedge rule exactness up to (3, 11) (36 points), andIntegrator::processvolumes of a right prism and a tapered wedge.