| ... | ... | @@ -125,6 +125,48 @@ time = 48:00:00 |
|
|
|
|
|
|
|
<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><summary>Plots</summary>
|
| ... | ... | @@ -135,11 +177,32 @@ time = 48:00:00 |
|
|
|
## Scaling test
|
|
|
|
|
|
|
|
Ran the above [Fiducial 4s BBH Inj](#fiducial-bbh-inj) with nCPU 32, 64, 128.
|
|
|
|
|
|
|
|
<details><summary>Plots</summary>
|
|
|
|
|
|
|
|
|
|
|
|
| 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"> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
</details>
|
| ... | ... | |
