From 99f42744833cbd7f9d97ae5657c3df2553150d5c Mon Sep 17 00:00:00 2001
From: julient31 <julien.tranchida1@gmail.com>
Date: Thu, 24 May 2018 06:59:40 -0600
Subject: [PATCH] Commit JT 052418 rm comments
---
src/SPIN/compute_spin.cpp | 1 -
src/SPIN/fix_langevin_spin.cpp | 7 -------
2 files changed, 8 deletions(-)
diff --git a/src/SPIN/compute_spin.cpp b/src/SPIN/compute_spin.cpp
index 47d88b6c64..ab3bab2487 100644
--- a/src/SPIN/compute_spin.cpp
+++ b/src/SPIN/compute_spin.cpp
@@ -132,7 +132,6 @@ void ComputeSpin::compute_vector()
magtot[2] *= scale;
magtot[3] = sqrt((magtot[0]*magtot[0])+(magtot[1]*magtot[1])+(magtot[2]*magtot[2]));
spintemperature = hbar*tempnumtot;
- //spintemperature /= (2.0*kb*tempdenomtot);
spintemperature /= (kb*tempdenomtot);
vector[0] = magtot[0];
diff --git a/src/SPIN/fix_langevin_spin.cpp b/src/SPIN/fix_langevin_spin.cpp
index b4e7f67653..d12ec11237 100644
--- a/src/SPIN/fix_langevin_spin.cpp
+++ b/src/SPIN/fix_langevin_spin.cpp
@@ -134,10 +134,8 @@ void FixLangevinSpin::init()
double hbar = force->hplanck/MY_2PI; // eV/(rad.THz)
double kb = force->boltz; // eV/K
D = (MY_2PI*alpha_t*gil_factor*kb*temp);
- //D = (alpha_t*gil_factor*kb*temp);
D /= (hbar*dts);
sigma = sqrt(2.0*D);
- //sigma = sqrt(D);
}
/* ---------------------------------------------------------------------- */
@@ -176,11 +174,6 @@ void FixLangevinSpin::add_temperature(double fmi[3])
double rx = sigma*(2.0*random->uniform() - 1.0);
double ry = sigma*(2.0*random->uniform() - 1.0);
double rz = sigma*(2.0*random->uniform() - 1.0);
- //printf("test rd : %g \n",2.0*random->uniform() - 1.0);
- //printf("test gaussian : %g \n", random->gaussian());
- //double rx = sigma*(random->gaussian());
- //double ry = sigma*(random->gaussian());
- //double rz = sigma*(random->gaussian());
// adding the random field
--
GitLab