ScaLAPACK 2.1  2.1 ScaLAPACK: Scalable Linear Algebra PACKage
psrptseptst.f
Go to the documentation of this file.
1 *
2 *
3  PROGRAM psrptseptst
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
22  \$ mattype, n, nb, npcol, nprocs, nprow
23  REAL abstol, thresh
24 * ..
25 * .. Local Arrays ..
26  INTEGER desca( dlen_ ), iclustr( maxn ), ifail( maxn ),
27  \$ iseed( 4 ), iwork( liwork )
28  REAL 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, psseptst
38 * ..
39 * .. Executable Statements ..
40 *
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.1175494351e-37
65  thresh = .350000e+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 psseptst( desca, uplo, n, mattype, subtests, thresh, n,
79  \$ abstol, iseed, a, copya, z, lda, win, wnew, ifail,
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 PSRPTSEPTST
99 *
100  END
psseptst
subroutine psseptst(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: psseptst.f:6
descinit
subroutine descinit(DESC, M, N, MB, NB, IRSRC, ICSRC, ICTXT, LLD, INFO)
Definition: descinit.f:3
psrptseptst
program psrptseptst
Definition: psrptseptst.f:3