*> \brief \b ZLAT2C converts a double complex triangular matrix to a complex triangular matrix.
*
* =========== DOCUMENTATION ===========
*
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
*> Download ZLAT2C + dependencies
*>
*> [TGZ]
*>
*> [ZIP]
*>
*> [TXT]
*> \endhtmlonly
*
* Definition:
* ===========
*
* SUBROUTINE ZLAT2C( UPLO, N, A, LDA, SA, LDSA, INFO )
*
* .. Scalar Arguments ..
* CHARACTER UPLO
* INTEGER INFO, LDA, LDSA, N
* ..
* .. Array Arguments ..
* COMPLEX SA( LDSA, * )
* COMPLEX*16 A( LDA, * )
* ..
*
*
*> \par Purpose:
* =============
*>
*> \verbatim
*>
*> ZLAT2C converts a COMPLEX*16 triangular matrix, SA, to a COMPLEX
*> triangular matrix, A.
*>
*> RMAX is the overflow for the SINGLE PRECISION arithmetic
*> ZLAT2C checks that all the entries of A are between -RMAX and
*> RMAX. If not the convertion is aborted and a flag is raised.
*>
*> This is an auxiliary routine so there is no argument checking.
*> \endverbatim
*
* Arguments:
* ==========
*
*> \param[in] UPLO
*> \verbatim
*> UPLO is CHARACTER*1
*> = 'U': A is upper triangular;
*> = 'L': A is lower triangular.
*> \endverbatim
*>
*> \param[in] N
*> \verbatim
*> N is INTEGER
*> The number of rows and columns of the matrix A. N >= 0.
*> \endverbatim
*>
*> \param[in] A
*> \verbatim
*> A is COMPLEX*16 array, dimension (LDA,N)
*> On entry, the N-by-N triangular coefficient matrix A.
*> \endverbatim
*>
*> \param[in] LDA
*> \verbatim
*> LDA is INTEGER
*> The leading dimension of the array A. LDA >= max(1,N).
*> \endverbatim
*>
*> \param[out] SA
*> \verbatim
*> SA is COMPLEX array, dimension (LDSA,N)
*> Only the UPLO part of SA is referenced. On exit, if INFO=0,
*> the N-by-N coefficient matrix SA; if INFO>0, the content of
*> the UPLO part of SA is unspecified.
*> \endverbatim
*>
*> \param[in] LDSA
*> \verbatim
*> LDSA is INTEGER
*> The leading dimension of the array SA. LDSA >= max(1,M).
*> \endverbatim
*>
*> \param[out] INFO
*> \verbatim
*> INFO is INTEGER
*> = 0: successful exit.
*> = 1: an entry of the matrix A is greater than the SINGLE
*> PRECISION overflow threshold, in this case, the content
*> of the UPLO part of SA in exit is unspecified.
*> \endverbatim
*
* Authors:
* ========
*
*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
*> \date September 2012
*
*> \ingroup complex16OTHERauxiliary
*
* =====================================================================
SUBROUTINE ZLAT2C( UPLO, N, A, LDA, SA, LDSA, INFO )
*
* -- LAPACK auxiliary routine (version 3.4.2) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* September 2012
*
* .. Scalar Arguments ..
CHARACTER UPLO
INTEGER INFO, LDA, LDSA, N
* ..
* .. Array Arguments ..
COMPLEX SA( LDSA, * )
COMPLEX*16 A( LDA, * )
* ..
*
* =====================================================================
*
* .. Local Scalars ..
INTEGER I, J
DOUBLE PRECISION RMAX
LOGICAL UPPER
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE, DIMAG
* ..
* .. External Functions ..
REAL SLAMCH
LOGICAL LSAME
EXTERNAL SLAMCH, LSAME
* ..
* .. Executable Statements ..
*
RMAX = SLAMCH( 'O' )
UPPER = LSAME( UPLO, 'U' )
IF( UPPER ) THEN
DO 20 J = 1, N
DO 10 I = 1, J
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR.
$ ( DBLE( A( I, J ) ).GT.RMAX ) .OR.
$ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR.
$ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN
INFO = 1
GO TO 50
END IF
SA( I, J ) = A( I, J )
10 CONTINUE
20 CONTINUE
ELSE
DO 40 J = 1, N
DO 30 I = J, N
IF( ( DBLE( A( I, J ) ).LT.-RMAX ) .OR.
$ ( DBLE( A( I, J ) ).GT.RMAX ) .OR.
$ ( DIMAG( A( I, J ) ).LT.-RMAX ) .OR.
$ ( DIMAG( A( I, J ) ).GT.RMAX ) ) THEN
INFO = 1
GO TO 50
END IF
SA( I, J ) = A( I, J )
30 CONTINUE
40 CONTINUE
END IF
50 CONTINUE
*
RETURN
*
* End of ZLAT2C
*
END