Changes
Page history
Update O4 Review
authored
Mar 30, 2023
by
Avi Vajpeyi
Show whitespace changes
Inline
Side-by-side
O4-Review.md
View page @
8fe0a240
...
...
@@ -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>
...
...
...
...
