ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bug fix in original read_dump for add keep atom count

Browse Source
This commit is contained in:
Steve Plimpton
2018-11-29 15:54:14 -07:00
parent a2671f8943
commit 4e9b987fec
1 changed files with 8 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/read_dump.cpp
View File

@ -1050,12 +1050,17 @@ void ReadDump::process_atoms()
} }
} }
// if addflag = YESADD // 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 consistent IDs to new snapshot atoms across all procs // assign consistent IDs to new snapshot atoms across all procs
if (addflag == YESADD) { 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) if (atom->natoms < 0 || atom->natoms >= MAXBIGINT)
error->all(FLERR,"Too many total atoms"); error->all(FLERR,"Too many total atoms");
if (atom->tag_enable) atom->tag_extend(); if (atom->tag_enable) atom->tag_extend();
@ -1069,7 +1074,6 @@ void ReadDump::process_atoms()
memory->destroy(updateflag); memory->destroy(updateflag);
memory->destroy(newflag); memory->destroy(newflag);
} }
/* ---------------------------------------------------------------------- /* ----------------------------------------------------------------------