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

Allow string periodic for periodic BCs #1347

Merged
merged 1 commit into from
Nov 15, 2024
Merged

Conversation

marchdf
Copy link
Contributor

@marchdf marchdf commented Nov 15, 2024

Summary

Allow string periodic for periodic BCs.

Basically it's weird to get an error message that says:

std::runtime_error: Setting periodic BC on a periodic face xlo is not allowed

Pull request type

Please check the type of change introduced:

  • Bugfix
  • Feature
  • Code style update (formatting, renaming)
  • Refactoring (no functional changes, no api changes)
  • Build related changes
  • Documentation content changes
  • Other (please describe):

Checklist

The following is included:

  • new unit-test(s)
  • new regression test(s)
  • documentation for new capability

This PR was tested by running:

  • the unit tests
    • on GPU
    • on CPU
  • the regression tests
    • on GPU
    • on CPU

@marchdf marchdf requested a review from mbkuhn November 15, 2024 16:09
@mbkuhn
Copy link
Contributor

mbkuhn commented Nov 15, 2024

so, the geom.is_periodic still has to be correct, but this now allows the user to put a periodic argument for that boundary face instead of having to remove the boundary face argument. What happens if a user puts periodic for the face but has geom.is_periodic = 0 for that direction?

@marchdf
Copy link
Contributor Author

marchdf commented Nov 15, 2024

the geom.is_periodic still has to be correct, but this now allows the user to put a periodic argument for that boundary face instead of having to remove the boundary face argument.

yes or (currently) it has to be xlo.type = null, which is weird.

What happens if a user puts periodic for the face but has geom.is_periodic = 0 for that direction?

If I do xlo.type=periodic geometry.is_periodic=0 1 0 then it aborts with:

Invalid BC specification for non-periodic boundary = periodic

Copy link
Contributor

@mbkuhn mbkuhn left a comment

Choose a reason for hiding this comment

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

works for me!

@marchdf marchdf merged commit 9ecbeae into Exawind:main Nov 15, 2024
14 of 15 checks passed
@marchdf marchdf deleted the tweak-bc-check branch November 15, 2024 17:15
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