Zinv_triang Subroutine

subroutine Zinv_triang(A, uplo, diag)

Arguments

Type IntentOptional Attributes Name
complex(kind=8), dimension(:,:) :: A
character(len=1), optional :: uplo
character(len=1), optional :: diag

Calls

proc~~zinv_triang~2~~CallsGraph proc~zinv_triang~2 Zinv_triang ztrtri ztrtri proc~zinv_triang~2->ztrtri

Source Code

subroutine Zinv_triang(A,uplo,diag)
  complex(8),dimension(:,:)           :: A
  character(len=1),optional        :: uplo,diag
  character(len=1)                 :: uplo_
  character(len=1)                 :: diag_
  integer                          :: ndim1,ndim2
  integer                          :: n,lda,info
  uplo_="U";if(present(uplo))uplo_=uplo
  diag_="N";if(present(diag))diag_=diag !not a unit triangular matrix
  lda = max(1,size(A,1))
  n   = size(A,2)
  call ztrtri(uplo_,diag_,n,A,lda,info)
  if(info/=0)stop "Error MATRIX/D_mat_invertTRIANG: ztrtri"
end subroutine Zinv_triang