Better mesh size and maximal time-step size for all the EAGLE examples.

examples/EAGLE_100/eagle_100.yml

@@ -23,7 +23,7 @@ TimeIntegration:
   dt_max:     1e-3  # The maximal time-step size of the simulation (in internal units).
 
 Scheduler:
-  max_top_level_cells: 80
+  max_top_level_cells: 32
 
 # Parameters governing the snapshots
 Snapshots:
@@ -44,6 +44,7 @@ Gravity:
   theta:                  0.85      # Opening angle (Multipole acceptance criterion)
   comoving_softening:     0.0026994 # Comoving softening length (in internal units).
   max_physical_softening: 0.0007    # Physical softening length (in internal units).
+  mesh_side_length:       256
 
 # Parameters for the hydrodynamics scheme
 SPH:

examples/EAGLE_12/eagle_12.yml

@@ -20,7 +20,7 @@ TimeIntegration:
   time_begin: 0.    # The starting time of the simulation (in internal units).
   time_end:   1e-2  # The end time of the simulation (in internal units).
   dt_min:     1e-10 # The minimal time-step size of the simulation (in internal units).
-  dt_max:     1e-3  # The maximal time-step size of the simulation (in internal units).
+  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).
   
 # Parameters governing the snapshots
 Snapshots:

examples/EAGLE_25/eagle_25.yml

@@ -30,10 +30,7 @@ TimeIntegration:
   time_begin: 0.    # The starting time of the simulation (in internal units).
   time_end:   1e-2  # The end time of the simulation (in internal units).
   dt_min:     1e-10 # The minimal time-step size of the simulation (in internal units).
-  dt_max:     1e-3  # The maximal time-step size of the simulation (in internal units).
-
-Scheduler:
-  max_top_level_cells:    16
+  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).
 
 # Parameters governing the snapshots
 Snapshots:
@@ -54,7 +51,7 @@ Gravity:
   theta:                  0.7     # Opening angle (Multipole acceptance criterion)
   comoving_softening:     0.0026994 # Comoving softening length (in internal units).
   max_physical_softening: 0.0007    # Physical softening length (in internal units).
-  mesh_side_length:       32
+  mesh_side_length:       64
 
 # Parameters for the hydrodynamics scheme
 SPH:

examples/EAGLE_50/eagle_50.yml

@@ -20,10 +20,10 @@ TimeIntegration:
   time_begin: 0.    # The starting time of the simulation (in internal units).
   time_end:   1e-2  # The end time of the simulation (in internal units).
   dt_min:     1e-10 # The minimal time-step size of the simulation (in internal units).
-  dt_max:     1e-3  # The maximal time-step size of the simulation (in internal units).
+  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).
 
 Scheduler:
-  max_top_level_cells: 20
+  max_top_level_cells: 16
 
 # Parameters governing the snapshots
 Snapshots:
@@ -44,7 +44,7 @@ Gravity:
   theta:                  0.7      # Opening angle (Multipole acceptance criterion)
   comoving_softening:     0.0026994 # Comoving softening length (in internal units).
   max_physical_softening: 0.0007    # Physical softening length (in internal units).
-  mesh_side_length:       64
+  mesh_side_length:       128
 
 # Parameters for the hydrodynamics scheme
 SPH:

src/mesh_gravity.c

@@ -291,7 +291,8 @@ void pm_mesh_compute_potential(struct pm_mesh* mesh, const struct space* s,
   const double dim[3] = {s->dim[0], s->dim[1], s->dim[2]};
 
   if (r_s <= 0.) error("Invalid value of a_smooth");
-  if (mesh->dim[0] != dim[0] || mesh->dim[1] != dim[1] || mesh->dim[2] != dim[2])
+  if (mesh->dim[0] != dim[0] || mesh->dim[1] != dim[1] ||
+      mesh->dim[2] != dim[2])
     error("Domain size does not match the value stored in the space.");
 
   /* Some useful constants */