[Einstein Toolkit] #949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ---------------------------------------------+------------------------------ Reporter: bmundim | Owner: Type: defect | Status: new Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Keywords: RotatingSymmetry90 Weyl scalars | ---------------------------------------------+------------------------------ RotatingSymmetry90 is aborting when evaluating tensor type for Weyl scalars. It doesn't take into account the ManualCartesian type for example. The attached patch avoids that.
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ------------------------------------+--------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: review Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: | Keywords: RotatingSymmetry90 Weyl scalars ------------------------------------+--------------------------------------- Changes (by bmundim):
* status: new => review
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ------------------------------------+--------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: closed Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: invalid | Keywords: RotatingSymmetry90 Weyl scalars ------------------------------------+--------------------------------------- Changes (by bmundim):
* status: review => closed * resolution: => invalid
Comment:
On a second thought, I should actually comment out the Weyl scalars from my parameter file, since it doesn't make any sense to measure them in spherical symmetry!
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ------------------------------------+--------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: closed Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: invalid | Keywords: RotatingSymmetry90 Weyl scalars ------------------------------------+---------------------------------------
Comment (by rhaas):
I am not quite sure if one does not have to be a bit more careful with ManualCartesian quantities. If I define Psi = x*y^2 then this would have ManualCartesian parities of (-,+,+). If I understand RotatingSymmetry90 correctly it would copy the value psi(i=10,j=1,k=0) to point (i=-1,j=10,k=0). However the Psi I defined would have (among other changes) changed sign as x (or i) changes sign. So a simple copy would not seem to be correct.
My understanding of RotatingSymmetry is limited and that of Weyl scalars very spotty at best. So don't take this objection too serious :-).
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ------------------------------------+--------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: closed Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: invalid | Keywords: RotatingSymmetry90 Weyl scalars ------------------------------------+---------------------------------------
Comment (by hinder):
The manual Cartesian parities specify the sign changes of the function under ''reflections'' in the planes perpendicular to the axes. These can therefore be used directly by ReflectonSymmetry. With RotatingSymmetry180, the symmetry operation is fortuitously equivalent to a reflection in the x direction composed with a reflection in the y direction, so if you know how a function transforms under reflections in the axes, you can determine how it transforms under the 180 degree rotation. RotatingSymmetry180 contains the logic to apply the parities that have been chosen. However, I don't think that a 90 degree rotation can be written as a composition of reflections, so the manual Cartesian parities are not useful here. In any case, RotatingSymmetry90 does not take account of these parities.
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ------------------------------------+--------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: closed Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: invalid | Keywords: RotatingSymmetry90 Weyl scalars ------------------------------------+---------------------------------------
Comment (by eschnett):
This patch only removes a check, it does not implement actually looking at these parities.
It does make sense to look at them -- e.g. Psi_2 contains valuable information about mass and spin also in spherical symmetry.
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ----------------------+----------------------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: reopened Priority: major | Milestone: Component: Cactus | Version: Resolution: | Keywords: RotatingSymmetry90 ----------------------+----------------------------------------------------- Changes (by tbode):
* keywords: RotatingSymmetry90 Weyl scalars => RotatingSymmetry90 * resolution: invalid => * status: closed => reopened * component: EinsteinToolkit thorn => Cactus
Comment:
This is a duplicate of #884, which I will now close as this thread has more content. While it may not be of use to you in your specific parameter file and project, this is not strictly invalid. WeylScal4 calculates Weyl scalars besides Psi4, which may be of interest to those studying systems for which RotatingSymmetry90 is applicable. Regardless, a symmetry thorn should support all tensortypealiases available. I am reopening, but downgrading to minor as I expect this won't be of interest to many and it does abort with a helpful error. I'll get around to it when I can.
#949: RotatingSymmetry90: abort when evaluating tensor type for Weyl scalars ----------------------+----------------------------------------------------- Reporter: bmundim | Owner: Type: defect | Status: reopened Priority: major | Milestone: Component: Cactus | Version: Resolution: | Keywords: RotatingSymmetry90 ----------------------+-----------------------------------------------------
Comment (by bmundim):
Thanks for pointing that out! I see that it has its value for the other Weyl scalars other than Psi4. Erik: thanks for pointing out the issue with the patch. I guess I rushed into it and didn't actually implement the correct solution. I will take a deeper look later and see if I can come up with a better patch.
Cheers!
trac@lists.einsteintoolkit.org
-
Einstein Toolkit