From 0ae40c596799aa2c69c12e93a606d07aecc99572 Mon Sep 17 00:00:00 2001 From: Viktor Klippenstein Date: Wed, 22 Jan 2020 17:13:10 +0100 Subject: [PATCH] Add unit conversion in setup() --- src/fix_langevin.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/fix_langevin.cpp b/src/fix_langevin.cpp index f8f85fcf30..4ec34a3631 100644 --- a/src/fix_langevin.cpp +++ b/src/fix_langevin.cpp @@ -339,7 +339,7 @@ void FixLangevin::setup(int vflag) if (rmass) { for (int i = 0; i < nlocal; i++) if (mask[i] & groupbit) { - dtfm = 0.5 * dt / rmass[i]; + dtfm = force->ftm2v * 0.5 * dt / rmass[i]; v[i][0] -= dtfm * f[i][0]; v[i][1] -= dtfm * f[i][1]; v[i][2] -= dtfm * f[i][2]; @@ -355,7 +355,7 @@ void FixLangevin::setup(int vflag) } else { for (int i = 0; i < nlocal; i++) if (mask[i] & groupbit) { - dtfm = 0.5 * dt / mass[type[i]]; + dtfm = force->ftm2v * 0.5 * dt / mass[type[i]]; v[i][0] -= dtfm * f[i][0]; v[i][1] -= dtfm * f[i][1]; v[i][2] -= dtfm * f[i][2]; @@ -389,7 +389,7 @@ void FixLangevin::setup(int vflag) if (rmass) { for (int i = 0; i < nlocal; i++) if (mask[i] & groupbit) { - dtfm = 0.5 * dt / rmass[i]; + dtfm = force->ftm2v * 0.5 * dt / rmass[i]; v[i][0] += dtfm * f[i][0]; v[i][1] += dtfm * f[i][1]; v[i][2] += dtfm * f[i][2]; @@ -401,7 +401,7 @@ void FixLangevin::setup(int vflag) } else { for (int i = 0; i < nlocal; i++) if (mask[i] & groupbit) { - dtfm = 0.5 * dt / mass[type[i]]; + dtfm = force->ftm2v * 0.5 * dt / mass[type[i]]; v[i][0] += dtfm * f[i][0]; v[i][1] += dtfm * f[i][1]; v[i][2] += dtfm * f[i][2];