SCALAPACK 2.2.2
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches
pzrptseptst.f
Go to the documentation of this file.
1*
2*
3 PROGRAM pzrptseptst
4*
5* -- ScaLAPACK routine (version 1.7) --
6* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
7* and University of California, Berkeley.
8* November 15, 1997
9*
10* Repeat parallel Hermitian eigenproblem test
11* .. Parameters ..
12*
13 INTEGER dlen_
14 parameter( dlen_ = 9 )
15 INTEGER maxn, lwork, liwork
16 parameter( maxn = 200, lwork = 500000,
17 $ liwork = 6*maxn+4 )
18* ..
19* .. Local Scalars ..
20 CHARACTER subtests, uplo
21 INTEGER context, iam, info, ipostpad, iprepad, lda,
22 $ mattype, n, nb, npcol, nprocs, nprow
23 DOUBLE PRECISION abstol, thresh
24* ..
25* .. Local Arrays ..
26 INTEGER desca( dlen_ ), iclustr( maxn ), ifail( maxn ),
27 $ iseed( 4 ), iwork( liwork )
28 DOUBLE PRECISION gap( maxn ), rwork( lwork ), win( maxn ),
29 $ wnew( maxn )
30 COMPLEX*16 a( maxn*maxn ), copya( maxn*maxn ),
31 $ work( lwork ), z( maxn*maxn )
32* ..
33*
34*
35* .. External Subroutines ..
36*
37 EXTERNAL blacs_exit, blacs_get, blacs_gridinit,
38 $ blacs_pinfo, blacs_setup, descinit, pzseptst
39* ..
40* .. Executable Statements ..
41*
42 iprepad = 3
43 ipostpad = 3
44 lda = maxn
45*
46* These lines should be replaced by the output from pxSEPdriver
47*
48*
49 iseed( 1 ) = 2312
50 iseed( 2 ) = 3709
51 iseed( 3 ) = 666
52 iseed( 4 ) = 3371
53 uplo = 'U'
54 subtests = 'Y'
55 n = 33
56 nprow = 2
57 npcol = 2
58 nb = 4
59 mattype = 9
60* note: the printout often makes a mess of ABSTOL
61 abstol = 0.1175494351d-37
62 thresh = .350000d+01
63*
64 CALL blacs_pinfo( iam, nprocs )
65 IF( nprocs.LT.1 ) THEN
66*
67 nprocs = nprow*npcol
68 CALL blacs_setup( iam, nprocs )
69 END IF
70 CALL blacs_get( -1, 0, context )
71 CALL blacs_gridinit( context, 'R', nprow, npcol )
72*
73 CALL descinit( desca, n, n, nb, nb, 0, 0, context, lda, info )
74*
75 CALL pzseptst( desca, uplo, n, mattype, subtests, thresh, n,
76 $ abstol, iseed, a, copya, z, lda, win, wnew, ifail,
77 $ iclustr, gap, iprepad, ipostpad, work,
78 $ lwork-iprepad-ipostpad, rwork,
79 $ lwork-iprepad-ipostpad, iwork,
80 $ liwork-iprepad-ipostpad, 6, info )
81*
82*
83*
84* Uncomment this line on SUN systems to avoid the useless print out
85*
86* CALL IEEE_FLAGS( 'clear', 'exception', 'underflow', '')
87*
88*
89*
90*
91 CALL blacs_exit( 0 )
92 stop
93*
94*
95*
96* End of PZRPTSEPTST
97*
98 END
subroutine descinit(desc, m, n, mb, nb, irsrc, icsrc, ictxt, lld, info)
Definition descinit.f:3
program pzrptseptst
Definition pzrptseptst.f:3
subroutine pzseptst(desca, uplo, n, mattype, subtests, thresh, order, abstol, iseed, a, copya, z, lda, win, wnew, ifail, iclustr, gap, iprepad, ipostpad, work, lwork, rwork, lrwork, iwork, liwork, nout, info)
Definition pzseptst.f:8