[ET Trac] [Einstein Toolkit] #254: tensorparity=-1 quantities in ReflectionSymmetry
Einstein Toolkit
trac-noreply at einsteintoolkit.org
Mon Jan 31 20:12:12 CST 2011
#254: tensorparity=-1 quantities in ReflectionSymmetry
--------------------------------------------+-------------------------------
Reporter: roland.haas@… | Type: defect
Status: new | Priority: major
Milestone: | Component: Cactus
Version: | Keywords:
--------------------------------------------+-------------------------------
I would like to suggest applying the attached patch to ReflectionSymmetry
to fix (I believe) its handling of the tensorparity attribute.
Explanation/Rationale:
In ReflectionSymmetry/interpolate.c parities[dir] is set to the behaviour
under reflection across the "dir"=0 plane of an object of type
tensortypealis. It does not take tensorparity into account. Therefore
check_dir[dir] must take tensorparity into account explicitly.
The second change is "parity = 1" instead of "parity = tensorparity" which
states that if a point is never reflected then the value that was on the
grid is to be used as-is.
I have tested these changes with the magnetic field created by a wire
along the x,y,z axis and reflection across the various planes and with
RotatingSymmetry180 and *90. Please see ticket #149
https://trac.einsteintoolkit.org/ticket/149
for sample parameter files and test scripts.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/254>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
More information about the Trac
mailing list