--- branches/mmeineke/OOPSE/libmdtools/StatWriter.cpp 2003/03/21 17:42:12 377 +++ trunk/OOPSE/libmdtools/StatWriter.cpp 2003/10/28 16:03:37 829 @@ -1,4 +1,6 @@ -#include +#define _FILE_OFFSET_BITS 64 + +#include #include #include @@ -16,7 +18,7 @@ StatWriter::StatWriter( SimInfo* the_entry_plug ){ strcpy( outName, entry_plug->statusName ); - std::cerr << "Opening " << outName << " for stat\n"; + //std::cerr << "Opening " << outName << " for stat\n"; outFile.open(outName, ios::out | ios::trunc ); @@ -30,7 +32,7 @@ StatWriter::StatWriter( SimInfo* the_entry_plug ){ } //outFile.setf( ios::scientific ); - outFile << "#time(fs)\ttot_E\tpotential\tkinetic\ttemperature\n"; + outFile << "#time(fs)\ttot_E\tpotential\tkinetic\ttemperature\tpressure\tvolume\tconserved quantity\n"; @@ -61,13 +63,16 @@ void StatWriter::writeStat( double currentTime ){ void StatWriter::writeStat( double currentTime ){ - double totE, potE, kinE, temp; + double totE, potE, kinE, temp, press, vol; + double conservedQuantity; totE = tStats->getTotalE(); potE = tStats->getPotential(); kinE = tStats->getKinetic(); temp = tStats->getTemperature(); - + press = tStats->getPressure(); + vol = tStats->getVolume(); + conservedQuantity = entry_plug->the_integrator->getConservedQuantity(); #ifdef IS_MPI if(worldRank == 0 ){ #endif // is_mpi @@ -78,7 +83,11 @@ void StatWriter::writeStat( double currentTime ){ << totE << "\t" << potE << "\t" << kinE << "\t" - << temp << "\n"; + << temp << "\t" + << press << "\t" + << vol << "\t" + << conservedQuantity << "\n"; + outFile.flush(); #ifdef IS_MPI