*> \brief \b ZLAG2C
*
* =========== DOCUMENTATION ===========
*
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
*> Download ZLAG2C + dependencies
*>
*> [TGZ]
*>
*> [ZIP]
*>
*> [TXT]
*> \endhtmlonly
*
* Definition:
* ===========
*
* SUBROUTINE ZLAG2C( M, N, A, LDA, SA, LDSA, INFO )
*
* .. Scalar Arguments ..
* INTEGER INFO, LDA, LDSA, M, N
* ..
* .. Array Arguments ..
* COMPLEX SA( LDSA, * )
* COMPLEX*16 A( LDA, * )
* ..
*
*
*> \par Purpose:
* =============
*>
*> \verbatim
*>
*> ZLAG2C converts a COMPLEX*16 matrix, SA, to a COMPLEX matrix, A.
*>
*> RMAX is the overflow for the SINGLE PRECISION arithmetic
*> ZLAG2C 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] M
*> \verbatim
*> M is INTEGER
*> The number of lines of the matrix A. M >= 0.
*> \endverbatim
*>
*> \param[in] N
*> \verbatim
*> N is INTEGER
*> The number of columns of the matrix A. N >= 0.
*> \endverbatim
*>
*> \param[in] A
*> \verbatim
*> A is COMPLEX*16 array, dimension (LDA,N)
*> On entry, the M-by-N coefficient matrix A.
*> \endverbatim
*>
*> \param[in] LDA
*> \verbatim
*> LDA is INTEGER
*> The leading dimension of the array A. LDA >= max(1,M).
*> \endverbatim
*>
*> \param[out] SA
*> \verbatim
*> SA is COMPLEX array, dimension (LDSA,N)
*> On exit, if INFO=0, the M-by-N coefficient matrix SA; if
*> INFO>0, the content 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 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 November 2011
*
*> \ingroup complex16OTHERauxiliary
*
* =====================================================================
SUBROUTINE ZLAG2C( M, N, A, LDA, SA, LDSA, INFO )
*
* -- LAPACK auxiliary routine (version 3.4.0) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* November 2011
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDSA, M, N
* ..
* .. Array Arguments ..
COMPLEX SA( LDSA, * )
COMPLEX*16 A( LDA, * )
* ..
*
* =====================================================================
*
* .. Local Scalars ..
INTEGER I, J
DOUBLE PRECISION RMAX
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE, DIMAG
* ..
* .. External Functions ..
REAL SLAMCH
EXTERNAL SLAMCH
* ..
* .. Executable Statements ..
*
RMAX = SLAMCH( 'O' )
DO 20 J = 1, N
DO 10 I = 1, M
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 30
END IF
SA( I, J ) = A( I, J )
10 CONTINUE
20 CONTINUE
INFO = 0
30 CONTINUE
RETURN
*
* End of ZLAG2C
*
END