particlesampler.h

Go to the documentation of this file.

#ifndef PARTICLESAMPLER_H
#define PARTICLESAMPLER_H


#include <cassert>
#include <string>
#include <iostream>
#include <fstream>
using namespace std;

#include <histogram.h>
#include <particle.h>
#include <typedefs.h>


/*
\brief Periodically sample points in the particle simulation.

Use:
  Configure the state to sample different data.
   eg 
   0 No sampling
   1 Samples the magnitude of overall veclocity
   2 Samples the speed in the x-axis
   3 Samples the speed in the y-axis

The samplecountermax is the number of times this object
 is iterated over before it samples.
 The filecountermax is then the number of times that
 sampling occures before writing to the file.
*/
class particlesampler
{
  uintc & state;

  uint samplecounter;
  uint filecounter;
  histogram<double> hist;
  string filename; 
         
  particle const * & pi;
  uintc & numParticles;

  void sample();
  void writefile();
public:

  uint samplecountermax;   

  uint filecountermax;

  particlesampler
  (
    uintc & state_,
    uintc samplecountermax_,
    uintc filecountermax_,
    histogram<double> const hist_,
    string const filename_,
    particle const * & pi_,
    uintc & numParticles_
  );

  void set(histogram<double> const hist_)
    { hist=hist_; }
 
  void reset();

  void operator ++ ();
//    { ++samplecounter; if ((samplecounter%samplecountermax)==0) sample(); }

};



#endif