tryrand.cpp

Go to the documentation of this file.

#define WANT_STREAM
#include "StdAfx.h"

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

#include "include.h"
#include "newran.h"

#ifdef use_namespace
using namespace NEWRAN;
#endif

void test1(int);
void test2(int);
void test3(int);
void test4(int);
void Histogram(Random*, int);          // draw histogram with n obsv
void Hist(Random* rx, int n , const int cnt );          // draw histogram with n obsv

int main( int argc , char **argv )
{

        Random::Set(0.46875);
        long n = 1000;
    Uniform u , u1 , u2 , u3 , u4 , u5;

        SumRandom       *sr;

        sr = new SumRandom( 2 * u * 2 * u );
        Hist( sr , 1000 , 42 );
        delete sr;

   {

      SumRandom sr1 = (u * u) + (1 - (u * u));
      Histogram(&sr1,n);
      cout << "Mean and variance should be 0.25 and 0.048611\n";
      
      SumRandom sr2 = u + 2 * u * u;
      Histogram(&sr2,n);
      cout << "Mean and variance should be 0.25 and 0.048611\n";
      
      SumRandom sr3 = (u-0.5) * (u-0.5) + u;
          Histogram(&sr3,n);
      cout << "Mean and variance should be 0.0 and 0.006944\n";

   }

   Real* s1; Real* s2; Real* s3; Real* s4;
   cout << "\nBegin test\n";   // Forces cout to allocate memory at beginning
   { s1 = new Real[8000]; delete [] s1; }
   { s3 = new Real; delete s3;}
   {

      Real* A = new Real[3750];

      long n_large = 1000000;

/*
      Uniform u;
      SumRandom sr1 = u*u;
      SumRandom sr2 = (u-0.5)*(u-0.5);
      Histogram(&sr1,n);
      cout << "Mean and variance should be 0.25 and 0.048611" << endl;
      Histogram(&sr2,n);
      cout << "Mean and variance should be 0.0 and 0.006944" << endl;
*/   
                //      RING    
#if 1
/*
                {
      static Real probs[]={.05 , 0.25, .025 , 0 , 0.025 ,.3 , .025 , 0 , 0.025 , .25 , .05};
      DiscreteGen discrete(11,probs); Normal nn;
      SumRandom sr = discrete + (nn * 0.25)(2) + 10.0;
      Histogram(&discrete,n);
      Histogram(&sr,n);
   }
*/
                {
      static Real probs[]={.05 , 0.25, .02 , 0.01 , 0.02 ,.3 , .02 , 0.01 , 0.02 , .25 , .05};
      DiscreteGen discrete(11,probs); Normal nn;
      SumRandom sr = discrete + (nn * 0.25)(2) + 10.0;
      Hist(&discrete,n,50);
      Hist(&sr,n,50);
                }

                {
                        
                static Real probs[]={ 0.4 , 0.000 ,0.000 ,0.12 
                        ,0.000 ,0.000 ,0.000
                        ,0.12 
                        ,0.000 ,0.000 ,0.000,0.000 ,0.000 ,0.000,0.000 ,0.000
                        ,0.12 , 0.12 , 0.12 };
      DiscreteGen discrete(19,probs); 
      Hist(&discrete,n,19);
                Uniform nn;     //      Normal nn;
      {
                        SumRandom sr = discrete + (nn * 0.1)( 2 );
//                      Hist(&sr,n,7);
                }
      {
                        SumRandom sr = discrete + (nn * 0.25)( 2 );
                        Hist(&sr,n,19);
                }
      {
                        SumRandom sr = discrete + (nn * 0.5)( 2 );
//                      Hist(&sr,n,7);
                }
      {
                        SumRandom sr = discrete + (nn * 0.75)( 2 );
//                      Hist(&sr,n,7);
                }
   }
#endif

//      test1(n);
                test2(n);
//      test3(n_large);
//      test4(n);

      cout << "\nEnd of tests\n";

      delete [] A;
   }

 /*
 { s2 = new Real[8000]; delete [] s2; }
   cout << "\n(The following memory checks are probably not valid with all\n";
   cout << "compilers - see documentation)\n";
   cout << "\nChecking for lost memory: "
      << (unsigned long)s1 << " " << (unsigned long)s2 << " ";
   if (s1 != s2) cout << " - error\n"; else cout << " - ok\n\n";
   { s4 = new Real; delete s4;}
   cout << "\nChecking for lost memory: "
      << (unsigned long)s3 << " " << (unsigned long)s4 << " ";
   if (s3 != s4) cout << " - error\n"; else cout << " - ok\n\n";
*/
   return 0;
}

---