diff options
Diffstat (limited to 'src/c/linearAlgebra/lu/dlua.c')
| -rw-r--r-- | src/c/linearAlgebra/lu/dlua.c | 9 |
1 files changed, 2 insertions, 7 deletions
diff --git a/src/c/linearAlgebra/lu/dlua.c b/src/c/linearAlgebra/lu/dlua.c index d3cf6d11..a5c4f5ff 100644 --- a/src/c/linearAlgebra/lu/dlua.c +++ b/src/c/linearAlgebra/lu/dlua.c @@ -21,12 +21,9 @@ extern double dgetrf_ (int* , int* , double* , int* , int* , int* ); void dlua(double*inp1, int size, double* out1, double* out2) { - char TRANS = 'N'; int INFO=3; int LDA = 3; - int LDB = 3; int N = 3; - int NRHS = 1; int IPIV[2] ; /*double A[9] = @@ -37,15 +34,13 @@ void dlua(double*inp1, int size, double* out1, double* out2) };*/ //void LAPACK_dgetrf( lapack_int* m, lapack_int* n, double* a, lapack_int* lda, lapack_int* ipiv, lapack_int *info ); -// dgetrf_(&N,&N,inp1,&LDA,IPIV,&INFO); + dgetrf_(&N,&N,inp1,&LDA,IPIV,&INFO); - dgetrf_(6,2,inp1,6,IPIV,&INFO); - // checks INFO, if INFO != 0 something goes wrong, for more information see the MAN page of dgetrf. for(int i = 0; i < 3*3; i++) printf("%lf \n ", inp1[i]); - printf("PIVOTTTTTT \n "); + printf("PIVOTTTTTT \n "); for(int i = 0; i < 2; i++) printf("%d \n ", IPIV[i]); |