testVectorSave.cc File Reference

#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_VECTOR_ENGINE_RESTORE
#define	VERBOSER
#define	VERBOSER2

Functions
std::ofstream	output ("testVectorSave.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 >
std::vector< unsigned long >	vectorRestore1 (int n, std::vector< double > &v)
template<>
std::vector< unsigned long >	vectorRestore1< NonRandomEngine > (int n, std::vector< double > &v)
template<class E >
int	vectorRestore2 (const std::vector< unsigned long > state, 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 11 of file testVectorSave.cc.

TEST_VECTOR_ENGINE_RESTORE

#define TEST_VECTOR_ENGINE_RESTORE

Definition at line 23 of file testVectorSave.cc.

VERBOSER

#define VERBOSER

Definition at line 25 of file testVectorSave.cc.

VERBOSER2

#define VERBOSER2

Definition at line 26 of file testVectorSave.cc.

Function Documentation

aSequence()

std::vector< double > aSequence

(

int

n

)

Definition at line 40 of file testVectorSave.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 vectorTest64< NonRandomEngine >().

equals()

bool equals	(	double	a,
		double	b
	)

Definition at line 34 of file testVectorSave.cc.

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

equals01()

bool equals01

(

const std::vector< double > &

ab

)

Definition at line 31 of file testVectorSave.cc.

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

Referenced by equals().

main()

int main

(

)

Definition at line 206 of file testVectorSave.cc.

           {
  int stat = 0;
 
#ifdef TEST_VECTOR_ENGINE_RESTORE
  output << "\n=================================\n";
  output << "         Part IX \n";
  output << "Save/restore of engines to vectors\n";
  output << "=================================\n\n";
 
  stat |= vectorRestore<DualRand>(113);
  stat |= vectorRestore<DRand48Engine>(114);
  stat |= vectorRestore<Hurd160Engine>(115);
  stat |= vectorRestore<Hurd288Engine>(116);
  stat |= vectorRestore<HepJamesRandom>(117);
  stat |= vectorRestore<MTwistEngine>(118);
  stat |= vectorRestore<RanecuEngine>(139);
  stat |= vectorRestore<Ranlux64Engine>(119);
  stat |= vectorRestore<RanluxEngine>(120);
  stat |= vectorRestore<RanshiEngine>(121);
  stat |= vectorRestore<TripleRand>(122);
  stat |= vectorRestore<NonRandomEngine>(123);
  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

(

"testVectorSave.cout"

)

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

vectorRestore()

template<class E >

int vectorRestore

(

int

n

)

Definition at line 191 of file testVectorSave.cc.

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

vectorRestore1()

template<class E >

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

Definition at line 103 of file testVectorSave.cc.

                                                                      {
  output << "Vector restore for " << E::engineName() << "\n";
  E e(97538466);                    
  double r=0;                       
  for (int i=0; i<n; i++) r += e.flat();        
  std::vector<unsigned long> state = e.put();       
  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 state;
}

vectorRestore1< NonRandomEngine >()

template<>

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

Definition at line 103 of file testVectorSave.cc.

                                                               {
#ifdef VERBOSER2
  output << "Vector restore for " << NonRandomEngine::engineName() << "\n";
#endif
  std::vector<double> nonRand = aSequence(500);
  NonRandomEngine e; 
  e.setRandomSequence(&nonRand[0], (int)nonRand.size());
  double __attribute__ ((unused)) r = 0;
  for (int i=0; i<n; i++) r += e.flat();
  std::vector<unsigned long> state = e.put();       
  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 state;
}

vectorRestore2()

template<class E >

int vectorRestore2	(	const std::vector< unsigned long >	state,
		const std::vector< double > &	v
	)

Definition at line 138 of file testVectorSave.cc.

                                        {
  int stat = 0;
  std::vector<double> k;
  HepRandomEngine * a;
  a = HepRandomEngine::newEngine(state);
  if (!a) {
      std::cout << "???? could not restore engine state from vector for " 
            << E::engineName() << "\n"; 
      #ifdef CLEAN_OUTPUT
      output << "???? could not restore engine state from vector for " 
             << E::engineName() << "\n"; 
      #endif
      stat |= 1048576;
      return stat;
  }
  if (a->name() != E::engineName()) {
      #ifdef CLEAN_OUTPUT
      std::cout << "???? restored engine state from vector for " 
            << E::engineName() << "to different type of engine: "
        << a->name() << "\n"
    << "There is probably a clash in CRC hashes for these two names!\n"; 
      output << "???? restored engine state from vector for " 
            << E::engineName() << "to different type of engine: "
        << a->name() << "\n"
    << "There is probably a clash in CRC hashes for these two names!\n"; 
      #endif
      stat |= 1048576;
      return stat;
  }
  for (int j=0; j<25; j++) k.push_back(a->flat());
  delete a;
#ifdef VERBOSER2
  output << "First four of k are: " 
        << k[0] << ",  " << k[1] << ",  " << k[2] << ",  " << k[3] << "\n";
#endif
  for (int m=0; m<25; m++) {
    if ( v[m] != k[m] ) {
      std::cout << "???? Incorrect vector restored value for anonymous engine: " 
            << E::engineName() << "\n"; 
      #ifdef CLEAN_OUTPUT
      output << "???? Incorrect vector restored value for anonymous engine: " 
            << E::engineName() << "\n"; 
      #endif
      stat |= 1048576;
      return stat;
    }
  }
  return stat;       
}

vectorTest64()

template<class E >

int vectorTest64

(

int

n

)

Definition at line 56 of file testVectorSave.cc.

                        {
  output << "Vector restore 64bit test for " << E::engineName() << "\n";
 
  E e;                  
  double x = 0; 
  for (int i=0; i<n; i++) x += e.flat();        
  std::vector<unsigned long> v = e.put();
  x = e.flat();
  output << "x = " << x << std::endl;
 
  E f;
  v[0] &= 0xffffffffUL;
  f.get(v);
  double y = f.flat();
  output << "y = " << y << std::endl;
  if( x != y ) return n;
 
  return 0;
}

vectorTest64< NonRandomEngine >()

template<>

int vectorTest64< NonRandomEngine >

(

int

n

)

Definition at line 77 of file testVectorSave.cc.

                                         {
  output << "Vector restore 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();
  x = e.flat();
  output << "x = " << x << std::endl;
 
  NonRandomEngine f;
  f.setRandomSequence(&nonRand[0], (int)nonRand.size());
 
  v[0] &= 0xffffffffUL;
  f.get(v);
  double y = f.flat();
  output << "y = " << y << std::endl;
  if( x != y ) return n;
 
  return 0;
}