2 /// a portable way to pause a program
4 std::cout << "Press enter to continue...";
5 std::getline(std::cin, dummy);
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);
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];
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);
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];
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);
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]<<" ";