Re: Conflict with FOF and Starformation

From: Volker Springel <vspringel_at_MPA-Garching.MPG.DE>
Date: Thu, 15 Apr 2021 14:36:24 +0200

Dear Leonard,

Thanks for pointing out this problem. I agree with your analysis. The fix to simply call domain_rearrange_particle_sequence() can have some problematic side effects in some situations however (it invalidates the SPH neighbor search tree, but this may still be used before its next reconstruction). I think it should be safe to simply do a full domain decomposition before a snapshot/FOF is done when star formation is active.

I have added such thing in lines 525-530 of src/main/ If this doesn't cure the problem for you, please speak up.


> On 12. Apr 2021, at 18:22, Leonard Romano <> wrote:
> Dear Gadget list members,
> I am running cosmological simulations with star formation and experience crashes whenever gas particles are converted into stars in the same step as a snapshot is written. In particular as soon as "particle_exchange_based_on_PS" is called it crashes with the error "Bummer2".
> My guess would be that the particle list has not been reordered yet as it would be right before a domain decomposition and so a call to "domain_rearrange_particle_sequence" seems to be missing.
> I propose promoting "domain_rearrange_particle_sequence" to a public method and calling it right before the I/O in order to avoid similar problems.
> Best,
> Leonard
> --
> ============================================
> Leonard Romano, B.Sc.(レオナルド・ロマノ)
> Physics Department
> Technical University of Munich (TUM), Germany
> Theoretical Astrophysics Group
> Department of Earth and Space Science
> Graduate School of Science, Osaka University, Japan
> he / him / his
> ============================================
> -----------------------------------------------------------
> If you wish to unsubscribe from this mailing, send mail to
> with a subject of: unsubscribe gadget-list
> A web-archive of this mailing list is available here:
Received on 2021-04-15 14:36:24

This archive was generated by hypermail 2.3.0 : 2023-01-10 10:01:32 CET