Implement reduced quintic triangle element#89
Conversation
![github-advanced-security[bot]](https://avatars.githubusercontent.com/in/57789?s=60&v=4){kind=small}
| } | ||
|
|
||
| // Solve A*X = B using our LU solver (row-major) | ||
| int info = solveLU_internal(20, 18, &A[0][0], 20, &B[0][0], 18); |
|
@jacobmerson @Joshua-Kloepfer This is ready for view when you have a chance. |
| ReducedQuinticTriangleShape>) { | ||
| assert(elemCoeffs.data() != nullptr); | ||
| auto coeffSlice = Kokkos::subview(elemCoeffs, ent, Kokkos::ALL()); | ||
| return shapeFn.getValues(localCoord, coeffSlice); |
There was a problem hiding this comment.
@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.
cwsmith
left a comment
cwsmith
left a comment
There was a problem hiding this comment.
Hi @Sichao25,
Thanks for adding this.
#91 was just merged into main and changed the way parametric coordinates are handled and their ordering for triangles and tets. I'm guessing that will break things here, but I haven't dug deep to confirm.
Some comments are below from a quick scan of the code.
A few of the comments are just reminders for you, @jacobmerson, @Joshua-Kloepfer, and I to discuss.
4 participants
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.