[Users] POWER code issues

[Roland Haas]

Hello Anuj,

Thank you or your email.

Those sound like bugs indeed. Could you create an ET ticket at 

https://bitbucket.org/einsteintoolkit/tickets/issues/new

and I will try and see if I can incorporate the fixes you suggested
quickly.

Can you provide your replacement fitting code that you refer to in:

> 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.

Yours,
Roland

> 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



-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://keys.gnupg.net.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
URL: <http://lists.einsteintoolkit.org/pipermail/users/attachments/20230831/e9bf0590/attachment-0001.sig>