-
Notifications
You must be signed in to change notification settings - Fork 213
Add Ricci Scalar compute tags to EvolveGeneralizedHarmonic.hpp #3576
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
Conversation
|
@markscheel could you take an expert-look at this one? :) |
| GeneralizedHarmonic::Tags::Phi<volume_dim, frame>, | ||
| Tags::deriv<GeneralizedHarmonic::Tags::Phi<volume_dim, frame>, | ||
| tmpl::size_t<3>, frame>>; |
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.
(edited) Earlier I thought this Tag should go away, but I see now that it should remain. I was confused about the frame. Note that this will only work if frame is Frame::Inertial.
I think everything in this type alias should have Frame::Inertial, even if the horizon is found in a different frame.
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.
Sounds good, so I'm okay to leave this as is? I tried to compile without it but I kept getting an error that said it was required if I was trying to call the Ricci scalar tag.
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.
Yes, please leave this as is, except please change frame to Frame::Inertial in interpolator_source_vars. Let me know if you get an error if you do that.
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.
Okay, that compiled
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.
Why wouldn't the tag work for other frames? @markscheel
| StrahlkorperTags::ExtrinsicCurvatureCompute<frame>, | ||
| StrahlkorperGr::Tags::SpinFunctionCompute<frame>>, | ||
| StrahlkorperGr::Tags::SpinFunctionCompute<frame>, | ||
| StrahlkorperTags::RicciScalarCompute<frame>>, |
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.
(edited) This looks ok. This is the function that computes the (2D) Ricci scalar of the Strahlkorper, right? Not the (3D) Ricci scalar of the volume.
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.
I'm not sure. @geoffrey4444, would you be able to clarify if this computes (2D) or (3D)?
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.
Yes, it's 2D
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.
Ok, sounds good.
|
@MarloMo I took another look at this; I was confused by names and frames. Just one thing, which is changing everything to Inertial in interpolator_source_vars. |
74c7227 to
e54bcad
Compare
markscheel
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.
Please squash.
e54bcad to
3c3991c
Compare
Proposed changes
Adding the Ricci scalar compute tags to
EvolveGeneralizedHarmonic.hppthat computes the Min and Max Ricci scalar of a Strahlkorper. This PR also adds arguments tags that are dependent on calculating the Ricci scalar quantities.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