文件列表:

Makefile.am (792, 2023-12-02)
__init__.py (27150, 2023-12-02)
config.h (3888, 2023-12-02)
configure.ac (1200, 2023-12-02)
ctype_structs.py (36499, 2023-12-02)
figures/ (0, 2023-12-02)
figures/example.jpg (161396, 2023-12-02)
figures/example2.jpg (175835, 2023-12-02)
interface.c (60314, 2023-12-02)
main.c (37438, 2023-12-02)
mcmc.py (43583, 2023-12-02)
mcmc_bilby.py (5405, 2023-12-02)
mkconf.sh (398, 2023-12-02)
mode_collector.py (331, 2023-12-02)
myAlloc.c (3338, 2023-12-02)
myAlloc.h (775, 2023-12-02)
myDatatypes.h (5300, 2023-12-02)
myFileIO.c (19356, 2023-12-02)
myFileIO.h (924, 2023-12-02)
myLog.c (9806, 2023-12-02)
myLog.h (2430, 2023-12-02)
myOptparser.c (17327, 2023-12-02)
myOptparser.h (697, 2023-12-02)
mySpherical.c (25408, 2023-12-02)
mySpherical.h (1184, 2023-12-02)
myUtils.c (32067, 2023-12-02)
myUtils.h (7555, 2023-12-02)
newFactorizedWaveform.c (112162, 2023-12-02)
newFactorizedWaveform.h (3692, 2023-12-02)
newFactorizedWaveformPrec.c (113410, 2023-12-02)
newFactorizedWaveformPrec.h (3052, 2023-12-02)
pBHRingdown.c (290950, 2023-12-02)
pBHRingdown.h (3701, 2023-12-02)
pCore.c (403104, 2023-12-02)
pCore.h (28328, 2023-12-02)
pDebug.c (6787, 2023-12-02)
pDebug.h (170, 2023-12-02)
pEnergyFlux.c (61614, 2023-12-02)
... ...

# pyWaveformGenerator calculate GW waveform Installation && compile ``` ./mkconf.sh ./configure make ``` sometimes the lib cannot use in Apple core macbook. Run in python ```python from pyWaveformGenerator.waveform import calculate_waveform m1 = m2 = 10 chi1x = chi1y = chi1z = 0 chi2x = chi2y = chi2z = 0 e0 = 0 dL = 100 zeta_rad = iota_rad = beta_rad = Phic_rad = 0 fMin = 20 waveform, dynamics = calculate_waveform((m1, m2, chi1x, chi1y, chi1z, chi2x, chi2y, chi2z, e0, dL, zeta_rad, iota_rad, beta_rad, Phic_rad), fMin) ``` see more details and description in waveform.py Cite this repo via [arXiv:2102.08614](https://arxiv.org/abs/2102.08614) and [arXiv:2310.04552](https://arxiv.org/abs/2310.04552). Generate eccentric-precession waveform, run ```python from pyWaveformGenerator.waveform import calculate_waveform_ep m1 = m2 = 10 chi1x = 0.2 chi1y = chi1z = 0 chi2y = -0.3 chi2x = chi2z = 0 e0 = 0.2 dL = 100 zeta_rad = 0 iota_rad = beta_rad = Phic_rad = 0 fMin = 20 Mf_ref = 0.003 waveform, dynamics = calculate_waveform_ep((m1, m2, chi1x, chi1y, chi1z, chi2x, chi2y, chi2z, e0, dL, zeta_rad, iota_rad, beta_rad, Phic_rad), fMin, Mf_ref = Mf_ref, srate=16384, is_coframe=False) ``` In this case one have to make sure that one of $\chi_{1x,1y,2x,2y}\neq0$. Then we can see the waveform ```Python import matplotlib.pyplot as plt h22 = waveform.h22 fig = plt.figure(figsize = (12, 6)) ax1 = fig.add_subplot(3, 1, 1) ax2 = fig.add_subplot(3, 1, 2) ax3 = fig.add_subplot(3, 1, 3) ax1.plot(h22.time, h22.real) ax1.set_ylabel(r'$\Re[h_{22}]$') ax2.plot(h22.time, h22.amp) ax2.set_ylabel(r'${\rm Amp}[h_{22}]$') ax3.plot(h22.time, h22.phase) ax3.set_ylabel(r'${\rm Phase}[h_{22}]$') ax3.set_xlabel(r'$t[M]$') plt.show() ``` which will looks like  ## Anomaly and reference frequency ```python from pyWaveformGenerator.waveform import calculate_waveform m1 = m2 = 10 chi1x = chi1y = chi1z = 0 chi2x = chi2y = chi2z = 0 e0 = 0.3 dL = 100 iota_rad = beta_rad = Phic_rad = 0 zeta_rad = 30 * np.pi/180 fMin = 20 Mf_ref = 0.003 waveform, dynamics = calculate_waveform((m1, m2, chi1x, chi1y, chi1z, chi2x, chi2y, chi2z, e0, dL, zeta_rad, iota_rad, beta_rad, Phic_rad), fMin, Mf_ref = Mf_ref) ``` for spin-precession case (on test) ```python from pyWaveformGenerator.waveform import calculate_waveform_ep q = 1.3 mtot = 20 m1 = q * mtot / (1. + q) m2 = mtot / (1. + q) chi1x = 0.4 chi1y = chi1z = 0 chi2y = -0.3 chi2x = chi2z = 0 e0 = 0.2 dL = 100 zeta_rad = 60 * np.pi/180 iota_rad = beta_rad = Phic_rad = 0 fMin = 20 Mf_ref = 0.002 waveform, dynamics = calculate_waveform_ep((m1, m2, chi1x, chi1y, chi1z, chi2x, chi2y, chi2z, e0, dL, zeta_rad, iota_rad, beta_rad, Phic_rad), fMin, Mf_ref = Mf_ref, srate=16384, is_coframe=False) ``` The waveform looks like ```python import matplotlib.pyplot as plt h22 = waveform.h22 h21 = waveform.h21 fig = plt.figure(figsize = (12, 6)) ax1 = fig.add_subplot(3, 1, 1) ax2 = fig.add_subplot(3, 1, 2) ax3 = fig.add_subplot(3, 1, 3) ax1.plot(h22.time, h22.real) ax1.set_ylabel(r'$\Re[h_{22}]$') ax2.plot(h21.time, h21.real) ax2.set_ylabel(r'$\Re[h_{21}]$') ax3.plot(h22.time, h22.amp) ax3.plot(h21.time, h21.amp) ax3.set_ylabel(r'${\rm Amp}[h_{2m}]$') ax3.set_xlabel(r'$t[M]$') plt.savefig(pwd / 'example2.jpg',dpi = 200) plt.show() ``` 

近期下载者：

相关文件：

评论：[我要评论] [举报此文件]

收藏者：