Change numerical evidence to log space.

Author: alicjapolanska

examples/rosenbrock.py

@@ -205,7 +205,7 @@ def run_example(
     """
     Set up and run multiple simulations
     """
-    n_realisations = 1
+    n_realisations = 100
     ln_evidence_inv_summary = np.zeros((n_realisations, 5))
     for i_realisation in range(n_realisations):
         if n_realisations > 1:
@@ -453,7 +453,7 @@ def run_example(
     if n_realisations > 1:
         save_name = (
             save_name_start
-            + "_rosenbrock_evidence_inv_T"
+            + "_rosenbrock_evidence_log_inv_T"
             + str(temperature)
             + "_realisations.dat"
         )
@@ -464,8 +464,8 @@ def run_example(
 
         if ndim == 2:
             evidence_inv_analytic_summary = np.zeros(1)
-            evidence_inv_analytic_summary[0] = 1.0 / evidence_numerical_integration
-            save_name = save_name_start + "_rosenbrock_evidence_inv" + "_analytic.dat"
+            evidence_inv_analytic_summary[0] = -np.log(evidence_numerical_integration)
+            save_name = save_name_start + "_rosenbrock_evidence_log_inv" + "_analytic.dat"
             np.savetxt(
                 save_name,
                 evidence_inv_analytic_summary,