Skip to content

Update tnt.py#330

Draft
jmv2009 wants to merge 21 commits intomscroggs:mainfrom
jmv2009:patch-1
Draft

Update tnt.py#330
jmv2009 wants to merge 21 commits intomscroggs:mainfrom
jmv2009:patch-1

Conversation

@jmv2009
Copy link
Contributor

@jmv2009 jmv2009 commented Jan 2, 2026

Add 1st order 0-Form elements
Add interval, triangle, tetrahedron, prism and pyramid for 0-Forms, modified from Cockburn & Fu
Modified the internal degrees of freedom using integration by parts to Laplacian.
Shifted order of curl for matching with 0-Forms in complex

  • There appears to be a problem for 3rd order prism on the backend.

Todo:

  • Hexahedron 1st
  • Eliminate curl kernel of 3d elements (tetrahedron already works?)
  • Pyramid (visualization is defelement issue)
  • Hdiv
  • DPC pyramid for inclusion in complex (serendipity type)

Add 1st order 0-Form elements
Add prism and pyramid for 0-Forms, modified from Cockburn & Fu
Fix minors
clean up
@jmv2009 jmv2009 marked this pull request as draft January 3, 2026 08:37
@jmv2009
Copy link
Contributor Author

jmv2009 commented Jan 17, 2026

Complexes should now be virtually complete.
Lowest order may not always work yet.
Pyramid required DPC for highest forms (3-Forms), others discontinuous tensor product.

continuity = "L2"
value_type = "scalar"
last_updated = "2023.07.1"
last_updated = "2026.1"
Copy link
Owner

Choose a reason for hiding this comment

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

while you're developing, you probably want to add cache = False here to stop it locally caching the element matrices

@jmv2009
Copy link
Contributor Author

jmv2009 commented Jan 28, 2026

Thanks! (commit some high order stuff pending).

@jmv2009
Copy link
Contributor Author

jmv2009 commented Jan 29, 2026

Got all elements up to k=4 complex (lagrange=5).

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.

2 participants