ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/NVT.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/NVT.cpp (file contents):
Revision 1125 by gezelter, Mon Apr 19 22:13:01 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";
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;
# Line 261 | 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  

Diff Legend

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