parallelize numberDensity() with MPIWhirlPool