info=LAPACKE_dgetrf(LAPACK_ROW_MAJOR,n,n,&(A[0]),n,&pivots[0]);info2=LAPACKE_dgetrs(LAPACK_ROW_MAJOR,'N',n,1,&A[0],n,&pivots[0],&b[0],1);Results: info=-4;info2=-5; that means problems on A. If I change the A initialization to:A[i*n+j]=//values (plus is ...
The only place in my code where I call LAPACKE_dgetrf and LAPACKE_dgetri is indet.c (which I also attach). Please note that although the functions are calledatlasdet, atlasinv, .. etc, only mkl routines are used. (The naming convention isbecause ...
lapack_int LAPACKE_dgetrf (int matrix_layout , lapack_int m , lapack_int n , double * a , lapack_int lda , lapack_int * ipiv ); lapack_int LAPACKE_dgetri (int matrix_layout , lapack_int n , double * a , lapack_int lda , const lapack_int ...