From 088d870e2062e7613b60f317eaccccd1dfff41cf Mon Sep 17 00:00:00 2001 From: jrgissing Date: Tue, 30 May 2023 17:17:54 -0400 Subject: [PATCH] check remaining per-atom initializations --- src/atom.cpp | 13 ++++++++----- 1 file changed, 8 insertions(+), 5 deletions(-) diff --git a/src/atom.cpp b/src/atom.cpp index d5f8b64a52..3e9e177ae0 100644 --- a/src/atom.cpp +++ b/src/atom.cpp @@ -2131,13 +2131,16 @@ void Atom::add_molecule_atom(Molecule *onemol, int iatom, int ilocal, tagint off for (int i = 0; i < nivector; ++i) if (ivname[i] != nullptr) ivector[i][ilocal] = 0; - for (int i = 0; i < ndvector; ++i) dvector[i][ilocal] = 0.0; + for (int i = 0; i < ndvector; ++i) + if (dvname[i] != nullptr) dvector[i][ilocal] = 0.0; for (int i = 0; i < niarray; ++i) - for (int j = 0; j < icols[i]; ++j) - iarray[i][ilocal][j] = 0; + if (ianame[i] != nullptr) + for (int j = 0; j < icols[i]; ++j) + iarray[i][ilocal][j] = 0; for (int i = 0; i < ndarray; ++i) - for (int j = 0; j < dcols[i]; ++j) - darray[i][ilocal][j] = 0.0; + if (daname[i] != nullptr) + for (int j = 0; j < dcols[i]; ++j) + darray[i][ilocal][j] = 0.0; if (molecular != Atom::MOLECULAR) return;