min_image fix; Saves long simulations from crashing; Addresses https://matsci.org/t/fix-rigid-atoms-have-moved-too-far-apart/53926

PolyachenkoYA · jhenin · commit 151869f8857b · 2024-12-16T16:14:15.000+01:00
diff --git a/lammps/src/COLVARS/colvarproxy_lammps.cpp b/lammps/src/COLVARS/colvarproxy_lammps.cpp
@@ -205,7 +205,7 @@ cvm::rvector colvarproxy_lammps::position_distance(cvm::atom_pos const &pos1,
   double xtmp = pos2.x - pos1.x;
   double ytmp = pos2.y - pos1.y;
   double ztmp = pos2.z - pos1.z;
-  _lmp->domain->minimum_image(xtmp,ytmp,ztmp);
+  _lmp->domain->minimum_image_big(xtmp,ytmp,ztmp);
   return {xtmp, ytmp, ztmp};
 }
 

Original file line number	Diff line number	Diff line change
`@@ -205,7 +205,7 @@ cvm::rvector colvarproxy_lammps::position_distance(cvm::atom_pos const &pos1,`
`205`	`205`	`double xtmp = pos2.x - pos1.x;`
`206`	`206`	`double ytmp = pos2.y - pos1.y;`
`207`	`207`	`double ztmp = pos2.z - pos1.z;`
`208`		`- _lmp->domain->minimum_image(xtmp,ytmp,ztmp);`
	`208`	`+ _lmp->domain->minimum_image_big(xtmp,ytmp,ztmp);`
`209`	`209`	`return {xtmp, ytmp, ztmp};`
`210`	`210`	`}`
`211`	`211`