Skip to content

Refspace fixes - #4538

Open
roystgnr wants to merge 7 commits into
libMesh:develfrom
roystgnr:refspace_fixes
Open

Refspace fixes#4538
roystgnr wants to merge 7 commits into
libMesh:develfrom
roystgnr:refspace_fixes

Conversation

@roystgnr

@roystgnr roystgnr commented Sep 2, 2026

Copy link
Copy Markdown
Member

In one of @rochi00's branches his new test coverage discovered some utterly broken master point coordinate definitions for Pyramid13, Pyramid14, and especially Pyramid18. I fixed those, added some tests for the fix, and it turns out Tri7/Prism20/Prism21 also had issues. All fixed and in test coverage now.

I don't think this is a slow one but it's worth logging.
I'm not sure how we had the wrong data here.
Rochi found bugs in our higher-order Pyramids reference elements, and
these tests catch them.
@moosebuild

Copy link
Copy Markdown

Job Coverage, step Generate coverage on c3bcaaa wanted to post the following:

Coverage

0a0a9d #4538 c3bcaa
Total Total +/- New
Rate 65.94% 65.91% -0.02% 100.00%
Hits 79392 79272 -120 3
Misses 41014 40997 -17 0

Diff coverage report

Full coverage report

This comment will be updated on new commits.

@roystgnr

roystgnr commented Sep 2, 2026

Copy link
Copy Markdown
Member Author

I think we may have a threading race condition problem with ReferenceElem::get(), despite doing a bunch of work in reference_elem.C to avoid that...

@roystgnr roystgnr mentioned this pull request Sep 2, 2026
Comment thread src/geom/cell_pyramid.C
Comment on lines +47 to +55
{-0.5, -0.5, 0.5},
{0.5, -0.5, 0.5},
{0.5, 0.5, 0.5},
{-0.5, 0.5, 0.5},
{0, 0, 0},
{0, -Real(2)/3, Real(1)/3},
{ Real(2)/3, 0, Real(1)/3},
{0, Real(2)/3, Real(1)/3},
{-Real(2)/3, 0, Real(1)/3}

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Hmm... I don't think I realized that the base class's _master_points array was meant to have entries for all possible subtypes that might exist? Maybe we should document that better somehow?

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