ScaLAPACK 2.1  2.1 ScaLAPACK: Scalable Linear Algebra PACKage
pcrptseptst.f
Go to the documentation of this file.
1  PROGRAM pcrptseptst
2 *
3 * -- ScaLAPACK routine (version 1.7) --
4 * University of Tennessee, Knoxville, Oak Ridge National Laboratory,
5 * and University of California, Berkeley.
6 * August 14, 2001
7 *
8 * Repeat parallel Hermitian eigenproblem test
9 * .. Parameters ..
10 *
11  INTEGER dlen_
12  parameter( dlen_ = 9 )
13  INTEGER maxn, lwork, liwork
14  parameter( maxn = 200, lwork = 500000,
15  \$ liwork = 6*maxn+4 )
16 * ..
17 * .. Local Scalars ..
18  CHARACTER subtests, uplo
20  \$ mattype, n, nb, npcol, nprocs, nprow
21  REAL abstol, thresh
22 * ..
23 * .. Local Arrays ..
24  INTEGER desca( dlen_ ), iclustr( maxn ), ifail( maxn ),
25  \$ iseed( 4 ), iwork( liwork )
26  REAL gap( maxn ), rwork( lwork ), win( maxn ),
27  \$ wnew( maxn )
28  COMPLEX a( maxn*maxn ), copya( maxn*maxn ),
29  \$ work( lwork ), z( maxn*maxn )
30 * ..
31 *
32 *
33 * .. External Subroutines ..
34 *
35  EXTERNAL blacs_exit, blacs_get, blacs_gridinit,
36  \$ blacs_pinfo, blacs_setup, descinit, pcseptst
37 * ..
38 * .. Executable Statements ..
39 *
42  lda = maxn
43 *
44 * These lines should be replaced by the output from pxSEPdriver
45 *
46 *
47  iseed( 1 ) = 2312
48  iseed( 2 ) = 3709
49  iseed( 3 ) = 666
50  iseed( 4 ) = 3371
51  uplo = 'U'
52  subtests = 'Y'
53  n = 33
54  nprow = 2
55  npcol = 2
56  nb = 4
57  mattype = 9
58 * note: the printout often makes a mess of ABSTOL
59  abstol = 0.1175494351e-37
60  thresh = .350000e+01
61 *
62  CALL blacs_pinfo( iam, nprocs )
63  IF( nprocs.LT.1 ) THEN
64 *
65  nprocs = nprow*npcol
66  CALL blacs_setup( iam, nprocs )
67  END IF
68  CALL blacs_get( -1, 0, context )
69  CALL blacs_gridinit( context, 'R', nprow, npcol )
70 *
71  CALL descinit( desca, n, n, nb, nb, 0, 0, context, lda, info )
72 *
73  CALL pcseptst( desca, uplo, n, mattype, subtests, thresh, n,
74  \$ abstol, iseed, a, copya, z, lda, win, wnew, ifail,
79 *
80 *
81 *
82 * Uncomment this line on SUN systems to avoid the useless print out
83 *
84 * CALL IEEE_FLAGS( 'clear', 'exception', 'underflow', '')
85 *
86 *
87 *
88 *
89  CALL blacs_exit( 0 )
90  stop
91 *
92 *
93 *
94 * End of PCRPTSEPTST
95 *
96  END
pcseptst
subroutine pcseptst(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: pcseptst.f:8
pcrptseptst
program pcrptseptst
Definition: pcrptseptst.f:1
descinit
subroutine descinit(DESC, M, N, MB, NB, IRSRC, ICSRC, ICTXT, LLD, INFO)
Definition: descinit.f:3