[Users] Bug in QuasiLocalMeasures

[Marco Brito]

Good afternoon,

I found a bug in the QuasiLocalMeasures thorn, in the computation of the ADM momentum. In the source code file qlm_analyse.F90, in line 361 there is a plus sign instead of a minus sign, contradicting the ADM momentum formula: P^i_{ADM} = \frac{1}{8\pi}\lim_{r\to\infty} \oint_S (K^i_j - \delta^i_j K)dS^j [1]

As expected for slices where K \neq 0  one was getting a wrong result. The correction is easy, just replace the plus sign by a minus sign. In any case you can find a patch file in attachment. Should I in any case file a ticket or a Pull Request somewhere regarding this?


Best regards,

Marco Brito.

[1] M. Alcubierre, Introduction to 3+1 Numerical Relativity, International Series of Monographs on Physics (Oxford University Press, Oxford, 2008).

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20260402/566674a6/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: qlm_patch
Type: application/octet-stream
Size: 512 bytes
Desc: qlm_patch
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20260402/566674a6/attachment.obj>