testEngineCopy.cc File Reference

#include "CLHEP/Units/GlobalPhysicalConstants.h"
#include "CLHEP/Random/Randomize.h"
#include "CLHEP/Random/NonRandomEngine.h"
#include "CLHEP/Random/defs.h"
#include <iostream>
#include <iomanip>
#include <vector>

Go to the source code of this file.

Macros
#define	CLEAN_OUTPUT
#define	TEST_ENGINE_COPY
#define	VERBOSER
#define	VERBOSER2

Functions
std::ofstream	output ("testEngineCopy.cout")
bool	equals01 (const std::vector< double > &ab)
bool	equals (double a, double b)
std::vector< double >	aSequence (int n)
template<class E >
int	vectorTest64 (int n)
template<>
int	vectorTest64< NonRandomEngine > (int n)
template<class E >
E	vectorRestore1 (int n, std::vector< double > &v)
template<>
NonRandomEngine	vectorRestore1< NonRandomEngine > (int n, std::vector< double > &v)
template<class E >
int	vectorRestore2 (E &f, const std::vector< double > &v)
template<class E >
int	vectorRestore (int n)
int	main ()

Macro Definition Documentation

CLEAN_OUTPUT

#define CLEAN_OUTPUT

Definition at line 12 of file testEngineCopy.cc.

TEST_ENGINE_COPY

#define TEST_ENGINE_COPY

Definition at line 24 of file testEngineCopy.cc.

VERBOSER

#define VERBOSER

Definition at line 26 of file testEngineCopy.cc.

VERBOSER2

#define VERBOSER2

Definition at line 27 of file testEngineCopy.cc.

Function Documentation

aSequence()

std::vector< double > aSequence

(

int

n

)

Definition at line 47 of file testEngineCopy.cc.

                                   {
  std::vector<double> v;
  DualRand e(13542);
  RandFlat f(e);
  for (int i=0; i<n; i++) {
    v.push_back(f()); 
  }
  return v;
}

Referenced by vectorRestore1< NonRandomEngine >(), and vectorTest64< NonRandomEngine >().

equals()

bool equals	(	double	a,
		double	b
	)

Definition at line 41 of file testEngineCopy.cc.

                                {
  std::vector<double> ab(2);
  ab[0]=a;  ab[1]=b;
  return (equals01(ab));
}

Referenced by vectorTest64().

equals01()

bool equals01

(

const std::vector< double > &

ab

)

Definition at line 38 of file testEngineCopy.cc.

                                           {
  return ab[1]==ab[0];
}  

Referenced by equals().

main()

int main

(

)

Definition at line 193 of file testEngineCopy.cc.

           {
  int stat = 0;
 
#ifdef TEST_ENGINE_COPY
  output << "\n=================================\n";
  output << "         Part IX \n";
  output << "    Copy test of engines\n";
  output << "=================================\n\n";
 
  stat |= vectorRestore<DualRand>(113);
  // copies of DRand48Engine are not allowed
  //stat |= vectorRestore<DRand48Engine>(114);
  stat |= vectorRestore<Hurd160Engine>(115);
  stat |= vectorRestore<Hurd288Engine>(116);
  stat |= vectorRestore<HepJamesRandom>(117);
  stat |= vectorRestore<MixMaxRng>(149);
  stat |= vectorRestore<MTwistEngine>(118);
  stat |= vectorRestore<RanecuEngine>(139);
  stat |= vectorRestore<Ranlux64Engine>(119);
  stat |= vectorRestore<RanluxEngine>(120);
  stat |= vectorRestore<RanluxppEngine>(119);
  stat |= vectorRestore<RanshiEngine>(121);
  stat |= vectorRestore<TripleRand>(122);
  stat |= vectorRestore<NonRandomEngine>(123);
  // anonymous engines are not copyable
  //stat |= vectorRestore<RandEngine>(129);
#endif
 
  output << "\n=============================================\n\n";
 
  if (stat != 0) {
     std::cout << "One or more problems detected: stat = " << stat << "\n";
     output << "One or more problems detected: stat = " << stat << "\n";
  }  else {
     output << "ranRestoreTest passed with no problems detected.\n";    
  }
 
  if (stat == 0) return 0;
  if (stat > 0) return -(stat|1);
  return stat|1;
}   

output()

std::ofstream output

(

"testEngineCopy.cout"

)

Referenced by main(), vectorRestore1(), vectorRestore1< NonRandomEngine >(), vectorRestore2(), vectorTest64(), and vectorTest64< NonRandomEngine >().

vectorRestore()

template<class E >

int vectorRestore

(

int

n

)

Definition at line 178 of file testEngineCopy.cc.

                         {
  std::vector<double> v;
  int status1 = vectorTest64<E>(n);
  E f = vectorRestore1<E>(n,v);
  int status2 = vectorRestore2<E>(f, v);  
  return (status1 | status2);  
}

vectorRestore1()

template<class E >

E vectorRestore1	(	int	n,
		std::vector< double > &	v
	)

Definition at line 118 of file testEngineCopy.cc.

                                               {
  output << "Copy for " << E::engineName() << "\n";
  E e(97538466);                    
  double r=0;
  for (int i=0; i<n; i++) r += e.flat();        
  E f(e);    
  for (int j=0; j<25; j++) v.push_back(e.flat());   
#ifdef VERBOSER2
  output << "First four of v are: " 
        << v[0] << ",  " << v[1] << ",  " << v[2] << ",  " << v[3] << "\n";
#endif
  return f;
}

vectorRestore1< NonRandomEngine >()

template<>

NonRandomEngine vectorRestore1< NonRandomEngine >	(	int	n,
		std::vector< double > &	v
	)

Definition at line 133 of file testEngineCopy.cc.

                                                               {
#ifdef VERBOSER2
  output << "Copy for " << NonRandomEngine::engineName() << "\n";
#endif
  std::vector<double> nonRand = aSequence(500);
  NonRandomEngine e; 
  e.setRandomSequence(&nonRand[0], (int)nonRand.size());
  double r=0;
  for (int i=0; i<n; i++) r += e.flat();
  NonRandomEngine f(e);
  for (int j=0; j<25; j++) v.push_back(e.flat()); 
#ifdef VERBOSER2
  output << "First four of v are: " 
        << v[0] << ",  " << v[1] << ",  " << v[2] << ",  " << v[3] << "\n";
#endif
  return f;
}

vectorRestore2()

template<class E >

int vectorRestore2	(	E &	f,
		const std::vector< double > &	v
	)

Definition at line 153 of file testEngineCopy.cc.

                                                       {
  int stat = 0;
  std::vector<double> k;
  for (int j=0; j<25; j++) k.push_back(f.flat());
#ifdef VERBOSER2
  output << "First four of k are: " 
        << k[0] << ",  " << k[1] << ",  " << k[2] << ",  " << k[3] << "\n";
#endif
  for (int m1=0; m1<25; m1++) {
    if ( v[m1] != k[m1] ) {
      std::cout << "???? Incorrect copy restored value for engine: " 
            << E::engineName() << "\n"; 
      #ifdef CLEAN_OUTPUT
      output << "???? Incorrect copy restored value for engine: " 
            << E::engineName() << "\n"; 
      #endif
      stat |= 1048576;
      return stat;
    }
  }
  return stat;       
}

vectorTest64()

template<class E >

int vectorTest64

(

int

n

)

Definition at line 61 of file testEngineCopy.cc.

                        {
  output << "Copy 64bit test for " << E::engineName() << "\n";
 
  E e;                  
  double x = 0; 
  for (int i=0; i<n; i++) x += e.flat();        
  E f( e );
  x = e.flat();
  output << "x = " << x << std::endl;
 
  double y = f.flat();
  output << "y = " << y << std::endl;
  if( x != y ) return n;
  
  for( int i=0; i<1000; ++i ) {
      x = e.flat();
      y = f.flat();
      if( !equals(x,y) ) {
          output << "i = " << i << " x, y " << x << " " << y
             << " vectorTest64 problem: e != f \n";
      return n+i;
      }
  }
 
  return 0;
}

vectorTest64< NonRandomEngine >()

template<>

int vectorTest64< NonRandomEngine >

(

int

n

)

Definition at line 89 of file testEngineCopy.cc.

                                         {
  output << "Copy 64bit test for " << NonRandomEngine::engineName() << "\n";
 
  std::vector<double> nonRand = aSequence(500);
  NonRandomEngine e; 
  e.setRandomSequence(&nonRand[0], (int)nonRand.size());
 
  double x = 0; 
  for (int i=0; i<n; i++) x += e.flat();        
  std::vector<unsigned long> v = e.put();
  NonRandomEngine f(e);
  x = e.flat();
  output << "x = " << x << std::endl;
 
  double y = f.flat();
  output << "y = " << y << std::endl;
  if( x != y ) return n;
  
  for( int i=0; i<300; ++i ) {
      if( e.flat() != f.flat() ) {
          output << "i = " << i << " vectorTest64 for NonRandomEngine problem: e != f \n";
      return n+i;
      }
  }
 
  return 0;
}