diff --git a/examples/EvrardCollapse_3D/evrard.yml b/examples/EvrardCollapse_3D/evrard.yml
index 72543ae40c5a886faddf44563948b54b7ff635c8..db696230448683d1b8cbee6706999a7eda8feabc 100644
--- a/examples/EvrardCollapse_3D/evrard.yml
+++ b/examples/EvrardCollapse_3D/evrard.yml
@@ -32,14 +32,10 @@ SPH:
 Gravity:
   eta:                   0.025    # Constant dimensionless multiplier for time integration.
   epsilon:               0.001      # Softening length (in internal units).
-  theta:                 0.1
+  theta:                 0.9
   a_smooth:              1.25     # (Optional) Smoothing scale in top-level cell sizes to smooth the long-range forces over (this is the default value).
   r_cut:                 4.5      # (Optional) Cut-off in number of top-level cells beyond which no FMM forces are computed (this is the default value).
 
 # Parameters related to the initial conditions
 InitialConditions:
   file_name:  ./evrard.hdf5       # The file to read
-
-Scheduler:
-  tasks_per_cell: 400
-  cell_split_size:       20 # or so
diff --git a/examples/EvrardCollapse_3D/makeIC.py b/examples/EvrardCollapse_3D/makeIC.py
index 15886f88eaf3faf75d607d984e05f8aa62354303..beb133cca8ee3f6a15e6ed5d202b07199a682ff4 100644
--- a/examples/EvrardCollapse_3D/makeIC.py
+++ b/examples/EvrardCollapse_3D/makeIC.py
@@ -45,7 +45,7 @@ pos[:,2] = r * cos_theta
 
 pos += array([R, R, R])
 
-h = ones(numPart) * R / numPart**(1. / 3.)
+h = ones(numPart) * 2. * R / numPart**(1. / 3.)
 
 # Generate extra arrays
 v = zeros((numPart, 3))