Computes the inverse of a triangular matrix stored in the Rectangular Full Packed (RFP) format.
FORTRAN 77:
call stftri( transr, uplo, diag, n, a, info )
call dtftri( transr, uplo, diag, n, a, info )
call ctftri( transr, uplo, diag, n, a, info )
call ztftri( transr, uplo, diag, n, a, info )
C:
lapack_int LAPACKE_<?>tftri( int matrix_order, char transr, char uplo, char diag, lapack_int n, <datatype>* a );
The FORTRAN 77 interfaces are specified in the mkl_lapack.fi and mkl_lapack.h include files, the Fortran 95 interfaces are specified in the lapack.f90 include file, and the C interfaces are specified in the mkl_lapacke.h include file.
Computes the inverse of a triangular matrix A stored in the Rectangular Full Packed (RFP) format. For the description of the RFP format, see Matrix Storage Schemes.
This is the block version of the algorithm, calling Level 3 BLAS.
The data types are given for the Fortran interface. A <datatype> placeholder, if present, is used for the C interface data types in the C interface section above. See C Interface Conventions for the C interface principal conventions and type defintions.
CHARACTER*1. Must be 'N', 'T' (for real data) or 'C' (for complex data).
If transr = 'N', the Normal transr of RFP A is stored.
If transr = 'T', the Transpose transr of RFP A is stored.
If transr = 'C', the Conjugate-Transpose transr of RFP A is stored.
CHARACTER*1. Must be 'U' or 'L'.
Indicates whether the upper or lower triangular part of RFP A is stored:
If uplo = 'U', the array a stores the upper triangular part of the matrix A.
If uplo = 'L', the array a stores the lower triangular part of the matrix A.
CHARACTER*1. Must be 'N' or 'U'.
If diag = 'N', then A is not a unit triangular matrix.
If diag = 'U', A is unit triangular: diagonal elements of A are assumed to be 1 and not referenced in the array a.
INTEGER. The order of the matrix A; n ≥ 0.
REAL for stftri
DOUBLE PRECISION for dtftri
COMPLEX for ctftri
DOUBLE COMPLEX for ztftri.
Array, DIMENSION (n*(n+1)/2). The array a contains the matrix A in the RFP format.
The (triangular) inverse of the original matrix in the same storage format.
INTEGER. If info=0, the execution is successful.
If info = -i, the i-th parameter had an illegal value.
If info = i, A(i,i) is exactly zero. The triangular matrix is singular and its inverse cannot be computed.
Copyright © 1994 - 2010, Intel Corporation. All rights reserved.