-
Notifications
You must be signed in to change notification settings - Fork 213
Only interpolate for horizon finds from certain blocks #6611
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
Only interpolate for horizon finds from certain blocks #6611
Conversation
d36f635 to
9008704
Compare
| "outflow-type boundary condition, you must use that."); | ||
| } | ||
| using domain::BoundaryConditions::is_periodic; | ||
| if (boundary_condition_ != nullptr and is_periodic(boundary_condition_)) { |
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.
Only the important domains ;) :D
src/ParallelAlgorithms/ApparentHorizonFinder/InterpolationTarget.cpp
Outdated
Show resolved
Hide resolved
src/ParallelAlgorithms/ApparentHorizonFinder/InterpolationTarget.hpp
Outdated
Show resolved
Hide resolved
src/ParallelAlgorithms/ApparentHorizonFinder/InterpolationTarget.hpp
Outdated
Show resolved
Hide resolved
src/ParallelAlgorithms/Interpolation/Actions/InitializeInterpolator.hpp
Outdated
Show resolved
Hide resolved
|
@nilsdeppe Pushed fixups |
nilsdeppe
left a comment
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.
LGTM, please squash :)
| std::make_index_sequence< | ||
| std::tuple_size_v<std::decay_t<decltype(all_horizon_options)>>>{}); | ||
| expand_pack(append_to_result( | ||
| tmpl::at<option_tags<Metavariables>, tmpl::size_t<Is + 1>>::name(), |
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.
you can use at_c and then don't need the tmpl::size_t
These are currently taken as options from the input file. Eventually we could have this be more automatic, e.g. the domain creators know which blocks could potentially hold a horizon, but that would be a future optimization
86d81d2 to
3d620a2
Compare
Proposed changes
These are currently taken as options from the input file. Eventually
we could have this be more automatic, e.g. the domain creators know
which blocks could potentially hold a horizon, but that would be
a future optimization.
This is the first step in a lot of backend work to combine the horizon finders and get rid of the interpolator group component.
I've tested this on a single BH while interpolating only some blocks and all blocks. And also tried on a BBH
Upgrade instructions
Each entry under the
ApparentHorizons:block of the input file gets a new optionBlocksForInterpolation:. This could either be a list of block groups orAllto indicate all blocks in the doman.Code review checklist
make docto generate the documentation locally intoBUILD_DIR/docs/html.Then open
index.html.code review guide.
bugfixornew featureif appropriate.Further comments