name wgedi(3f) - [m_la] computes the determinant and inverse of a matrix using the factors computed by wgeco(3f) or wgefa(3f). synopsis subroutine ml_wgedi(ar,ai,lda,n,ipvt,detr,deti,workr,worki,job)
integer(kind=4) :: lda real(kind=8) :: ar(lda,*) real(kind=8) :: ai(lda,*) integer(kind=4) :: n integer(kind=4) :: ipvt(*) real(kind=8) :: detr(2) real(kind=8) :: deti(2) real(kind=8) :: workr(*) real(kind=8) :: worki(*) integer(kind=4) :: jobdescription wgedi(3f) computes the determinant and inverse of a matrix using the factors computed by wgeco(3f) or wgefa(3f).
on entry
a double-complex(lda, n) the output from wgeco or wgefa.on returnlda integer the leading dimension of the array a.
n integer the order of the matrix a.
ipvt integer(n) the pivot vector from wgeco(3f) or wgefa(3f).
work double-complex(n) work vector. contents destroyed.
job integer
= 11 both determinant and inverse. = 01 inverse only. = 10 determinant only.
a inverse of original matrix if requested. otherwise unchanged.error conditiondet double-complex(2) determinant of original matrix if requested. otherwise not referenced.
determinant = det(1) * 10.0**det(2) with 1.0 .le. cabs1(det(1) .lt. 10.0 or det(1) .eq. 0.0 .
a division by zero will occur if the input factor contains a zero on the diagonal and the inverse is requested. it will not occur if the subroutines are called correctly and if wgeco(3f) has set rcond .gt. 0.0 or wgefa(3f) has set info .eq. 0 .subroutines and functionslinpack. this version dated 07/01/79 . cleve moler, university of new mexico, argonne national lab.
blas waxpy,mat_wscal,mat_wswap fortran dabs,mod
Nemo Release 3.1 | wgedi (3) | June 29, 2025 |