[ViewVC] Diff of: group/trunk/OOPSE/libmdtools/NVT.cpp

Comparing trunk/OOPSE/libmdtools/NVT.cpp (file contents):
Revision 1097 by gezelter, Mon Apr 12 20:32:20 2004 UTC vs.
Revision 1129 by tim, Thu Apr 22 03:29:30 2004 UTC

#	Line 50 \| Line 50 \| *template<typename T> NVT<T>::NVT ( SimInfo theInfo, F**
50		}
51		}
52
53	–
54	–	std::cerr << "building oldVel with \t" << integrableObjects.size() << "\n";
53		oldVel = new double[3*integrableObjects.size()];
54		oldJi = new double[3*integrableObjects.size()];
55		}
#	Line 120 \| Line 118 \| template<typename T> void NVT<T>::moveA() {
118		// Finally, evolve chi a half step (just like a velocity) using
119		// temperature at time t, not time t+dt/2
120
121	+	//std::cerr << "targetTemp = " << targetTemp << " instTemp = " << instTemp << " tauThermostat = " << tauThermostat << " integral of Chi = " << integralOfChidt << "\n";
122	+
123		chi += dt2 * ( instTemp / targetTemp - 1.0) / (tauThermostat*tauThermostat);
124		integralOfChidt += chi*dt2;
125
#	Line 259 \| Line 259 \| template<typename T> double NVT<T>::getConservedQuanti
259		double thermostat_kinetic;
260		double thermostat_potential;
261
262	<	fkBT = (double)(info->getNDF() ) * kB * targetTemp;
262	>	fkBT = (double)(info->ndf) * kB * targetTemp;
263
264		Energy = tStats->getTotalE();
265

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines