Added check for null boxes in link cells

plumed · Jan 20, 2025 · e7a1618 · e7a1618
1 parent c86000e
commit e7a1618
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/tools/LinkCells.cpp b/src/tools/LinkCells.cpp
@@ -48,6 +48,10 @@ void LinkCells::buildCellLists( const std::vector<Vector>& pos, const std::vecto
 
   // Must be able to check that pbcs are not nonsensical in some way?? -- GAT
 
+  double determinant = pbc.getBox().determinant();
+
+  plumed_assert(determinant > epsilon) <<"Cell lists cannot be built when passing a box with null volume. Volume is "<<determinant;
+
   // Setup the pbc object by copying it from action
   mypbc.setBox( pbc.getBox() );