Fix XYZ trajectories - prevent final frame with 0,0,0 coords

Thanks to brockdyer03 https://discuss.avogadro.cc/t/opening-trajectory-files/5570/11 Signed-off-by: Geoff Hutchison <[email protected]>
OpenChemistry · Sep 29, 2024 · 93893d4 · 93893d4
1 parent 8563422
commit 93893d4
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/avogadro/io/xyzformat.cpp b/avogadro/io/xyzformat.cpp
@@ -138,6 +138,7 @@ bool XyzFormat::read(std::istream& inStream, Core::Molecule& mol)
     getline(inStream, buffer); // Skip the blank
     mol.setCoordinate3d(mol.atomPositions3d(), 0);
     int coordSet = 1;
+    bool done = false;
     while (numAtoms == numAtoms2) {
       Array<Vector3> positions;
       positions.reserve(numAtoms);
@@ -146,6 +147,7 @@ bool XyzFormat::read(std::istream& inStream, Core::Molecule& mol)
         getline(inStream, buffer);
         if (inStream.eof()) {
           numAtoms2 = 0;
+          done = true;
           break; // break this inner loop
         }
 
@@ -160,7 +162,8 @@ bool XyzFormat::read(std::istream& inStream, Core::Molecule& mol)
         positions.push_back(pos);
       }
 
-      mol.setCoordinate3d(positions, coordSet++);
+      if (!done)
+        mol.setCoordinate3d(positions, coordSet++);
 
       if (getline(inStream, buffer)) {
         if (inStream.eof()) {