151 SUBROUTINE zdrvab( DOTYPE, NM, MVAL, NNS,
152 $ nsval, thresh, nmax, a, afac, b,
153 $ x, work, rwork, swork, iwork, nout )
161 INTEGER NM, NMAX, NNS, NOUT
162 DOUBLE PRECISION THRESH
166 INTEGER MVAL( * ), NSVAL( * ), IWORK( * )
167 DOUBLE PRECISION RWORK( * )
169 COMPLEX*16 A( * ), AFAC( * ), B( * ),
176 DOUBLE PRECISION ZERO
177 parameter ( zero = 0.0d+0 )
179 parameter ( ntypes = 11 )
181 parameter ( ntests = 1 )
185 CHARACTER DIST, TRANS,
TYPE, XTYPE
187 INTEGER I, IM, IMAT, INFO, IOFF, IRHS,
188 $ izero, kl, ku, lda, m, mode, n,
189 $ nerrs, nfail, nimat, nrhs, nrun
190 DOUBLE PRECISION ANORM, CNDNUM
193 INTEGER ISEED( 4 ), ISEEDY( 4 )
194 DOUBLE PRECISION RESULT( ntests )
204 INTRINSIC dcmplx, dble, max, min, sqrt
212 COMMON / infoc / infot, nunit, ok, lerr
213 COMMON / srnamc / srnamt
216 DATA iseedy / 2006, 2007, 2008, 2009 /
223 path( 1: 1 ) =
'Zomplex precision'
229 iseed( i ) = iseedy( i )
242 IF( m.LE.0 .OR. n.LE.0 )
245 DO 100 imat = 1, nimat
249 IF( .NOT.dotype( imat ) )
254 zerot = imat.GE.5 .AND. imat.LE.7
255 IF( zerot .AND. n.LT.imat-4 )
261 CALL zlatb4( path, imat, m, n,
TYPE, KL, KU, ANORM, MODE,
265 CALL zlatms( m, n, dist, iseed,
TYPE, RWORK, MODE,
266 $ cndnum, anorm, kl, ku,
'No packing', a, lda,
272 CALL alaerh( path,
'ZLATMS', info, 0,
' ', m, n, -1,
273 $ -1, -1, imat, nfail, nerrs, nout )
283 ELSE IF( imat.EQ.6 )
THEN
286 izero = min( m, n ) / 2 + 1
288 ioff = ( izero-1 )*lda
294 CALL zlaset(
'Full', m, n-izero+1, dcmplx(zero),
295 $ dcmplx(zero), a( ioff+1 ), lda )
307 CALL zlarhs( path, xtype,
' ', trans, n, n, kl,
308 $ ku, nrhs, a, lda, x, lda, b,
315 CALL zlacpy(
'Full', m, n, a, lda, afac, lda )
317 CALL zcgesv( n, nrhs, a, lda, iwork, b, lda, x, lda,
318 $ work, swork, rwork, iter, info)
321 CALL zlacpy(
'Full', m, n, afac, lda, a, lda )
327 IF( info.NE.izero )
THEN
329 IF( nfail.EQ.0 .AND. nerrs.EQ.0 )
330 $
CALL alahd( nout, path )
333 IF( info.NE.izero .AND. izero.NE.0 )
THEN
334 WRITE( nout, fmt = 9988 )
'ZCGESV',info,
337 WRITE( nout, fmt = 9975 )
'ZCGESV',info,
349 CALL zlacpy(
'Full', n, nrhs, b, lda, work, lda )
351 CALL zget08( trans, n, n, nrhs, a, lda, x, lda, work,
352 $ lda, rwork, result( 1 ) )
366 IF ((thresh.LE.0.0e+00)
367 $ .OR.((iter.GE.0).AND.(n.GT.0)
368 $ .AND.(result(1).GE.sqrt(dble(n))))
369 $ .OR.((iter.LT.0).AND.(result(1).GE.thresh)))
THEN
371 IF( nfail.EQ.0 .AND. nerrs.EQ.0 )
THEN
372 WRITE( nout, fmt = 8999 )
'DGE'
373 WRITE( nout, fmt =
'( '' Matrix types:'' )' )
374 WRITE( nout, fmt = 8979 )
375 WRITE( nout, fmt =
'( '' Test ratios:'' )' )
376 WRITE( nout, fmt = 8960 )1
377 WRITE( nout, fmt =
'( '' Messages:'' )' )
380 WRITE( nout, fmt = 9998 )trans, n, nrhs,
381 $ imat, 1, result( 1 )
391 IF( nfail.GT.0 )
THEN
392 WRITE( nout, fmt = 9996 )
'ZCGESV', nfail, nrun
394 WRITE( nout, fmt = 9995 )
'ZCGESV', nrun
396 IF( nerrs.GT.0 )
THEN
397 WRITE( nout, fmt = 9994 )nerrs
400 9998
FORMAT(
' TRANS=''', a1,
''', N =', i5,
', NRHS=', i3,
', type ',
401 $ i2,
', test(', i2,
') =', g12.5 )
402 9996
FORMAT( 1x, a6,
': ', i6,
' out of ', i6,
403 $
' tests failed to pass the threshold' )
404 9995
FORMAT( /1x,
'All tests for ', a6,
405 $
' routines passed the threshold ( ', i6,
' tests run)' )
406 9994
FORMAT( 6x, i6,
' error messages recorded' )
410 9988
FORMAT(
' *** ', a6,
' returned with INFO =', i5,
' instead of ',
411 $ i5, /
' ==> M =', i5,
', type ',
416 9975
FORMAT(
' *** Error code from ', a6,
'=', i5,
' for M=', i5,
418 8999
FORMAT( / 1x, a3,
': General dense matrices' )
419 8979
FORMAT( 4x,
'1. Diagonal', 24x,
'7. Last n/2 columns zero', / 4x,
420 $
'2. Upper triangular', 16x,
421 $
'8. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
422 $
'3. Lower triangular', 16x,
'9. Random, CNDNUM = 0.1/EPS',
423 $ / 4x,
'4. Random, CNDNUM = 2', 13x,
424 $
'10. Scaled near underflow', / 4x,
'5. First column zero',
425 $ 14x,
'11. Scaled near overflow', / 4x,
426 $
'6. Last column zero' )
427 8960
FORMAT( 3x, i2,
': norm_1( B - A * X ) / ',
428 $
'( norm_1(A) * norm_1(X) * EPS * SQRT(N) ) > 1 if ITERREF',
429 $ / 4x,
'or norm_1( B - A * X ) / ',
430 $
'( norm_1(A) * norm_1(X) * EPS ) > THRES if DGETRF' )
subroutine alahd(IOUNIT, PATH)
ALAHD
subroutine alaerh(PATH, SUBNAM, INFO, INFOE, OPTS, M, N, KL, KU, N5, IMAT, NFAIL, NERRS, NOUT)
ALAERH
subroutine zlarhs(PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B, LDB, ISEED, INFO)
ZLARHS
subroutine zlacpy(UPLO, M, N, A, LDA, B, LDB)
ZLACPY copies all or part of one two-dimensional array to another.
subroutine zlatb4(PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE, CNDNUM, DIST)
ZLATB4
subroutine zdrvab(DOTYPE, NM, MVAL, NNS, NSVAL, THRESH, NMAX, A, AFAC, B, X, WORK, RWORK, SWORK, IWORK, NOUT)
ZDRVAB
subroutine zlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
ZLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
subroutine zlatms(M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX, KL, KU, PACK, A, LDA, WORK, INFO)
ZLATMS
subroutine zget08(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
ZGET08
subroutine zcgesv(N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK, SWORK, RWORK, ITER, INFO)
ZCGESV computes the solution to system of linear equations A * X = B for GE matrices (mixed precisio...