- removed comment
BVec in HydroBase/interface.ccl should be an axial vector
Issue #149
resolved
the magnetic field is a pseudovector so should have tensorparity=-1 in interface.ccl.
CCTK_REAL Bvec[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" tensorparity=-1 interpolator="matter"' "Magnetic field components B^i"
Parities are used by the Reflection thorn.
Keyword:
Comments (6)
-
-
- changed status to open
- assigned issue to
- removed comment
Roland, can you create an apply the corresponding patch?
-
reporter - removed comment
provided patch and test case. Please have a look.
-
- changed status to open
- assigned issue to
- removed comment
-
reporter - removed comment
Replying to [comment:4 knarf]: Frank: the patch to interface.ccl only was applied (by me, taking Erik's comment2 and the discussion we had about this on the mailing list as a review) in r36. The tests require patches to ReflectionSymmetry (ticket
#254). -
reporter - changed status to resolved
- removed comment
- Log in to comment
Yes, this is one way to specify the symmetries of the magnetic fields, as axial vector. tensorparity=-1 is implemented and should be supported, but since magnetic fields are one of the first use cases this may require some careful debugging. Please add some test cases to ensure this is working.
The other way would be to treat it as antisymmetric tensor, which is also defined by 3 of its components. However, I believe that antisymmetric tensors are not yet supported by the symmetry thorns -- but that could be added.