Changeset 543ce4 for molecuilder/src/atom_trajectoryparticle.cpp

Timestamp:

Nov 4, 2009, 7:56:04 PM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

4ef101, aa8542

Parents:

ec70ec

Message:

Huge change from ofstream * (const) out --> Log().

• first shift was done via regular expressions
• then via error messages from the code
• note that class atom, class element and class molecule kept in parts their output stream, was they print to file.
• make check runs fine
• MISSING: Verbosity is not fixed for everything (i.e. if no endl; is present and next has Verbose(0) ...)

Signed-off-by: Frederik Heber <heber@…>

File:

• molecuilder/src/atom_trajectoryparticle.cpp (modified) (4 diffs)

molecuilder/src/atom_trajectoryparticle.cpp

@@ -10 +10 @@
 #include "config.hpp"
 #include "element.hpp"
+#include "log.hpp"
 #include "parser.hpp"
 #include "verbose.hpp"
@@ -71 +72 @@
 {
   if (Trajectory.R.size() <= (unsigned int)(MaxSteps)) {
-    //cout << "Increasing size for trajectory array of " << keyword << " to " << (MaxSteps+1) << "." << endl;
+    //Log() << Verbose(0) << "Increasing size for trajectory array of " << keyword << " to " << (MaxSteps+1) << "." << endl;
     Trajectory.R.resize(MaxSteps+1);
     Trajectory.U.resize(MaxSteps+1);
@@ -120 +121 @@
 //      out << ")\t(";
 //      for (int d=0;d<NDIM;d++)
-//        cout << Trajectory.U.at(NextStep).x[d] << " ";          // next step
+//        Log() << Verbose(0) << Trajectory.U.at(NextStep).x[d] << " ";          // next step
 //      out << ")" << endl;
 };
@@ -197 +198 @@
     // throw a dice to determine whether it gets hit by a heat bath particle
     if (((((rand()/(double)RAND_MAX))*configuration->TempFrequency) < 1.)) {
-      cout << Verbose(3) << "Particle " << *this << " was hit (sigma " << sigma << "): " << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << " -> ";
+      Log() << Verbose(3) << "Particle " << *this << " was hit (sigma " << sigma << "): " << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << " -> ";
       // pick three random numbers from a Boltzmann distribution around the desired temperature T for each momenta axis
       for (int d=0; d<NDIM; d++) {
         U[d] = gsl_ran_gaussian (r, sigma);
       }
-      cout << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl;
+      Log() << Verbose(0) << sqrt(U[0]*U[0]+U[1]*U[1]+U[2]*U[2]) << endl;
     }
     for (int d=0; d<NDIM; d++)