stattest1.cpp

Free C++ toolkit to facilitate Monte-Carlo simulation. This is a library covered under the LGPL. "MC

#include <iostream>

#include "../include/mcsLibreCollectStats.h"

int main(void)
{
  mcsLibreCollectStats stuff;

  /*  
      The following data set comes from Octave's "rand" command.
      There are 100 uniformly distributed values here. 

      Here the goal is to compare Octave's output to MCS-Libre's
      output to determine if the "mcsLibreCollectStats" class is
      working correctly.

  */

  stuff.add_value( 0.773395299911499);
  stuff.add_value( 0.654982209205627);
  stuff.add_value( 0.528605937957764);
  stuff.add_value( 0.383064240217209);
  stuff.add_value( 0.698346138000488);
  stuff.add_value( 0.703774750232697);
  stuff.add_value( 0.920608162879944);
  stuff.add_value( 0.14160367846489);
  stuff.add_value( 0.647386312484741);
  stuff.add_value( 0.532176375389099);
  stuff.add_value( 0.115888319909573);
  stuff.add_value( 0.799268782138824);
  stuff.add_value( 0.421339303255081);
  stuff.add_value( 0.888470828533173);
  stuff.add_value( 0.626602292060852);
  stuff.add_value( 0.160668984055519);
  stuff.add_value( 0.134032368659973);
  stuff.add_value( 0.946596384048462);
  stuff.add_value( 0.59737890958786);
  stuff.add_value( 0.0461915396153927);
  stuff.add_value( 0.371142834424973);
  stuff.add_value( 0.880790531635284);
  stuff.add_value( 0.161894351243973);
  stuff.add_value( 0.258866459131241);
  stuff.add_value( 0.481668800115585);
  stuff.add_value( 0.298350423574448);
  stuff.add_value( 0.762927234172821);
  stuff.add_value( 0.62574166059494);
  stuff.add_value( 0.0963747426867485);
  stuff.add_value( 0.0324103720486164);
  stuff.add_value( 0.542961061000824);
  stuff.add_value( 0.820049464702606);
  stuff.add_value( 0.63988870382309);
  stuff.add_value( 0.663719415664673);
  stuff.add_value( 0.307454794645309);
  stuff.add_value( 0.75092214345932);
  stuff.add_value( 0.486743360757828);
  stuff.add_value( 0.321563184261322);
  stuff.add_value( 0.677633285522461);
  stuff.add_value( 0.962463319301605);
  stuff.add_value( 0.167001754045486);
  stuff.add_value( 0.224868014454842);
  stuff.add_value( 0.276210963726044);
  stuff.add_value( 0.322929173707962);
  stuff.add_value( 0.538146734237671);
  stuff.add_value( 0.418469220399857);
  stuff.add_value( 0.273543417453766);
  stuff.add_value( 0.555981934070587);
  stuff.add_value( 0.179065570235252);
  stuff.add_value( 0.395357549190521);
  stuff.add_value( 0.209782853722572);
  stuff.add_value( 0.83413153886795);
  stuff.add_value( 0.407693028450012);
  stuff.add_value( 0.410890281200409);
  stuff.add_value( 0.71085649728775);
  stuff.add_value( 0.0032923580147326);
  stuff.add_value( 0.824960172176361);
  stuff.add_value( 0.793672323226929);
  stuff.add_value( 0.621214807033539);
  stuff.add_value( 0.896124541759491);
  stuff.add_value( 0.143671855330467);
  stuff.add_value( 0.682688474655151);
  stuff.add_value( 0.337428867816925);
  stuff.add_value( 0.998228967189789);
  stuff.add_value( 0.717763364315033);
  stuff.add_value( 0.600896418094635);
  stuff.add_value( 0.457094132900238);
  stuff.add_value( 0.43616458773613);
  stuff.add_value( 0.348047077655792);
  stuff.add_value( 0.812190592288971);
  stuff.add_value( 0.19325764477253);
  stuff.add_value( 0.773584723472595);
  stuff.add_value( 0.813496530056);
  stuff.add_value( 0.664782106876373);
  stuff.add_value( 0.96471244096756);
  stuff.add_value( 0.35065758228302);
  stuff.add_value( 0.952796697616577);
  stuff.add_value( 0.200670138001442);
  stuff.add_value( 0.235314786434174);
  stuff.add_value( 0.213651105761528);
  stuff.add_value( 0.726191401481628);
  stuff.add_value( 0.855210065841675);
  stuff.add_value( 0.754241824150085);
  stuff.add_value( 0.197466731071472);
  stuff.add_value( 0.793697774410248);
  stuff.add_value( 0.887583196163177);
  stuff.add_value( 0.482992380857468);
  stuff.add_value( 0.521092057228088);
  stuff.add_value( 0.779822528362274);
  stuff.add_value( 0.12359619140625);
  stuff.add_value( 0.0952867269515991);
  stuff.add_value( 0.583913087844849);
  stuff.add_value( 0.754240393638611);
  stuff.add_value( 0.856806576251984);
  stuff.add_value( 0.0383913069963455);
  stuff.add_value( 0.107383340597153);
  stuff.add_value( 0.55610066652298);
  stuff.add_value( 0.0410389825701714);
  stuff.add_value( 0.802560925483704);
  stuff.add_value( 0.488462775945663);

  cout << endl;
  cout << "Octave mean:  0.51271" << endl;
  cout << "MCS-Libre mean:  " << stuff.get_mean_value() << endl;
  cout << "Octave Stdev:  0.27863" << endl; 
  cout << "MCS-Libre Stdev:  " << stuff.get_stdev_value() << endl;
  cout << "Octave Sample Number:  100" << endl;
  cout << "MCS-Libre Sample NUmber:  " << stuff.get_number_samples() << endl;
  cout << "Octave Sum:  51.271" << endl;
  cout << "MCS-Libre Sum:  " << stuff.get_sum_of_values() << endl;
  cout << "Octave Max:  0.99823" << endl;
  cout << "MCS-Libre Max:  " << stuff.get_max_value() << endl;
  cout << "Octave Min:  0.0032924" << endl;
  cout << "MCS-Libre Min:  " << stuff.get_min_value() << endl;
  cout << endl;
  return 0;
}