Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Some assembly refactors #26269

Merged
merged 23 commits into from Dec 19, 2023
Merged

Conversation

lindsayad
Copy link
Member

Work done in support of issue #25666 and its PR #26073

@lindsayad lindsayad force-pushed the assembly-refactors-25666 branch 2 times, most recently from 12f8844 to 9a5907f Compare December 7, 2023 02:09
@lindsayad lindsayad self-assigned this Dec 7, 2023
@moosebuild
Copy link
Contributor

moosebuild commented Dec 7, 2023

Job Documentation on 7d80d46 wanted to post the following:

View the site here

This comment will be updated on new commits.

@lindsayad
Copy link
Member Author

@rui-hu @lingzou @lingzou-anl the tests that are exodiffing were previously using an incorrect calculation of JxW. Could someone check and see what you think of the new results?

@moosebuild
Copy link
Contributor

Job Coverage on 7d80d46 wanted to post the following:

Framework coverage

d450a5 #26269 7d80d4
Total Total +/- New
Rate 85.01% 85.01% - 100.00%
Hits 97252 97252 - 48
Misses 17153 17153 - 0

Diff coverage report

Full coverage report

Modules coverage

Coverage did not change

Full coverage reports

Reports

This comment will be updated on new commits.

@lindsayad lindsayad marked this pull request as ready for review December 8, 2023 16:06
@lindsayad
Copy link
Member Author

@roystgnr I'd most like your review on 7d80d46. I think I would likely add a sentence to the assertion that says "you really should just be integrating with a single JxW, so for DG you should be using the JxW from the element side" (with "element" being MOOSE jargon for "the governing element for an across-face element pair")

@loganharbour
Copy link
Member

Let me know if you need a review

@lindsayad
Copy link
Member Author

I do need one 😄

@lingzou
Copy link
Contributor

lingzou commented Dec 11, 2023

@lindsayad it is great that @loganharbour can help with the review. this MR is a bit beyond my capabilities.

@lindsayad
Copy link
Member Author

It would be good for a SAM dev to consider the SAM failures though and whether what you are now getting for output files reflect better solutions

@ke7kto
Copy link
Contributor

ke7kto commented Dec 13, 2023

This PR fixes issue #26328

@lingzou
Copy link
Contributor

lingzou commented Dec 16, 2023

It would be good for a SAM dev to consider the SAM failures though and whether what you are now getting for output files reflect better solutions

@lindsayad I am trying to understand what needs to be done on SAM side, like what has been changed in this MR, and what are needed in SAM; or maybe there is nothing needed in SAM, but just to make sure the failed cases should be updated and to accept the changes.

@lingzou
Copy link
Contributor

lingzou commented Dec 16, 2023

@lindsayad It seems to me that this MR may have introduced a bug somewhere. Please refer to one of the failed SAM test case: tests/components/SurfaceCoupling/GapHT.i

Note that this problem has analytical solutions (please see the test head of that problem). This MR produced results very far away from the analytical solutions, so this makes me think maybe

  1. there was a bug in this MR; or
  2. SAM should be updated following this MR;

unless I was completely wrong about the analytical solutions.

@lindsayad
Copy link
Member Author

It did not introduce a bug. It fixed a bug. I think that previously the JxW returned by the Assembly API was half of what it should have been. When FE::reinit is called with a point argument but not a weights argument, then dummy weights of 1 are used. However, for one Gauss point for integration, the weight should be 2.

@lindsayad
Copy link
Member Author

If this caused a diff from an analytical or reference solution, then it seems like maybe SAM was working around the bug

@lindsayad
Copy link
Member Author

On line 81 of SAMNodeFaceConstraint.C you multiply JxW by 2. Perhaps this was SAM working around the bug

@lindsayad lindsayad merged commit c3a7159 into idaholab:next Dec 19, 2023
43 of 44 checks passed
@lindsayad lindsayad deleted the assembly-refactors-25666 branch December 19, 2023 15:26
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
Development

Successfully merging this pull request may close these issues.

ADParsedFunctorMaterial evaluation attempting to access invalid qt elements
6 participants