Home Research Papers Posters Manuals Talks Download People Help Internal

util/ckstatistics.h

Go to the documentation of this file.

/*
(Copied out of) Orion's Standard Library
written by 
Orion Sky Lawlor, olawlor@acm.org, 11/25/2002
*/
#ifndef __UIUC_OSL_STATISTICS_H
#define __UIUC_OSL_STATISTICS_H

#include "math.h" //for sqrt
#include "stdio.h" //for fprintf
#include "pup.h" //for pup

template<class real,class ret>
class CkSampleT {
    real lo,hi; //smallest and largest values
    ret sum; //sum of values
    ret sq; //sum of squares of values
    int n; //number of values
public:
    CkSampleT(void) {
        lo=(real)1.0e20; hi=(real)-1.0e20;
        sum=sq=(ret)0;
        n=0;
    }
    inline void add(real r) {
        if (r<lo) lo=r;
        if (r>hi) hi=r;
        sum+=(ret)r;
        sq+=(ret)(r*r);
        n++;
    }
    inline void operator+=(real r) { add(r); }

    inline void add(const CkSampleT<real,ret> &errs) {
        if (errs.lo<lo) lo=errs.lo;
        if (errs.hi>hi) hi=errs.hi;
        sum+=errs.sum;
        sq+=errs.sq;
        n+=errs.n;
    }
    
    ret getMean(void) const {
        return sum/n;
    }
    ret getVariance(void) const {
        return (sq-sum*sum/n)/(n-1);
    }
    ret getStddev(void) const {
        return (ret)sqrt(getVariance());
    }
    real getMin(void) const {return lo;}
    real getMax(void) const {return hi;}
    int getCount(void) const {return n;}

    ret getRMS(void) const {return sqrt(sq/n);}
    
    void print(FILE *dest) {
        fprintf(dest,"ave= %g  stddev= %g  min= %g  max= %g  n= %d\n",
            (double)getMean(), (double)getStddev(), (double)getMin(), (double)getMax(), (int)getCount());
    }

        void printMinAveMax(FILE *dest) {
        fprintf(dest,"min= %g  ave= %g  max= %g \n",
            (double)getMin(), (double)getMean(), (double)getMax());
    }

    void print(void) {print(stdout);}
    
    void pup(PUP::er &p) {
        p|lo; p|hi;
        p|sum; p|sq;
        p|n;
    }
    friend inline void operator|(PUP::er &p,CkSampleT<real,ret> &s) {s.pup(p);}
};
typedef CkSampleT<double,double> CkSample;

#endif

---