283
282
if (rmass) sbuf[3] = rmass[j];
284
283
else sbuf[3] = mass[type[j]];
285
284
MPI_Sendrecv(sbuf,4,MPI_DOUBLE,all[1].proc,0,
286
rbuf,4,MPI_DOUBLE,all[1].proc,0,world,&status);
285
rbuf,4,MPI_DOUBLE,all[1].proc,0,world,MPI_STATUS_IGNORE);
287
286
vcm[0] = (sbuf[3]*sbuf[0] + rbuf[3]*rbuf[0]) / (sbuf[3] + rbuf[3]);
288
287
vcm[1] = (sbuf[3]*sbuf[1] + rbuf[3]*rbuf[1]) / (sbuf[3] + rbuf[3]);
289
288
vcm[2] = (sbuf[3]*sbuf[2] + rbuf[3]*rbuf[2]) / (sbuf[3] + rbuf[3]);