diff --git a/test/cxx_lapacke/src/example.cpp b/test/cxx_lapacke/src/example.cpp index 78a20b0..df9e384 100644 --- a/test/cxx_lapacke/src/example.cpp +++ b/test/cxx_lapacke/src/example.cpp @@ -1,6 +1,7 @@ #include #include -#if defined HAVE_MKL_LAPACK + +#if defined HAVE_MKL_LAPACK #include #include #pragma message "Using Intel MKL c-lapack interface, ." @@ -14,12 +15,12 @@ #pragma message "Using ATLAS/SYSTEM_NATIVE C lapack interface, ." #endif -void main(void) +int main() { - int i,j,n=3; - double a[n*n],b[n],small = 1.0e-12; - int ierr,ipiv[n]; - int roots_ok=0; + const int n = 3; + + double a[n*n]; + double b[n]; a[0] = 2.00; a[1] = 1.00; @@ -35,18 +36,28 @@ void main(void) b[1] = 3.00; b[2] = 5.00; + int ierr; + int ipiv[n]; + /* MKL, OpenBLAS */ #if defined HAVE_MKL_LAPACK || defined HAVE_OPENBLAS_LAPACK ierr = LAPACKE_dgesv(CblasColMajor, n, 1, a, n, ipiv, b, n); #else /* ATLAS, SYSTEM_NATIVE */ ierr = clapack_dgesv(CblasColMajor, n, 1, a, n, ipiv, b, n); #endif - if (ierr != 0) { fprintf(stderr, "\nC dgesv failure with error %d\n", ierr);} + if (ierr != 0) + { + fprintf(stderr, "\ndgesv failure with error %i\n", ierr); + } - if (abs(b[0] + 0.500) <= small && - abs(b[1] - 0.250) <= small && - abs(b[2] - 0.250) <= small) roots_ok=1; + const double small = 1.0e-12; - if ( roots_ok == 1 ) {fprintf(stdout,"PASSED"); } - else {fprintf(stderr,"\nC LAPACK dgesv failure!");} + if (abs(b[0] + 0.50) <= small && + abs(b[1] - 0.25) <= small && + abs(b[2] - 0.25) <= small) + { + printf("PASSED"); + } + + return 0; }