Commit a9f971db authored by Martin Perdacher's avatar Martin Perdacher

Update dataIo.cpp

parent 5b74f9c4
...@@ -9,11 +9,10 @@ void random_init(double *array, const int N, const int D){ ...@@ -9,11 +9,10 @@ void random_init(double *array, const int N, const int D){
seed[1]=(((tId*tId + 13) * 2)/3); seed[1]=(((tId*tId + 13) * 2)/3);
seed[2]=tId; seed[2]=tId;
#pragma omp for #pragma omp for collapse(2)
for (int i = 0; i < N; i++){ for (int i = 0; i < N; i++){
array[i*N+i] = erand48(seed) * N / 10.0 + N ; // ensure positive definiteness for (int j = 0; j < D; j++)
for (int j = 0; j < i; j++) array[i*N+j] = erand48(seed);
array[i*N+j] = array[j*N+i] = erand48(seed);
} }
} }
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment