Skip to content

Commit 62cd345

Browse files
volkerspringelstephenpcook
authored andcommitted
small bug fix: in case HIERARCHICAL_GRAVITY is not used, and the maximum used timestep sizes increases during a step, it could happen that for particles on the maximum timestep one gravity half-step is not applied (because HighestActiveTimeBin increases)
1 parent 285b495 commit 62cd345

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

src/amuse/community/arepo/src/time_integration/do_gravity_hydro.c

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -265,7 +265,7 @@ void find_gravity_timesteps_and_do_gravity_step_first_half(void)
265265
}
266266

267267
/* reconstruct list of active particles because it is used for other things too (i.e. wind particles) */
268-
timebin_make_list_of_active_particles_up_to_timebin(&TimeBinsGravity, All.HighestActiveTimeBin);
268+
timebin_make_list_of_active_particles_up_to_timebin(&TimeBinsGravity, All.HighestSynchronizedTimeBin);
269269
sumup_large_ints(1, &TimeBinsGravity.NActiveParticles, &TimeBinsGravity.GlobalNActiveParticles);
270270
#else /* #ifdef HIERARCHICAL_GRAVITY */
271271

@@ -276,7 +276,7 @@ void find_gravity_timesteps_and_do_gravity_step_first_half(void)
276276
timebin_make_list_of_active_particles_up_to_timebin(&TimeBinsGravity, TIMEBINS);
277277
else
278278
#endif /* #ifdef FORCE_EQUAL_TIMESTEPS */
279-
timebin_make_list_of_active_particles_up_to_timebin(&TimeBinsGravity, All.HighestActiveTimeBin);
279+
timebin_make_list_of_active_particles_up_to_timebin(&TimeBinsGravity, All.HighestSynchronizedTimeBin);
280280
sumup_large_ints(1, &TimeBinsGravity.NActiveParticles, &TimeBinsGravity.GlobalNActiveParticles);
281281

282282
mpi_printf("KICKS: 1st gravity for highest active timebin=%d: particles %lld\n", All.HighestActiveTimeBin,

0 commit comments

Comments
 (0)