Update O4 Review authored by Avi Vajpeyi's avatar Avi Vajpeyi
...@@ -125,6 +125,48 @@ time = 48:00:00 ...@@ -125,6 +125,48 @@ time = 48:00:00
<details><summary>Configs Contents</summary> <details><summary>Configs Contents</summary>
[Link to config](https://git.ligo.org/avi.vajpeyi/pbilby_o4_review/-/blob/main/injection_test/scaling_test/128_cpus_v1.ini)
```ini
nlive=2000
nact=50
sampling-seed=0
gaussian-noise=True
detectors=['H1', 'L1']
duration=4.0
generation-seed=1010
injection=True
injection-dict={'chirp_mass': 17.051544979894693, 'mass_ratio': 0.3183945489993522, 'a_1': 0.29526500202350264, 'a_2': 0.23262056301313416, 'tilt_1': 1.0264673717225983, 'tilt_2': 2.1701305583885513, 'phi_12': 5.0962562029664955, 'phi_jl': 2.518241237045709, 'luminosity_distance': 497.2983560174788, 'dec': 0.2205292600865073, 'ra': 3.952677097361719, 'theta_jn': 1.8795187965094322, 'psi': 2.6973435044499543, 'phase': 3.686990398567503, 'geocent_time':1238175800}
distance-marginalization=True
phase-marginalization=False
time-marginalization=True
jitter-time=True
trigger-time = 0
reference-frame=H1L1
time-reference=geocent
default-prior=BBHPriorDict
prior-dict={
chirp-mass: bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=14, maximum=20),
mass_ratio=bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1),
mass_1=Constraint(name='mass_1', minimum=1.001398, maximum=1000),
mass_2=Constraint(name='mass_2', minimum=1.001398, maximum=1000),
a_1 = Uniform(name='a_1', minimum=0, maximum=0.99),
a_2 = Uniform(name='a_2', minimum=0, maximum=0.99),
tilt_1 = Sine(name='tilt_1'),
tilt_2 = Sine(name='tilt_2'),
phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic'),
phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic'),
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3, unit='Mpc'),
dec: Cosine(name='dec'), ra: Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic'),
cos-theta-jn: Uniform(minimum=-1, maximum=1, name='cos_theta_jn'),
psi: Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic'),
phase: Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic'),
geocent_time = Uniform(minimum=-0.1, maximum=0.1, name=None, latex_label=None, unit=None, boundary=None)
}
waveform-approximant = IMRPhenomXPHM
```
</details> </details>
<details><summary>Plots</summary> <details><summary>Plots</summary>
...@@ -135,11 +177,32 @@ time = 48:00:00 ...@@ -135,11 +177,32 @@ time = 48:00:00
## Scaling test ## Scaling test
Ran the above [Fiducial 4s BBH Inj](#fiducial-bbh-inj) with nCPU 32, 64, 128. Ran the above [Fiducial 4s BBH Inj](#fiducial-bbh-inj) with nCPU 32, 64, 128.
<details><summary>Plots</summary> <details><summary>Plots</summary>
![scaling_1d_histograms](uploads/150257ce0d59c5211b6ef36476b0df68/scaling_1d_histograms.png) | ncpu | time | nsamples | log_evidence | log_evidence_err |
|-------:|--------:|-----------:|---------------:|-------------------:|
| 32 | 57307.8 | 14587 | -16332.8 | 0.108911 |
| 32 | 58007 | 14534 | -16332.8 | 0.108911 |
| 64 | 30241.6 | 15305 | -16332.8 | 0.109552 |
| 64 | 30356.7 | 15144 | -16332.8 | 0.109552 |
| 128 | 17208.8 | 15196 | -16332.7 | 0.109678 |
| 128 | 17072.3 | 15278 | -16332.7 | 0.109678 |
| Plots |
| ------ |
| <img src="uploads/c7e58df9f02943cc674e51337e018c61/scaling_metadata.png" height="500"> |
| <img src="uploads/81612c80fde5285764e850f795fd0e84/scaling_1d_histograms.png" height="1000"> |
![scaling_metadata](uploads/e4c9a6e86665457b85e9003ad269f5ad/scaling_metadata.png)
</details> </details>
... ...
......