tx258/lammps-gran-kokkos
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bug-fix for miscounting atoms for read_dump add keep

Browse Source
This commit is contained in:
Steve Plimpton
2018-11-29 15:37:58 -07:00
parent a4d21370bf
commit f92e61c58d
1 changed files with 7 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/read_dump.cpp
View File

@ -480,11 +480,16 @@ void ReadDump::atoms()
nread += nchunk;
}
// if addflag = YESADD or KEEPADD, update total atom count
if (addflag == YESADD || addflag == KEEPADD) {
bigint nblocal = atom->nlocal;
MPI_Allreduce(&nblocal,&atom->natoms,1,MPI_LMP_BIGINT,MPI_SUM,world);
}
// if addflag = YESADD, assign IDs to new snapshot atoms
if (addflag == YESADD) {
bigint nblocal = atom->nlocal;
MPI_Allreduce(&nblocal,&atom->natoms,1,MPI_LMP_BIGINT,MPI_SUM,world);
if (atom->natoms < 0 || atom->natoms >= MAXBIGINT)
error->all(FLERR,"Too many total atoms");
if (atom->tag_enable) atom->tag_extend();