search

bccp / nbodykit

Analysis kit for large-scale structure datasets, the massively parallel way
http://nbodykit.rtfd.io
GNU General Public License v3.0
111 stars 60 forks source link

fix to meshtools to avoid nans in legpoly at k=0 #622

Closed akrolewski closed 4 years ago

akrolewski commented 4 years ago

At kx=ky=kz=0, in slab.mu(los), the norm2 is 0 and thus mu is nan. This leads to a nan in legpoly which propagates into nans for all multipoles in the bin containing k=0 for multipoles other than the monopole. following suggestion of @michaelJwilson, modify SlabIterator.mu in meshtools to fix this issue.

rainwoodman commented 4 years ago

Thanks!