probeSG.cpp

Go to the documentation of this file.

#include <fstream>
#include "probeSG.h"

ProbeSG::ProbeSG( string _theName) :
  Sum_S(0), Sum_G(0), StartTime(0), lastTime(0), Probe(_theName, false)
{
}

void ProbeSG::Reset( void )
{
  Sum_S = 0;
  Sum_G = 0;
  Sum_S_part = 0;
  Sum_G_part = 0;
  StartTime = 0;
  lastTime = 0;
  Probe::Reset();
}

void ProbeSG::Zap( double time )
{
  Sum_S = 0;
  Sum_G = 0;
  Sum_S_part = 0;
  Sum_G_part = 0;
  StartTime = time;
  lastTime = time;
}

void ProbeSG::Observe_S( double value )
{
  Sum_S += value;
  Sum_S_part += value;
  Sum_G += value;
  Sum_G_part += value;
}

void ProbeSG::Observe_G( double value )
{
  Sum_G += value;
  Sum_G_part += value;
}

void ProbeSG::WriteToFile( double time )
{
  // output files are not left open, as they are when needed.
  // Append is a Probe's internal variable, false when the probe is created.
  // You can use it to choose if your file is to create from scratch or just appended.
  // Note that NePSing manages the probe renaming from one run to another.
  // theRealName.c_str() is the actual name of the output file, after ProbeManager unmangling.
  ofstream ofParam(theRealName.c_str(), Append ? ios::app : ios::out);
  
  Append = true;
  
  if (!ofParam) {
    throw ProbeException("Error opening file: " + theRealName);
  }
  
  ofParam << time << "\t";
  if( time == StartTime ) {
    // beware of the zero division ! 0/(0-0) is a bad idea even if IEEE numbers are active.
    ofParam << "0" << endl;
  }
  else {
    // feel free to use your own output format, but a gnuplot-like one is always good :-P
    ofParam << Sum_S/(time-StartTime) << "\t" << Sum_G/(time-StartTime);
    ofParam << "\t" << Sum_S_part/(time-lastTime) << "\t" << Sum_G_part/(time-lastTime) << endl;
  }
  
  Sum_S_part = 0;
  Sum_G_part = 0;
  lastTime = time;
}

---