floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Handle OPENBLAS_LOOPS for more stable results

This commit is contained in:
Martin Kroeker 2021-07-01 17:39:37 +02:00 committed by GitHub
parent 06e3b07ecb
commit dcfc5cf714
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 13 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
benchmark/linpack.c
View File

@ -72,17 +72,21 @@ int main(int argc, char *argv[]){
FLOAT *a, *b; FLOAT *a, *b;
blasint *ipiv; blasint *ipiv;
blasint m, i, j, info; blasint m, i, j, l, info;
blasint unit = 1; blasint unit = 1;
int from = 1; int from = 1;
int to = 200; int to = 200;
int step = 1; int step = 1;
int loops = 1;
FLOAT maxerr; FLOAT maxerr;
double time1, time2; double time1, time2, timeg1,timeg2;
char *p;
if ((p = getenv("OPENBLAS_LOOPS"))) loops=*p;
argc--;argv++; argc--;argv++;
if (argc > 0) { from = atol(*argv); argc--; argv++;} if (argc > 0) { from = atol(*argv); argc--; argv++;}
@ -110,9 +114,9 @@ int main(int argc, char *argv[]){
fprintf(stderr, " SIZE Residual Decompose Solve Total\n"); fprintf(stderr, " SIZE Residual Decompose Solve Total\n");
for(m = from; m <= to; m += step){ for(m = from; m <= to; m += step){
timeg1 = timeg2 = 0.;
fprintf(stderr, " %6d : ", (int)m); fprintf(stderr, " %6d : ", (int)m);
for (l = 0; l < loops; l++) {
for(j = 0; j < m; j++){ for(j = 0; j < m; j++){
for(i = 0; i < m * COMPSIZE; i++){ for(i = 0; i < m * COMPSIZE; i++){
a[(long)i + (long)j * (long)m * COMPSIZE] = ((FLOAT) rand() / (FLOAT) RAND_MAX) - 0.5; a[(long)i + (long)j * (long)m * COMPSIZE] = ((FLOAT) rand() / (FLOAT) RAND_MAX) - 0.5;
@ -138,7 +142,7 @@ int main(int argc, char *argv[]){
exit(1); exit(1);
} }
time1 = getsec(); timeg1 += getsec();
begin(); begin();
@ -151,8 +155,10 @@ int main(int argc, char *argv[]){
exit(1); exit(1);
} }
time2 = getsec(); timeg2 += getsec();
} //loops
time1=timeg1/(double)loops;
time2=timeg2/(double)loops;
maxerr = 0.; maxerr = 0.;
for(i = 0; i < m; i++){ for(i = 0; i < m; i++){