version 1.9, 2011/11/21 20:43:08
|
version 1.19, 2018/05/29 06:55:22
|
Line 2
|
Line 2
|
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
* Online html documentation available at |
* Online html documentation available at |
* http://www.netlib.org/lapack/explore-html/ |
* http://www.netlib.org/lapack/explore-html/ |
* |
* |
*> \htmlonly |
*> \htmlonly |
*> Download ZGEBAL + dependencies |
*> Download ZGEBAL + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgebal.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgebal.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgebal.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgebal.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgebal.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgebal.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO ) |
* SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* CHARACTER JOB |
* CHARACTER JOB |
* INTEGER IHI, ILO, INFO, LDA, N |
* INTEGER IHI, ILO, INFO, LDA, N |
Line 28
|
Line 28
|
* DOUBLE PRECISION SCALE( * ) |
* DOUBLE PRECISION SCALE( * ) |
* COMPLEX*16 A( LDA, * ) |
* COMPLEX*16 A( LDA, * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 83
|
Line 83
|
*> |
*> |
*> \param[out] ILO |
*> \param[out] ILO |
*> \verbatim |
*> \verbatim |
|
*> ILO is INTEGER |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[out] IHI |
*> \param[out] IHI |
*> \verbatim |
*> \verbatim |
|
*> IHI is INTEGER |
*> ILO and IHI are set to INTEGER such that on exit |
*> ILO and IHI are set to INTEGER such that on exit |
*> A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N. |
*> A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N. |
*> If JOB = 'N' or 'S', ILO = 1 and IHI = N. |
*> If JOB = 'N' or 'S', ILO = 1 and IHI = N. |
Line 116
|
Line 118
|
* Authors: |
* Authors: |
* ======== |
* ======== |
* |
* |
*> \author Univ. of Tennessee |
*> \author Univ. of Tennessee |
*> \author Univ. of California Berkeley |
*> \author Univ. of California Berkeley |
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date June 2017 |
* |
* |
*> \ingroup complex16GEcomputational |
*> \ingroup complex16GEcomputational |
* |
* |
Line 160
|
Line 162
|
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO ) |
SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.7.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2011 |
* June 2017 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOB |
CHARACTER JOB |
Line 189
|
Line 191
|
INTEGER I, ICA, IEXC, IRA, J, K, L, M |
INTEGER I, ICA, IEXC, IRA, J, K, L, M |
DOUBLE PRECISION C, CA, F, G, R, RA, S, SFMAX1, SFMAX2, SFMIN1, |
DOUBLE PRECISION C, CA, F, G, R, RA, S, SFMAX1, SFMAX2, SFMIN1, |
$ SFMIN2 |
$ SFMIN2 |
COMPLEX*16 CDUM |
|
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL DISNAN, LSAME |
LOGICAL DISNAN, LSAME |
INTEGER IZAMAX |
INTEGER IZAMAX |
DOUBLE PRECISION DLAMCH |
DOUBLE PRECISION DLAMCH, DZNRM2 |
EXTERNAL DISNAN, LSAME, IZAMAX, DLAMCH |
EXTERNAL DISNAN, LSAME, IZAMAX, DLAMCH, DZNRM2 |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL XERBLA, ZDSCAL, ZSWAP |
EXTERNAL XERBLA, ZDSCAL, ZSWAP |
* .. |
* .. |
* .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
INTRINSIC ABS, DBLE, DIMAG, MAX, MIN |
INTRINSIC ABS, DBLE, DIMAG, MAX, MIN |
* .. |
|
* .. Statement Functions .. |
|
DOUBLE PRECISION CABS1 |
|
* .. |
|
* .. Statement Function definitions .. |
|
CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( DIMAG( CDUM ) ) |
|
* .. |
|
* .. Executable Statements .. |
|
* |
* |
* Test the input parameters |
* Test the input parameters |
* |
* |
Line 324
|
Line 317
|
NOCONV = .FALSE. |
NOCONV = .FALSE. |
* |
* |
DO 200 I = K, L |
DO 200 I = K, L |
C = ZERO |
|
R = ZERO |
|
* |
* |
DO 150 J = K, L |
C = DZNRM2( L-K+1, A( K, I ), 1 ) |
IF( J.EQ.I ) |
R = DZNRM2( L-K+1, A( I, K ), LDA ) |
$ GO TO 150 |
|
C = C + CABS1( A( J, I ) ) |
|
R = R + CABS1( A( I, J ) ) |
|
150 CONTINUE |
|
ICA = IZAMAX( L, A( 1, I ), 1 ) |
ICA = IZAMAX( L, A( 1, I ), 1 ) |
CA = ABS( A( ICA, I ) ) |
CA = ABS( A( ICA, I ) ) |
IRA = IZAMAX( N-K+1, A( I, K ), LDA ) |
IRA = IZAMAX( N-K+1, A( I, K ), LDA ) |