Remnant-fits.md

@@ -276,6 +276,26 @@ and we run the following `summarypages` executable on commit 4eb7e4a (where para
 $ summarypages --webdir /home/charlie.hoy/public_html/projects/pesummary_review/remnant_fits/parallel_evolved/parallel --samples /home/charlie.hoy/projects/pesummary_review/remnant_fits/parallel_evolved/PhenomPv3HM_samples.hdf5 --evolve_spins --approximant IMRPhenomPv3HM --labels parallel --multi_process 10
 ```
 
+We then compare the samples generated by the two runs and compare highlight the maximum difference:
+
+```python
+import numpy as np
+
+no_parallel = np.genfromtxt("/home/charlie.hoy/public_html/projects/pesummary_review/remnant_fits/parallel_evolved/no_parallel/samples/no_parallelisation_pesummary.dat", names=True)
+parallel = np.genfromtxt("/home/charlie.hoy/public_html/projects/pesummary_review/remnant_fits/parallel_evolved/parallel/samples/parallel_pesummary.dat", names=True)
+
+assert len(set(no_parallel.dtype.names) - set(parallel.dtype.names)) == 0
+for i in no_parallel.dtype.names:
+    print("{} max difference: {}".format(
+        i, np.max(np.abs(no_parallel[i] - parallel[i]))
+    ))
+```
+
+with output:
+
+```bash
+
+```
 
 ## Review comments