276
275
if (rmass) sbuf[1] = rmass[ipos];
277
276
else sbuf[1] = mass[type[ipos]];
278
277
MPI_Sendrecv(sbuf,2,MPI_DOUBLE,all[1].proc,0,
279
rbuf,2,MPI_DOUBLE,all[1].proc,0,world,&status);
278
rbuf,2,MPI_DOUBLE,all[1].proc,0,world,MPI_STATUS_IGNORE);
280
279
vcm = (sbuf[1]*sbuf[0] + rbuf[1]*rbuf[0]) / (sbuf[1] + rbuf[1]);
281
280
v[ipos][vdim] = 2.0 * vcm - sbuf[0];
282
281
pswap += sbuf[1] * (vcm - sbuf[0]);