Changeset 1f2e46 for molecuilder/src/atom_trajectoryparticle.cpp
- Timestamp:
- Apr 22, 2010, 2:00:03 PM (16 years ago)
- Children:
- 423c7b
- Parents:
- c43766
- File:
-
- 1 edited
-
molecuilder/src/atom_trajectoryparticle.cpp (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
molecuilder/src/atom_trajectoryparticle.cpp
rc43766 r1f2e46 198 198 // throw a dice to determine whether it gets hit by a heat bath particle 199 199 if (((((rand()/(double)RAND_MAX))*configuration->TempFrequency) < 1.)) { 200 Log() << Verbose(3) << "Particle " << *this << " was hit (sigma " << sigma << "): " << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << " -> ";200 DoLog(3) && (Log() << Verbose(3) << "Particle " << *this << " was hit (sigma " << sigma << "): " << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << " -> "); 201 201 // pick three random numbers from a Boltzmann distribution around the desired temperature T for each momenta axis 202 202 for (int d=0; d<NDIM; d++) { 203 203 U[d] = gsl_ran_gaussian (r, sigma); 204 204 } 205 Log() << Verbose(2) << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl;205 DoLog(2) && (Log() << Verbose(2) << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl); 206 206 } 207 207 for (int d=0; d<NDIM; d++)
Note:
See TracChangeset
for help on using the changeset viewer.
