-
Notifications
You must be signed in to change notification settings - Fork 90
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
Update outer boundaries in PETSc 3D solver #2496
Conversation
clang-tidy review says "All clean, LGTM! 👍" |
Only affects boundary conditions if MXG > 1: First set of boundary cells are already set as part of the solve. Changes from always zero-gradient to extrapolating. If zero-gradient boundary conditions are set then the result will not change.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks @bendudson! This looks totally general, and consistent with boundary_standard.cxx
(with suggestions below). Nice how simple it is!
Is this one of things that can only be effectively tested with MMS? |
Probably. It's tricky because the extra boundary cells affected by this PR are outside the stencil of the forward operator, so aren't checked by the existing integrated tests. |
Co-authored-by: johnomotani <john.omotani@ukaea.uk>
Co-authored-by: johnomotani <john.omotani@ukaea.uk>
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me!
Only affects boundary conditions if MXG > 1: First set of boundary cells are already set as part of the solve.
Changes from always zero-gradient to extrapolating. If zero-gradient boundary conditions are set then the result will not change.