Skip to content
Snippets Groups Projects
Commit c2356d11 authored by Yifan Wang's avatar Yifan Wang
Browse files

update parity waveform without small phase expansion

parent e9a1ee98
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
......@@ -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',
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment