... | ... | @@ -286,7 +286,7 @@ for param in pesummary_samples.keys(): |
|
|
|
|
|
# Check posterior samples
|
|
|
not_included = []
|
|
|
for param in ["tilt1_isco", "tilt2_isco", "phi12_isco", "l_peak_evol_avg", "mf_evol_avg", "e_rad_evol_avg", "af_evol_avg"]:
|
|
|
for param in ["l_peak_evol_avg", "mf_evol_avg", "e_rad_evol_avg", "af_evol_avg"]:
|
|
|
if param in pesummary_dict.keys() and param != "e_rad_evol_avg":
|
|
|
comparison = np.round(pesummary_dict[param], 8) == np.round(lalinf_samples[param], 8)
|
|
|
difference = np.abs(pesummary_dict[param] - lalinf_samples[param])
|
... | ... | @@ -301,9 +301,6 @@ for param in ["tilt1_isco", "tilt2_isco", "phi12_isco", "l_peak_evol_avg", "mf_e |
|
|
With output:
|
|
|
|
|
|
```bash
|
|
|
Maximum difference for tilt1_isco: 0.00016935480523350144
|
|
|
Maximum difference for tilt2_isco: 0.00032114999367438557
|
|
|
Maximum difference for phi12_isco: 0.0008060677213324041
|
|
|
Maximum difference for l_peak_evol_avg: 1.0976890889136115e-05
|
|
|
Maximum difference for mf_evol_avg: 2.745710679619151e-05
|
|
|
Maximum difference for af_evol_avg: 9.770286555355057e-07
|
... | ... | |