Add code to compute tidal field in parallel for cosmoDC2

LSSTDESC / DC2_ia_mocks

Repository storing tools and scripts used to populate protoDC2 and DC2 with realistic alignments

MIT License

0 stars 1 forks source link

Add code to compute tidal field in parallel for cosmoDC2 #1

Open EiffL opened 6 years ago

EiffL commented 6 years ago

This issue is to cover the development of the tools necessary to compute the tidal field for DC2:

Bin a list of particles into a 3D array in parallel
Compute the Fourier transform
Compute tidal field
Inverse Fourier transform
Diagonalize locally
Return cube, or interpolate at galaxy positions

EiffL commented 5 years ago

@patricialarsen I'm reviving this issue to track the development of a tidal field code for DC2. I'll give a go using pfft https://github.com/rainwoodman/pfft-python

EiffL commented 5 years ago

Tidal field computation distributed in MPI is working! Thanks to @patricialarsen !

Now in the process of making run on cori to extract the tidal field at galaxy positions.

Code for this is being developed here: https://github.com/LSSTDESC/DC2_ia_mocks/tree/u/EiffL/tidal

EiffL commented 5 years ago

Hum... I don't understand, it's not working when trying to read the particles for the outter rim snapshot. I can submit MPI jobs on cori which work just fine to completion using the small particle file you sent me. But my jobs keep getting killed when loading the real snapshot. I've tried up to 16 Haswell nodes, I would assume this enough to hold the particles in memory.