-
Notifications
You must be signed in to change notification settings - Fork 213
Add DG FreeOutflow boundary condition to GRMHD #4196
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
Add DG FreeOutflow boundary condition to GRMHD #4196
Conversation
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! A few minor things you can squash immediately :)
| * | ||
| * - However, regardless of whether the normal component of the spatial | ||
| * velocity $n_iv^i$ is pointing outward or inward, the lorentz factor \f$W\f$ | ||
| * is copied into ghost zone without any modification. |
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.
Maybe write out (briefly) the justification we had for this?
| const Scalar<DataVector>& interior_lapse, | ||
| const tnsr::II<DataVector, 3, Frame::Inertial>& | ||
| interior_inv_spatial_metric) { | ||
| const size_t vector_size = get(interior_rest_mass_density).size(); |
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.
vector_size -> number_of_grid_points
| auto& exterior_divergence_cleaning_field = | ||
| get<::Tags::TempScalar<2>>(temp_buffer); | ||
|
|
||
| determinant_and_inverse(make_not_null(&sqrt_det_spatial_metric), |
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.
Is it not possible to get these two from the interior? If not, that's fine, just curious :)
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.
Turns out.. it can 😂
I pushed this change as a fixup since it comes with nontrivial change of codes
b9c2625 to
259243c
Compare
|
Squashed fixes to the first two comments and pushed a fixup for the third one. |
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.
Looks good, one comment, you can squash that in right away :)
tests/Unit/Evolution/Systems/GrMhd/ValenciaDivClean/BoundaryConditions/Test_FreeOutflow.cpp
Show resolved
Hide resolved
259243c to
72c61e0
Compare
Proposed changes
Upgrade instructions
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