pelahi / VELOCIraptor-STF

Galaxy/(sub)Halo finder for N-body simulations
MIT License
19 stars 26 forks source link

Fix for ordering of particles when returning Velociraptor group IDs to Swift #96

Closed jchelly closed 4 years ago

jchelly commented 4 years ago

When running on the fly from Swift and writing Velociraptor group IDs to the snapshot, the resulting dataset is scrambled. This happens because the particles are sorted by binding energy but the corresponding pfof array which stores the group membership is not reordered.

With this change the VELOCIraptorGroupIDs array in Swift outputs is correct when running on one MPI rank. On multiple ranks VELOCIraptorGroupIDs does not agree with the velociraptor output files so there's at least one more bug left.