zhetri2 (3p)

Name

zhetri2 - compute the inverse of a COMPLEX*16 Hermitian indefinite matrix A using the factorization A=U*D*U**T or A=L*D*L**T computed by ZHETRF

Synopsis

SUBROUTINE ZHETRI2(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


CHARACTER*1 UPLO

INTEGER INFO, LDA, LWORK, N

INTEGER IPIV(*)

DOUBLE COMPLEX A(LDA,*), WORK(*)


SUBROUTINE ZHETRI2_64(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


CHARACTER*1 UPLO

INTEGER*8 INFO, LDA, LWORK, N

INTEGER*8 IPIV(*)

DOUBLE COMPLEX A(LDA,*), WORK(*)


F95 INTERFACE
SUBROUTINE HETRI2(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


INTEGER :: N, LDA, LWORK, INFO

CHARACTER(LEN=1) :: UPLO

INTEGER, DIMENSION(:) :: IPIV

COMPLEX(8), DIMENSION(:,:) :: A

COMPLEX(8), DIMENSION(:) :: WORK


SUBROUTINE HETRI2_64(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


INTEGER(8) :: N, LDA, LWORK, INFO

CHARACTER(LEN=1) :: UPLO

INTEGER(8), DIMENSION(:) :: IPIV

COMPLEX(8), DIMENSION(:,:) :: A

COMPLEX(8), DIMENSION(:) :: WORK


C INTERFACE
#include <sunperf.h>

void  zhetri2  (char uplo, int n, doublecomplex *a, int lda, int *ipiv,
int *info);


void zhetri2_64 (char uplo, long n, doublecomplex *a,  long  lda,  long
*ipiv, long *info);

Description

Oracle Solaris Studio Performance Library                          zhetri2(3P)



NAME
       zhetri2  -  compute  the  inverse  of a COMPLEX*16 Hermitian indefinite
       matrix A using the factorization A=U*D*U**T or A=L*D*L**T  computed  by
       ZHETRF


SYNOPSIS
       SUBROUTINE ZHETRI2(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


       CHARACTER*1 UPLO

       INTEGER INFO, LDA, LWORK, N

       INTEGER IPIV(*)

       DOUBLE COMPLEX A(LDA,*), WORK(*)


       SUBROUTINE ZHETRI2_64(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


       CHARACTER*1 UPLO

       INTEGER*8 INFO, LDA, LWORK, N

       INTEGER*8 IPIV(*)

       DOUBLE COMPLEX A(LDA,*), WORK(*)


   F95 INTERFACE
       SUBROUTINE HETRI2(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


       INTEGER :: N, LDA, LWORK, INFO

       CHARACTER(LEN=1) :: UPLO

       INTEGER, DIMENSION(:) :: IPIV

       COMPLEX(8), DIMENSION(:,:) :: A

       COMPLEX(8), DIMENSION(:) :: WORK


       SUBROUTINE HETRI2_64(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)


       INTEGER(8) :: N, LDA, LWORK, INFO

       CHARACTER(LEN=1) :: UPLO

       INTEGER(8), DIMENSION(:) :: IPIV

       COMPLEX(8), DIMENSION(:,:) :: A

       COMPLEX(8), DIMENSION(:) :: WORK


   C INTERFACE
       #include <sunperf.h>

       void  zhetri2  (char uplo, int n, doublecomplex *a, int lda, int *ipiv,
                 int *info);


       void zhetri2_64 (char uplo, long n, doublecomplex *a,  long  lda,  long
                 *ipiv, long *info);


PURPOSE
       zhetri2  computes  the  inverse  of  a  COMPLEX*16 hermitian indefinite
       matrix A using the factorization A = U*D*U**T or A = L*D*L**T  computed
       by  ZHETRF.  ZHETRI2  set the LEADING DIMENSION of the workspace before
       calling ZHETRI2X that actually computes the inverse.


ARGUMENTS
       UPLO (input)
                 UPLO is CHARACTER*1
                 Specifies whether the details of the factorization are stored
                 as an upper or lower triangular matrix.
                 = 'U':  Upper triangular, form is A=U*D*U**T;
                 = 'L':  Lower triangular, form is A=L*D*L**T.


       N (input)
                 N is INTEGER
                 The order of the matrix A.  N >= 0.


       A (input/output)
                 A is COMPLEX*16 array, dimension (LDA,N)
                 On  entry,  the NB diagonal matrix D and the multipliers used
                 to obtain the factor U or L as computed by ZHETRF.
                 On exit, if INFO = 0, the (symmetric) inverse of the original
                 matrix.   If  UPLO  =  'U',  the upper triangular part of the
                 inverse is formed and the part of A below the diagonal is not
                 referenced;  if  UPLO  = 'L' the lower triangular part of the
                 inverse is formed and the part of A above the diagonal is not
                 referenced.


       LDA (input)
                 LDA is INTEGER
                 The leading dimension of the array A.
                 LDA >= max(1,N).


       IPIV (input)
                 IPIV is INTEGER array, dimension (N)
                 Details  of  the  interchanges  and  the NB structure of D as
                 determined by ZHETRF.


       WORK (output)
                 WORK is COMPLEX*16 array, dimension (N+NB+1)*(NB+3)


       LWORK (input)
                 LWORK is INTEGER
                 The dimension of the array WORK.
                 WORK is size >= (N+NB+1)*(NB+3)
                 If LWORK = -1, then a workspace query is assumed; the routine
                 calculates:
                 -  the  optimal size of the WORK array, returns this value as
                 the first entry of the WORK array,
                 - and no error message related to LWORK is issued by  XERBLA.


       INFO (output)
                 INFO is INTEGER
                 = 0: successful exit
                 < 0: if INFO = -i, the i-th argument had an illegal value
                 >  0: if INFO = i, D(i,i) = 0; the matrix is singular and its
                 inverse could not be computed.



                                  7 Nov 2015                       zhetri2(3P)