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

Fix a check for zero entries in the lumped mass matrix. #13211

Merged
merged 1 commit into from Jan 10, 2022

Conversation

drwells
Copy link
Member

@drwells drwells commented Jan 10, 2022

This is unrelated to #13179 but I noticed it while examining one of the failing tests (threading with oneAPI causes vectors to apparently be set to all zeros, which creates a floating point exception).

It doesn't make sense to check the values of inverse_lumped_diagonal here - that's what we are computing. The check should be on the values of the lumped mass matrix.

@drwells
Copy link
Member Author

drwells commented Jan 10, 2022

/rebuild

Copy link
Member

@kronbichler kronbichler left a comment

Choose a reason for hiding this comment

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

Great to fix this oversight. I notice that with the current code we should only rarely run into the if branch here because we already set the constrained entries to one before.

@bangerth bangerth merged commit d5995ec into dealii:master Jan 10, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

6 participants