`iread_gather_on_root` in `ConvolutionDensityMatricesFromWaveFunctions` hangs when executing in MPI. Probably because the same method (for `ConvolutionDensityMatricesFromFrequency`) is called by the main iterator.