[Users] POWER code issues

[Anuj Kankani]

Hi,

I was recently suggested the POWER (
https://git.ncsa.illinois.edu/elihu/Gravitational_Waveform_Extractor/) code
for some GW extrapolation work I was doing by Gabriele Bozzola but I
believe I found some bugs and he suggested I send an email here.

Running the script as is on the GW150914 gallery example gives me results
significantly different from what the plots on the ET website show. I found
replacing where the code currently fits the strain amplitude and phase to a
polynomial expression with numpy polyfit seems to fix the main problem.

I also noticed in the getCutoffFrequencyFromTwoPuncturesBBH function, lines
341 - 347 have an incorrect equation for angular momentum, it looks like a
typo between position and momentum. Right now its position *  momentum -
momentum^2.

I'm not sure if this code is being actively used, but I saw there were git
branches for the latest ET release.

Thanks,
Anuj Kankani
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20230830/1782bdf3/attachment.htm>