When running Palabos programs in parallel, there seems to be little or no mass transfer from one periodic domain to the next. For the shanchenMultiPhase example, segregation2D.cpp, all four quadrants of a four processor simulation have identical density distributions at every timestep. I end up with one droplet in each quadrant. For a droplet spreading simulation, the droplet will not spread into the next processor domain. Do you have any ideas as to what is causing this?
I get identical results using either mpich2 or Intel’s impi library. I’m compiling with MPIparallel = true and I do see a substantial speed up, so the mpi implementation seems to be working. I also compared v1.0 with v1.1, turned off optimization, and toggling SMPparallel with no change in results.