after the call to dgetri an array (CIC) already set to zero has nonzero components
cic=0.0d0
lda=10
lwork=10
matr_(1:10,1:10)=pij_(0:9,0:9)
!lda=10
!lwork=-1
call dgetrf( 10, 10, matr_, lda, ipiv, info )
!lwork=640
call dgetri( 10, matr_, lda, ipiv, work, lwork, info )
continue
is there any explanation?
Thanks a lot
Eddy