ScaLAPACK 2.1  2.1
ScaLAPACK: Scalable Linear Algebra PACKage
pdrptseptst.f
Go to the documentation of this file.
1 *
2 *
3  PROGRAM pdrptseptst
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 symmetric 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 hetero, 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 a( maxn*maxn ), copya( maxn*maxn ),
29  $ gap( maxn ), win( maxn ), wnew( maxn ),
30  $ work( lwork ), z( maxn*maxn )
31 * ..
32 *
33 *
34 * .. External Subroutines ..
35 *
36  EXTERNAL blacs_exit, blacs_get, blacs_gridinit,
37  $ blacs_pinfo, blacs_setup, descinit, pdseptst
38 * ..
39 * .. Executable Statements ..
40 *
41  iprepad = 3
42  ipostpad = 3
43  lda = maxn
44 *
45 * Set HETERO to 'Y' if you want to turn off the PxSYEV tests
46 *
47  hetero = 'N'
48 *
49 * These lines should be replaced by the output from pxSEPdriver
50 *
51 *
52  iseed( 1 ) = 2312
53  iseed( 2 ) = 3709
54  iseed( 3 ) = 666
55  iseed( 4 ) = 3371
56  uplo = 'U'
57  subtests = 'Y'
58  n = 33
59  nprow = 2
60  npcol = 2
61  nb = 4
62  mattype = 9
63 * note: the printout often makes a mess of ABSTOL
64  abstol = 0.1175494351d-37
65  thresh = .350000d+01
66 *
67  CALL blacs_pinfo( iam, nprocs )
68  IF( nprocs.LT.1 ) THEN
69 *
70  nprocs = nprow*npcol
71  CALL blacs_setup( iam, nprocs )
72  END IF
73  CALL blacs_get( -1, 0, context )
74  CALL blacs_gridinit( context, 'R', nprow, npcol )
75 *
76  CALL descinit( desca, n, n, nb, nb, 0, 0, context, lda, info )
77 *
78  CALL pdseptst( desca, uplo, n, mattype, subtests, thresh, n,
79  $ abstol, iseed, a, copya, z, lda, win, wnew, ifail,
80  $ iclustr, gap, iprepad, ipostpad, work,
81  $ lwork-iprepad-ipostpad, iwork,
82  $ liwork-iprepad-ipostpad, hetero, 6, info )
83 *
84 *
85 *
86 * Uncomment this line on SUN systems to avoid the useless print out
87 *
88 * CALL IEEE_FLAGS( 'clear', 'exception', 'underflow', '')
89 *
90 *
91 *
92 *
93  CALL blacs_exit( 0 )
94  stop
95 *
96 *
97 *
98 * End of PDRPTSEPTST
99 *
100  END
pdrptseptst
program pdrptseptst
Definition: pdrptseptst.f:3
pdseptst
subroutine pdseptst(DESCA, UPLO, N, MATTYPE, SUBTESTS, THRESH, ORDER, ABSTOL, ISEED, A, COPYA, Z, LDA, WIN, WNEW, IFAIL, ICLUSTR, GAP, IPREPAD, IPOSTPAD, WORK, LWORK, IWORK, LIWORK, HETERO, NOUT, INFO)
Definition: pdseptst.f:6
descinit
subroutine descinit(DESC, M, N, MB, NB, IRSRC, ICSRC, ICTXT, LLD, INFO)
Definition: descinit.f:3