Skip to content

Implement reduced quintic triangle element#89

Open
Sichao25 wants to merge 5 commits into
SCOREC:mainfrom
Sichao25:yus/quintic
Open

Implement reduced quintic triangle element#89
Sichao25 wants to merge 5 commits into
SCOREC:mainfrom
Sichao25:yus/quintic

Conversation

@Sichao25

Copy link
Copy Markdown
Contributor

The reduced quintic triangle element defines 6 DOFs per vertex. The element uses a local coordinate system rather than barycentric coordinates. The evaluation introduces a set of coefficients defined by the geometric parameters of each element.

}

// Solve A*X = B using our LU solver (row-major)
int info = solveLU_internal(20, 18, &A[0][0], 20, &B[0][0], 18);
@Sichao25 Sichao25 marked this pull request as ready for review June 12, 2026 21:52
@Sichao25

Copy link
Copy Markdown
Contributor Author

@jacobmerson @Joshua-Kloepfer This is ready for view when you have a chance.

Comment thread src/MeshField_Element.hpp
ReducedQuinticTriangleShape>) {
assert(elemCoeffs.data() != nullptr);
auto coeffSlice = Kokkos::subview(elemCoeffs, ent, Kokkos::ALL());
return shapeFn.getValues(localCoord, coeffSlice);

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@cwsmith @Joshua-Kloepfer the way we are doing this may warrant some discussion since the getValues requires additional geometric information.

This may come down to how we want to handle/formalize non-isoparametric elements.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants