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

correctly handle the case where there are no atoms in the fix group

Browse Source
This commit is contained in:
Axel Kohlmeyer
2022-08-12 10:46:07 -04:00
parent 12c28c393d
commit 09b3cbafdc
1 changed files with 11 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
src/REAXFF/fix_reaxff_species.cpp
View File

@ -26,6 +26,7 @@
#include "error.h"
#include "fix_ave_atom.h"
#include "force.h"
#include "group.h"
#include "memory.h"
#include "modify.h"
#include "neigh_list.h"
@ -395,9 +396,7 @@ void FixReaxFFSpecies::Output_ReaxFF_Bonds(bigint ntimestep, FILE * /*fp*/)
Nmole = Nspec = 0;
FindMolecule();
SortMolecule(Nmole);
FindSpecies(Nmole, Nspec);
vector_nmole = Nmole;
@ -522,16 +521,21 @@ void FixReaxFFSpecies::SortMolecule(int &Nmole)
hi = MAX(hi, nint(clusterID[n]));
}
int flagall;
MPI_Allreduce(&flag, &flagall, 1, MPI_INT, MPI_SUM, world);
if (flagall && me == 0)
error->warning(FLERR, "Atom with cluster ID = 0 included in fix reaxff/species group");
MPI_Allreduce(&lo, &idlo, 1, MPI_INT, MPI_MIN, world);
MPI_Allreduce(&hi, &idhi, 1, MPI_INT, MPI_MAX, world);
int nlen = idhi - idlo + 1;
if (nlen <= 0) { // no atoms in group
Nmole = 0;
return;
}
if (idlo == ntotal)
if (me == 0)
error->warning(FLERR, "Atom with cluster ID = maxmol included in fix reaxff/species group");
error->warning(FLERR, "Atom with cluster ID = maxmol included in fix reaxff/species group {}",group->names[igroup]);
MPI_Allreduce(&flag, &flagall, 1, MPI_INT, MPI_SUM, world);
if (flagall && me == 0)
error->warning(FLERR, "Atom with cluster ID = 0 included in fix reaxff/species group {}", group->names[igroup]);
int nlen = idhi - idlo + 1;
memory->create(molmap, nlen, "reaxff/species:molmap");
for (n = 0; n < nlen; n++) molmap[n] = 0;