dsytri2 (3p)

Name

dsytri2 - nite matrix A using the factorization A = U*D*U**T or A = L*D*L**T com- puted by DSYTRF

Synopsis

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


CHARACTER*1 UPLO

INTEGER INFO, LDA, LWORK, N

INTEGER IPIV(*)

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


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


CHARACTER*1 UPLO

INTEGER*8 INFO, LDA, LWORK, N

INTEGER*8 IPIV(*)

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


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


INTEGER :: N, LDA, LWORK, INFO

CHARACTER(LEN=1) :: UPLO

INTEGER, DIMENSION(:) :: IPIV

REAL(8), DIMENSION(:,:) :: A

REAL(8), DIMENSION(:) :: WORK


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


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

CHARACTER(LEN=1) :: UPLO

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

REAL(8), DIMENSION(:,:) :: A

REAL(8), DIMENSION(:) :: WORK


C INTERFACE
#include <sunperf.h>

void  dsytri2  (char  uplo,  int  n, double *a, int lda, int *ipiv, int
*info);


void dsytri2_64 (char uplo, long n, double *a, long  lda,  long  *ipiv,
long *info);

Description

Oracle Solaris Studio Performance Library                          dsytri2(3P)



NAME
       dsytri2  -  compute the inverse of a DOUBLE PRECISION symmetric indefi-
       nite matrix A using the factorization A = U*D*U**T or A = L*D*L**T com-
       puted by DSYTRF


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


       CHARACTER*1 UPLO

       INTEGER INFO, LDA, LWORK, N

       INTEGER IPIV(*)

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


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


       CHARACTER*1 UPLO

       INTEGER*8 INFO, LDA, LWORK, N

       INTEGER*8 IPIV(*)

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


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


       INTEGER :: N, LDA, LWORK, INFO

       CHARACTER(LEN=1) :: UPLO

       INTEGER, DIMENSION(:) :: IPIV

       REAL(8), DIMENSION(:,:) :: A

       REAL(8), DIMENSION(:) :: WORK


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


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

       CHARACTER(LEN=1) :: UPLO

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

       REAL(8), DIMENSION(:,:) :: A

       REAL(8), DIMENSION(:) :: WORK


   C INTERFACE
       #include <sunperf.h>

       void  dsytri2  (char  uplo,  int  n, double *a, int lda, int *ipiv, int
                 *info);


       void dsytri2_64 (char uplo, long n, double *a, long  lda,  long  *ipiv,
                 long *info);


PURPOSE
       dsytri2 computes the inverse of a DOUBLE PRECISION symmetric indefinite
       matrix A using the factorization A = U*D*U**T or A = L*D*L**T  computed
       by  DSYTRF.  DSYTRI2 sets the LEADING DIMENSION of the workspace before
       calling DSYTRI2X 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 DOUBLE PRECISION 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 DSYTRF.
                 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 DSYTRF.


       WORK (output)
                 WORK is DOUBLE PRECISION 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 LDWORK = -1, then a workspace query is assumed;  the  rou-
                 tine 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 LDWORK 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                       dsytri2(3P)