Completing the correction in fix langevin/spin.

2022-09-26 07:44:32 +02:00 · 2022-09-26 07:44:32 +02:00 · ac6cb062e6
parent 90f26f05fc
commit ac6cb062e6
1 changed files with 0 additions and 1 deletions
--- a/src/SPIN/fix_langevin_spin.cpp
+++ b/src/SPIN/fix_langevin_spin.cpp
@ -109,7 +109,6 @@ void FixLangevinSpin::init()
  double hbar = force->hplanck/MY_2PI;  // eV/(rad.THz)
  double kb = force->boltz;             // eV/K

-  //D = (alpha_t*gil_factor*kb*temp);
  D = (alpha_t*(1.0+(alpha_t)*(alpha_t))*kb*temp);
  D /= (hbar*dts);
  sigma = sqrt(2.0*D);