QMCStatistic.h

Go to the documentation of this file.

//            QMcBeaver
//
//         Constructed by 
//
//     Michael Todd Feldmann 
//              and 
//   David Randall "Chip" Kent IV
//
// Copyright 2000.  All rights reserved.
//
// drkent@users.sourceforge.net mtfeldmann@users.sourceforge.net

#ifndef QMCSTATISTIC_H
#define QMCSTATISTIC_H

#include <iostream>
#include <string>
#include <math.h>

#ifdef PARALLEL
#include <mpi.h>
#endif

using namespace std;


class QMCStatistic
{
private:
  long double sum;
  long double sum2;
  long double sum3;
  long double sum4;
  long double weights;
  unsigned long nsamples;

public:
  QMCStatistic();


  void zeroOut();


  unsigned long getNumberSamples() const;


  double getAverage() const;


  double getVariance() const;


  double getStandardDeviation() const;

  double getSkewness() const;
  double getKurtosis() const;

  void newSample(long double s, long double weight);
  void quickSample(long double s, long double weight);

  void operator = ( const QMCStatistic &rhs);

  QMCStatistic operator + (const QMCStatistic &rhs);

  void reWeight(double w);

  void toXML(ostream& strm);


  bool readXML(istream& strm);


  friend ostream& operator << (ostream& strm, const QMCStatistic & rhs);


#ifdef PARALLEL

private:
  static bool mpiTypeCreated;

  
  static void buildMpiType();

  static void buildMpiReduce();


  static void Reduce_Function(QMCStatistic *in, QMCStatistic *inout, 
                              int *len, MPI_Datatype *dptr);

public:

  static MPI_Datatype MPI_TYPE;


  static MPI_Op MPI_REDUCE;

#endif

};

#endif

---