add pmandel.h
Rewrite neutron-para with MPI_Reduce (by Marjorie)
FIXME: root isn't always zero
Add max_tasks to steal_work. It doesn't work (cf l. 284)