1 SUBROUTINE cqrt13( SCALE, M, N, A, LDA, NORMA, ISEED )
8 INTEGER lda, m, n, scale
52 parameter( one = 1.0e0 )
59 REAL clange, scasum, slamch
60 EXTERNAL clange, scasum, slamch
66 INTRINSIC cmplx,
REAL, sign
73 IF( m.LE.0 .OR. n.LE.0 )
79 CALL clarnv( 2, iseed, m, a( 1, j ) )
81 a( j, j ) = a( j, j ) + cmplx( sign( scasum( m, a( 1, j ),
82 $ 1 ),
REAL( A( J, J ) ) ) )
89 norma = clange(
'Max', m, n, a, lda, dummy )
90 smlnum = slamch(
'Safe minimum' )
92 CALL
slabad( smlnum, bignum )
93 smlnum = smlnum / slamch(
'Epsilon' )
100 CALL
clascl(
'General', 0, 0, norma, bignum, m, n, a, lda,
102 ELSE IF( scale.EQ.3 )
THEN
106 CALL
clascl(
'General', 0, 0, norma, smlnum, m, n, a, lda,
111 norma = clange(
'One-norm', m, n, a, lda, dummy )