From c2356d11af23a6c94551abc66f5eb508c9e90ae5 Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Wed, 20 Oct 2021 15:47:05 +0000
Subject: [PATCH] update parity waveform without small phase expansion
---
waveform/wrapper_mpv_nosmallexpansion/mpvnosmallwaveform.py | 4 ++--
waveform/wrapper_mpv_nosmallexpansion/setup.py | 2 +-
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/waveform/wrapper_mpv_nosmallexpansion/mpvnosmallwaveform.py b/waveform/wrapper_mpv_nosmallexpansion/mpvnosmallwaveform.py
index 2bc6315..4b0622c 100644
--- a/waveform/wrapper_mpv_nosmallexpansion/mpvnosmallwaveform.py
+++ b/waveform/wrapper_mpv_nosmallexpansion/mpvnosmallwaveform.py
@@ -27,8 +27,8 @@ def gen(**kwds):
zz = cosmology.redshift(kwds['distance'])
intz = integrate.quad(integrand_parityamu_mpvinverse, 0, zz)[0]
temp = kwds['parity_mpvinverse'] * intz / 1e9 / lal.QE_SI * (lal.H_SI / 2 / lal.PI) * lal.PI * lal.PI / lal.H0_SI
- expminus = np.exp(-1j*temp*hp.sample_frequencies**2)
- expplus = 1/expplus
+ expminus = numpy.exp(-1j*temp*hp.sample_frequencies**2)
+ expplus = 1/expminus
hp_parity = (hp+1j*hc)*expminus/2 + (hp-1j*hc)*expplus/2
hc_parity = (hp+1j*hc)*expminus/2j - (hp-1j*hc)*expplus/2j
diff --git a/waveform/wrapper_mpv_nosmallexpansion/setup.py b/waveform/wrapper_mpv_nosmallexpansion/setup.py
index 2e0c8ff..bcfd286 100644
--- a/waveform/wrapper_mpv_nosmallexpansion/setup.py
+++ b/waveform/wrapper_mpv_nosmallexpansion/setup.py
@@ -16,7 +16,7 @@ setup (
url = 'http://www.pycbc.org/',
#download_url = 'https://github.com/gwastro/revchirp/tarball/v%s' % VERSION,
keywords = ['pycbc', 'signal processing', 'gravitational waves'],
- py_modules = ['mpvnosmall'],
+ py_modules = ['mpvnosmallwaveform'],
entry_points = {"pycbc.waveform.fd":"mpvnosmall=mpvnosmallwaveform:gen"},
classifiers=[
'Programming Language :: Python',
--
GitLab