Debug

abkein · Jan 6, 2025 · 8dc560f · 8dc560f
1 parent 46d82f0
commit 8dc560f
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/fix_capture_vel.cpp b/src/fix_capture_vel.cpp
@@ -210,9 +210,9 @@ void FixCaptureVel::pre_force(int vflag)
       if (vm > sigmas[atom->type[i]]) {
         ++ncaptured[0];
         const double sigma = ::sqrt(sigmas[atom->type[i]] / nsigmasq);
-        v[i][0] *= v[i][0] / sigma;
-        v[i][1] = vrandom->gaussian() * sigma;
-        v[i][2] = vrandom->gaussian() * sigma;
+        v[i][0] *= sigma / v[i][0];
+        v[i][1] *= sigma / v[i][1];
+        v[i][2] *= sigma / v[i][2];
       }
       if (delete_overlap) {
         const double xtmp = x[i][0];