Dinv_triang Subroutine

subroutine Dinv_triang(A, uplo, diag)

Arguments

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

Calls

proc~~dinv_triang~2~~CallsGraph proc~dinv_triang~2 Dinv_triang dtrtri dtrtri proc~dinv_triang~2->dtrtri

Source Code

subroutine Dinv_triang(A,uplo,diag)
  real(8),dimension(:,:)           :: A
  character(len=1),optional        :: uplo,diag
  character(len=1)                 :: uplo_
  character(len=1)                 :: diag_
  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 dtrtri(uplo_,diag_,n,A,lda,info)
  if(info/=0)stop "Error MATRIX/D_mat_invertTRIANG: dtrtri"
end subroutine Dinv_triang