BVec in HydroBase/interface.ccl should be an axial vector

Create issue
Issue #149 closed
Roland Haas created an issue

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.


Comments (7)

  1. Erik Schnetter
    • removed 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.

  2. Erik Schnetter
    • changed status to open
    • assigned issue to
    • removed comment

    Roland, can you create an apply the corresponding patch?

  3. Roland Haas 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).

  4. Log in to comment