diff options
author | julie <julielangou@users.noreply.github.com> | 2011-10-06 06:53:11 +0000 |
---|---|---|
committer | julie <julielangou@users.noreply.github.com> | 2011-10-06 06:53:11 +0000 |
commit | e1d39294aee16fa6db9ba079b14442358217db71 (patch) | |
tree | 30e5aa04c1f6596991fda5334f63dfb9b8027849 /SRC/zsptri.f | |
parent | 5fe0466a14e395641f4f8a300ecc9dcb8058081b (diff) | |
download | lapack-e1d39294aee16fa6db9ba079b14442358217db71.tar.gz lapack-e1d39294aee16fa6db9ba079b14442358217db71.tar.bz2 lapack-e1d39294aee16fa6db9ba079b14442358217db71.zip |
Integrating Doxygen in comments
Diffstat (limited to 'SRC/zsptri.f')
-rw-r--r-- | SRC/zsptri.f | 149 |
1 files changed, 104 insertions, 45 deletions
diff --git a/SRC/zsptri.f b/SRC/zsptri.f index 3d41b61e..091e61a6 100644 --- a/SRC/zsptri.f +++ b/SRC/zsptri.f @@ -1,9 +1,111 @@ +*> \brief \b ZSPTRI +* +* =========== DOCUMENTATION =========== +* +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ +* +* Definition +* ========== +* +* SUBROUTINE ZSPTRI( UPLO, N, AP, IPIV, WORK, INFO ) +* +* .. Scalar Arguments .. +* CHARACTER UPLO +* INTEGER INFO, N +* .. +* .. Array Arguments .. +* INTEGER IPIV( * ) +* COMPLEX*16 AP( * ), WORK( * ) +* .. +* +* Purpose +* ======= +* +*>\details \b Purpose: +*>\verbatim +*> +*> ZSPTRI computes the inverse of a complex symmetric indefinite matrix +*> A in packed storage using the factorization A = U*D*U**T or +*> A = L*D*L**T computed by ZSPTRF. +*> +*>\endverbatim +* +* Arguments +* ========= +* +*> \param[in] UPLO +*> \verbatim +*> 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. +*> \endverbatim +*> +*> \param[in] N +*> \verbatim +*> N is INTEGER +*> The order of the matrix A. N >= 0. +*> \endverbatim +*> +*> \param[in,out] AP +*> \verbatim +*> AP is COMPLEX*16 array, dimension (N*(N+1)/2) +*> On entry, the block diagonal matrix D and the multipliers +*> used to obtain the factor U or L as computed by ZSPTRF, +*> stored as a packed triangular matrix. +*> \endverbatim +*> \verbatim +*> On exit, if INFO = 0, the (symmetric) inverse of the original +*> matrix, stored as a packed triangular matrix. The j-th column +*> of inv(A) is stored in the array AP as follows: +*> if UPLO = 'U', AP(i + (j-1)*j/2) = inv(A)(i,j) for 1<=i<=j; +*> if UPLO = 'L', +*> AP(i + (j-1)*(2n-j)/2) = inv(A)(i,j) for j<=i<=n. +*> \endverbatim +*> +*> \param[in] IPIV +*> \verbatim +*> IPIV is INTEGER array, dimension (N) +*> Details of the interchanges and the block structure of D +*> as determined by ZSPTRF. +*> \endverbatim +*> +*> \param[out] WORK +*> \verbatim +*> WORK is COMPLEX*16 array, dimension (N) +*> \endverbatim +*> +*> \param[out] INFO +*> \verbatim +*> 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. +*> \endverbatim +*> +* +* Authors +* ======= +* +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. +* +*> \date November 2011 +* +*> \ingroup complex16OTHERcomputational +* +* ===================================================================== SUBROUTINE ZSPTRI( UPLO, N, AP, IPIV, WORK, INFO ) * -* -- LAPACK routine (version 3.3.1) -- +* -- LAPACK computational routine (version 3.3.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* -- April 2011 -- +* November 2011 * * .. Scalar Arguments .. CHARACTER UPLO @@ -14,49 +116,6 @@ COMPLEX*16 AP( * ), WORK( * ) * .. * -* Purpose -* ======= -* -* ZSPTRI computes the inverse of a complex symmetric indefinite matrix -* A in packed storage using the factorization A = U*D*U**T or -* A = L*D*L**T computed by ZSPTRF. -* -* Arguments -* ========= -* -* UPLO (input) 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) INTEGER -* The order of the matrix A. N >= 0. -* -* AP (input/output) COMPLEX*16 array, dimension (N*(N+1)/2) -* On entry, the block diagonal matrix D and the multipliers -* used to obtain the factor U or L as computed by ZSPTRF, -* stored as a packed triangular matrix. -* -* On exit, if INFO = 0, the (symmetric) inverse of the original -* matrix, stored as a packed triangular matrix. The j-th column -* of inv(A) is stored in the array AP as follows: -* if UPLO = 'U', AP(i + (j-1)*j/2) = inv(A)(i,j) for 1<=i<=j; -* if UPLO = 'L', -* AP(i + (j-1)*(2n-j)/2) = inv(A)(i,j) for j<=i<=n. -* -* IPIV (input) INTEGER array, dimension (N) -* Details of the interchanges and the block structure of D -* as determined by ZSPTRF. -* -* WORK (workspace) COMPLEX*16 array, dimension (N) -* -* INFO (output) 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. -* * ===================================================================== * * .. Parameters .. |