cpp_thread_test/cpp_thread_safety_common.h

   1 inline void pauser(){
   2     /// a portable way to pause a program
   3     std::string dummy;
   4     std::cout << "Press enter to continue...";
   5     std::getline(std::cin, dummy);
   6 }
   7
   8 void FillMatrices(std::vector<std::vector<double>>& matBlock, std::mt19937_64& PRNG, std::uniform_real_distribution<double>& rngdist, const blasint randomMatSize, const uint32_t numConcurrentThreads, const uint32_t numMat){
   9         for(uint32_t i=0; i<numMat; i++){
  10                 for(uint32_t j = 0; j < static_cast<uint32_t>(randomMatSize*randomMatSize); j++){
  11                         matBlock[i][j] = rngdist(PRNG);
  12                 }
  13         }
  14         for(uint32_t i=numMat; i<(numConcurrentThreads*numMat); i+=numMat){
  15                 for(uint32_t j=0; j<numMat; j++){
  16                         matBlock[i+j] = matBlock[j];
  17                 }
  18         }
  19 }
  20
  21 void FillVectors(std::vector<std::vector<double>>& vecBlock, std::mt19937_64& PRNG, std::uniform_real_distribution<double>& rngdist, const blasint randomMatSize, const uint32_t numConcurrentThreads, const uint32_t numVec){
  22         for(uint32_t i=0; i<numVec; i++){
  23                 for(uint32_t j = 0; j < static_cast<uint32_t>(randomMatSize); j++){
  24                         vecBlock[i][j] = rngdist(PRNG);
  25                 }
  26         }
  27         for(uint32_t i=numVec; i<(numConcurrentThreads*numVec); i+=numVec){
  28                 for(uint32_t j=0; j<numVec; j++){
  29                         vecBlock[i+j] = vecBlock[j];
  30                 }
  31         }
  32 }
  33
  34 std::mt19937_64 InitPRNG(){
  35         std::random_device rd;
  36         std::mt19937_64 PRNG(rd()); //seed PRNG using /dev/urandom or similar OS provided RNG
  37         std::uniform_real_distribution<double> rngdist{-1.0, 1.0};
  38         //make sure the internal state of the PRNG is properly mixed by generating 10M random numbers
  39         //PRNGs often have unreliable distribution uniformity and other statistical properties before their internal state is sufficiently mixed
  40         for (uint32_t i=0;i<10000000;i++) rngdist(PRNG);
  41         return PRNG;
  42 }
  43
  44 void PrintMatrices(const std::vector<std::vector<double>>& matBlock, const blasint randomMatSize, const uint32_t numConcurrentThreads, const uint32_t numMat){
  45         for (uint32_t i=0;i<numConcurrentThreads*numMat;i++){
  46                 std::cout<<i<<std::endl;
  47                 for (uint32_t j = 0; j < static_cast<uint32_t>(randomMatSize); j++){
  48                         for (uint32_t k = 0; k < static_cast<uint32_t>(randomMatSize); k++){
  49                                 std::cout<<matBlock[i][j*randomMatSize + k]<<"  ";
  50                         }
  51                         std::cout<<std::endl;
  52                 }
  53                 std::cout<<std::endl;
  54         }
  55 }