stattest2.cpp

随机数 基于c++的Monte Carlo源代码

#include <iostream>

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

int main(void)
{
  mcsLibreCollectStats stuff;

  /*  
      The following data set comes from Octave's "normal_inv" command.
      There are 100 normally 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( 103.805978401296);
  stuff.add_value( 104.175277041221);
  stuff.add_value( 102.866054018321);
  stuff.add_value( 102.836651655618);
  stuff.add_value( 103.492279072354);
  stuff.add_value( 98.9036767491981);
  stuff.add_value( 100.297308478322);
  stuff.add_value( 102.341139611905);
  stuff.add_value( 99.1874498151906);
  stuff.add_value( 97.2033258189596);
  stuff.add_value( 101.626272564131);
  stuff.add_value( 98.5289766893555);
  stuff.add_value( 98.4778547064023);
  stuff.add_value( 102.234273368514);
  stuff.add_value( 101.980661966706);
  stuff.add_value( 98.1413193957097);
  stuff.add_value( 94.8490260419169);
  stuff.add_value( 104.42800859781);
  stuff.add_value( 98.7043280745783);
  stuff.add_value( 100.07812910571);
  stuff.add_value( 104.382970137208);
  stuff.add_value( 99.5659411337695);
  stuff.add_value( 101.448078610033);
  stuff.add_value( 100.174050453663);
  stuff.add_value( 98.5453905549026);
  stuff.add_value( 98.9973532155473);
  stuff.add_value( 100.988751823421);
  stuff.add_value( 98.4543456430133);
  stuff.add_value( 97.3791022591256);
  stuff.add_value( 102.053115300667);
  stuff.add_value( 95.8175182463874);
  stuff.add_value( 101.314426394675);
  stuff.add_value( 98.9704398945862);
  stuff.add_value( 101.41777048122);
  stuff.add_value( 100.155053192978);
  stuff.add_value( 97.4680239258081);
  stuff.add_value( 101.063252914462);
  stuff.add_value( 100.174801145391);
  stuff.add_value( 100.444140752729);
  stuff.add_value( 100.056793770736);
  stuff.add_value( 98.3216608033285);
  stuff.add_value( 99.1760897281355);
  stuff.add_value( 102.197523203382);
  stuff.add_value( 92.8616474281238);
  stuff.add_value( 98.2844982168141);
  stuff.add_value( 100.503459694797);
  stuff.add_value( 95.4041037987126);
  stuff.add_value( 104.189825070628);
  stuff.add_value( 96.2999839132369);
  stuff.add_value( 101.415622808945);
  stuff.add_value( 102.10130897974);
  stuff.add_value( 100.925268365182);
  stuff.add_value( 96.7398855894399);
  stuff.add_value( 100.725685579466);
  stuff.add_value( 108.767630461747);
  stuff.add_value( 98.9415955161789);
  stuff.add_value( 100.951066499067);
  stuff.add_value( 101.534965702031);
  stuff.add_value( 103.417168848783);
  stuff.add_value( 100.494754845197);
  stuff.add_value( 101.075059382415);
  stuff.add_value( 94.9145027842828);
  stuff.add_value( 103.431248716257);
  stuff.add_value( 102.455390928742);
  stuff.add_value( 101.195488581285);
  stuff.add_value( 102.778572865596);
  stuff.add_value( 100.926045394751);
  stuff.add_value( 99.5040871889738);
  stuff.add_value( 101.062023069284);
  stuff.add_value( 101.134131775524);
  stuff.add_value( 102.65536108468);
  stuff.add_value( 100.906351009869);
  stuff.add_value( 102.768735040896);
  stuff.add_value( 99.8497301261721);
  stuff.add_value( 101.266044118979);
  stuff.add_value( 92.4198407730632);
  stuff.add_value( 103.31431198593);
  stuff.add_value( 106.947654583544);
  stuff.add_value( 97.8374801993542);
  stuff.add_value( 96.6945682084629);
  stuff.add_value( 98.3069048200199);
  stuff.add_value( 98.2245594140177);
  stuff.add_value( 99.0084821567436);
  stuff.add_value( 99.7785792250738);
  stuff.add_value( 98.1577315551219);
  stuff.add_value( 99.6862044912687);
  stuff.add_value( 103.039350987945);
  stuff.add_value( 102.970156869157);
  stuff.add_value( 103.534475596802);
  stuff.add_value( 105.323191615961);
  stuff.add_value( 101.783774724164);
  stuff.add_value( 94.309217866653);
  stuff.add_value( 98.6975526469196);
  stuff.add_value( 102.559533103906);
  stuff.add_value( 98.3945273606385);
  stuff.add_value( 94.385165121521);
  stuff.add_value( 102.006254530308);
  stuff.add_value( 95.6332068941343);
  stuff.add_value( 97.7158825448676);
  stuff.add_value( 101.412879894476);

  cout << endl;
  cout << "Octave mean:  100.24" << endl;
  cout << "MCS-Libre mean:  " << stuff.get_mean_value() << endl;
  cout << "Octave Stdev:  2.8732" << 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:  10024" << endl;
  cout << "MCS-Libre Sum:  " << stuff.get_sum_of_values() << endl;
  cout << "Octave Max:  108.77" << endl;
  cout << "MCS-Libre Max:  " << stuff.get_max_value() << endl;
  cout << "Octave Min:  92.420" << endl;
  cout << "MCS-Libre Min:  " << stuff.get_min_value() << endl;
  cout << endl;
  return 0;
}