version 1.7, 2011/07/22 07:38:03
|
version 1.10, 2012/12/14 14:22:04
|
Line 1
|
Line 1
|
|
*> \brief \b ZTPMV |
|
* |
|
* =========== DOCUMENTATION =========== |
|
* |
|
* Online html documentation available at |
|
* http://www.netlib.org/lapack/explore-html/ |
|
* |
|
* Definition: |
|
* =========== |
|
* |
|
* SUBROUTINE ZTPMV(UPLO,TRANS,DIAG,N,AP,X,INCX) |
|
* |
|
* .. Scalar Arguments .. |
|
* INTEGER INCX,N |
|
* CHARACTER DIAG,TRANS,UPLO |
|
* .. |
|
* .. Array Arguments .. |
|
* COMPLEX*16 AP(*),X(*) |
|
* .. |
|
* |
|
* |
|
*> \par Purpose: |
|
* ============= |
|
*> |
|
*> \verbatim |
|
*> |
|
*> ZTPMV performs one of the matrix-vector operations |
|
*> |
|
*> x := A*x, or x := A**T*x, or x := A**H*x, |
|
*> |
|
*> where x is an n element vector and A is an n by n unit, or non-unit, |
|
*> upper or lower triangular matrix, supplied in packed form. |
|
*> \endverbatim |
|
* |
|
* Arguments: |
|
* ========== |
|
* |
|
*> \param[in] UPLO |
|
*> \verbatim |
|
*> UPLO is CHARACTER*1 |
|
*> On entry, UPLO specifies whether the matrix is an upper or |
|
*> lower triangular matrix as follows: |
|
*> |
|
*> UPLO = 'U' or 'u' A is an upper triangular matrix. |
|
*> |
|
*> UPLO = 'L' or 'l' A is a lower triangular matrix. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] TRANS |
|
*> \verbatim |
|
*> TRANS is CHARACTER*1 |
|
*> On entry, TRANS specifies the operation to be performed as |
|
*> follows: |
|
*> |
|
*> TRANS = 'N' or 'n' x := A*x. |
|
*> |
|
*> TRANS = 'T' or 't' x := A**T*x. |
|
*> |
|
*> TRANS = 'C' or 'c' x := A**H*x. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] DIAG |
|
*> \verbatim |
|
*> DIAG is CHARACTER*1 |
|
*> On entry, DIAG specifies whether or not A is unit |
|
*> triangular as follows: |
|
*> |
|
*> DIAG = 'U' or 'u' A is assumed to be unit triangular. |
|
*> |
|
*> DIAG = 'N' or 'n' A is not assumed to be unit |
|
*> triangular. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] N |
|
*> \verbatim |
|
*> N is INTEGER |
|
*> On entry, N specifies the order of the matrix A. |
|
*> N must be at least zero. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] AP |
|
*> \verbatim |
|
*> AP is COMPLEX*16 array of DIMENSION at least |
|
*> ( ( n*( n + 1 ) )/2 ). |
|
*> Before entry with UPLO = 'U' or 'u', the array AP must |
|
*> contain the upper triangular matrix packed sequentially, |
|
*> column by column, so that AP( 1 ) contains a( 1, 1 ), |
|
*> AP( 2 ) and AP( 3 ) contain a( 1, 2 ) and a( 2, 2 ) |
|
*> respectively, and so on. |
|
*> Before entry with UPLO = 'L' or 'l', the array AP must |
|
*> contain the lower triangular matrix packed sequentially, |
|
*> column by column, so that AP( 1 ) contains a( 1, 1 ), |
|
*> AP( 2 ) and AP( 3 ) contain a( 2, 1 ) and a( 3, 1 ) |
|
*> respectively, and so on. |
|
*> Note that when DIAG = 'U' or 'u', the diagonal elements of |
|
*> A are not referenced, but are assumed to be unity. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] X |
|
*> \verbatim |
|
*> X is (input/output) COMPLEX*16 array of dimension at least |
|
*> ( 1 + ( n - 1 )*abs( INCX ) ). |
|
*> Before entry, the incremented array X must contain the n |
|
*> element vector x. On exit, X is overwritten with the |
|
*> tranformed vector x. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] INCX |
|
*> \verbatim |
|
*> INCX is INTEGER |
|
*> On entry, INCX specifies the increment for the elements of |
|
*> X. INCX must not be zero. |
|
*> \endverbatim |
|
* |
|
* Authors: |
|
* ======== |
|
* |
|
*> \author Univ. of Tennessee |
|
*> \author Univ. of California Berkeley |
|
*> \author Univ. of Colorado Denver |
|
*> \author NAG Ltd. |
|
* |
|
*> \date November 2011 |
|
* |
|
*> \ingroup complex16_blas_level2 |
|
* |
|
*> \par Further Details: |
|
* ===================== |
|
*> |
|
*> \verbatim |
|
*> |
|
*> Level 2 Blas routine. |
|
*> The vector and matrix arguments are not referenced when N = 0, or M = 0 |
|
*> |
|
*> -- Written on 22-October-1986. |
|
*> Jack Dongarra, Argonne National Lab. |
|
*> Jeremy Du Croz, Nag Central Office. |
|
*> Sven Hammarling, Nag Central Office. |
|
*> Richard Hanson, Sandia National Labs. |
|
*> \endverbatim |
|
*> |
|
* ===================================================================== |
SUBROUTINE ZTPMV(UPLO,TRANS,DIAG,N,AP,X,INCX) |
SUBROUTINE ZTPMV(UPLO,TRANS,DIAG,N,AP,X,INCX) |
|
* |
|
* -- Reference BLAS level2 routine (version 3.4.0) -- |
|
* -- Reference BLAS is a software package provided by Univ. of Tennessee, -- |
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
|
* November 2011 |
|
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INCX,N |
INTEGER INCX,N |
CHARACTER DIAG,TRANS,UPLO |
CHARACTER DIAG,TRANS,UPLO |
* .. |
* .. |
* .. Array Arguments .. |
* .. Array Arguments .. |
DOUBLE COMPLEX AP(*),X(*) |
COMPLEX*16 AP(*),X(*) |
* .. |
* .. |
* |
* |
* Purpose |
|
* ======= |
|
* |
|
* ZTPMV performs one of the matrix-vector operations |
|
* |
|
* x := A*x, or x := A**T*x, or x := A**H*x, |
|
* |
|
* where x is an n element vector and A is an n by n unit, or non-unit, |
|
* upper or lower triangular matrix, supplied in packed form. |
|
* |
|
* Arguments |
|
* ========== |
|
* |
|
* UPLO - CHARACTER*1. |
|
* On entry, UPLO specifies whether the matrix is an upper or |
|
* lower triangular matrix as follows: |
|
* |
|
* UPLO = 'U' or 'u' A is an upper triangular matrix. |
|
* |
|
* UPLO = 'L' or 'l' A is a lower triangular matrix. |
|
* |
|
* Unchanged on exit. |
|
* |
|
* TRANS - CHARACTER*1. |
|
* On entry, TRANS specifies the operation to be performed as |
|
* follows: |
|
* |
|
* TRANS = 'N' or 'n' x := A*x. |
|
* |
|
* TRANS = 'T' or 't' x := A**T*x. |
|
* |
|
* TRANS = 'C' or 'c' x := A**H*x. |
|
* |
|
* Unchanged on exit. |
|
* |
|
* DIAG - CHARACTER*1. |
|
* On entry, DIAG specifies whether or not A is unit |
|
* triangular as follows: |
|
* |
|
* DIAG = 'U' or 'u' A is assumed to be unit triangular. |
|
* |
|
* DIAG = 'N' or 'n' A is not assumed to be unit |
|
* triangular. |
|
* |
|
* Unchanged on exit. |
|
* |
|
* N - INTEGER. |
|
* On entry, N specifies the order of the matrix A. |
|
* N must be at least zero. |
|
* Unchanged on exit. |
|
* |
|
* AP - COMPLEX*16 array of DIMENSION at least |
|
* ( ( n*( n + 1 ) )/2 ). |
|
* Before entry with UPLO = 'U' or 'u', the array AP must |
|
* contain the upper triangular matrix packed sequentially, |
|
* column by column, so that AP( 1 ) contains a( 1, 1 ), |
|
* AP( 2 ) and AP( 3 ) contain a( 1, 2 ) and a( 2, 2 ) |
|
* respectively, and so on. |
|
* Before entry with UPLO = 'L' or 'l', the array AP must |
|
* contain the lower triangular matrix packed sequentially, |
|
* column by column, so that AP( 1 ) contains a( 1, 1 ), |
|
* AP( 2 ) and AP( 3 ) contain a( 2, 1 ) and a( 3, 1 ) |
|
* respectively, and so on. |
|
* Note that when DIAG = 'U' or 'u', the diagonal elements of |
|
* A are not referenced, but are assumed to be unity. |
|
* Unchanged on exit. |
|
* |
|
* X - COMPLEX*16 array of dimension at least |
|
* ( 1 + ( n - 1 )*abs( INCX ) ). |
|
* Before entry, the incremented array X must contain the n |
|
* element vector x. On exit, X is overwritten with the |
|
* tranformed vector x. |
|
* |
|
* INCX - INTEGER. |
|
* On entry, INCX specifies the increment for the elements of |
|
* X. INCX must not be zero. |
|
* Unchanged on exit. |
|
* |
|
* Further Details |
|
* =============== |
|
* |
|
* Level 2 Blas routine. |
|
* The vector and matrix arguments are not referenced when N = 0, or M = 0 |
|
* |
|
* -- Written on 22-October-1986. |
|
* Jack Dongarra, Argonne National Lab. |
|
* Jeremy Du Croz, Nag Central Office. |
|
* Sven Hammarling, Nag Central Office. |
|
* Richard Hanson, Sandia National Labs. |
|
* |
|
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. Parameters .. |
* .. Parameters .. |
DOUBLE COMPLEX ZERO |
COMPLEX*16 ZERO |
PARAMETER (ZERO= (0.0D+0,0.0D+0)) |
PARAMETER (ZERO= (0.0D+0,0.0D+0)) |
* .. |
* .. |
* .. Local Scalars .. |
* .. Local Scalars .. |
DOUBLE COMPLEX TEMP |
COMPLEX*16 TEMP |
INTEGER I,INFO,IX,J,JX,K,KK,KX |
INTEGER I,INFO,IX,J,JX,K,KK,KX |
LOGICAL NOCONJ,NOUNIT |
LOGICAL NOCONJ,NOUNIT |
* .. |
* .. |