diff --git a/examples/EAGLE_12/eagle_12.yml b/examples/EAGLE_12/eagle_12.yml index d2192b68e1bee7721a892d832daa4cd0dadb3836..f56c330590ac25cc5b3fe8f68ed68aa1e94d6490 100644 --- a/examples/EAGLE_12/eagle_12.yml +++ b/examples/EAGLE_12/eagle_12.yml @@ -26,7 +26,7 @@ Statistics: # Parameters for the self-gravity scheme Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. - epsilon: 0.0000001 # Softening length (in internal units). + epsilon: 0.001 # Softening length (in internal units). theta: 0.7 # Opening angle (Multipole acceptance criterion) # Parameters for the hydrodynamics scheme diff --git a/src/gravity.c b/src/gravity.c index de42d58f87cb420de0e9ad9f3e47750568571d17..f3082622744e4d3bb5b1b3270029a3f3b84a5663 100644 --- a/src/gravity.c +++ b/src/gravity.c @@ -436,7 +436,7 @@ void gravity_exact_force_compute_mapper(void *map_data, int nr_gparts, a_grav[2] += f * dz; /* Apply Ewald correction for periodic BC */ - if (periodic) { + if (periodic && r > 1e-5 * hi) { double corr[3]; gravity_exact_force_ewald_evaluate(dx, dy, dz, corr);