1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
|
*
*
PROGRAM PZRPTSEPTST
*
* -- ScaLAPACK routine (version 1.7) --
* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
* and University of California, Berkeley.
* November 15, 1997
*
* Repeat parallel Hermitian eigenproblem test
* .. Parameters ..
*
INTEGER DLEN_
PARAMETER ( DLEN_ = 9 )
INTEGER MAXN, LWORK, LIWORK
PARAMETER ( MAXN = 200, LWORK = 500000,
$ LIWORK = 6*MAXN+4 )
* ..
* .. Local Scalars ..
CHARACTER SUBTESTS, UPLO
INTEGER CONTEXT, IAM, INFO, IPOSTPAD, IPREPAD, LDA,
$ MATTYPE, N, NB, NPCOL, NPROCS, NPROW
DOUBLE PRECISION ABSTOL, THRESH
* ..
* .. Local Arrays ..
INTEGER DESCA( DLEN_ ), ICLUSTR( MAXN ), IFAIL( MAXN ),
$ ISEED( 4 ), IWORK( LIWORK )
DOUBLE PRECISION GAP( MAXN ), RWORK( LWORK ), WIN( MAXN ),
$ WNEW( MAXN )
COMPLEX*16 A( MAXN*MAXN ), COPYA( MAXN*MAXN ),
$ WORK( LWORK ), Z( MAXN*MAXN )
* ..
*
*
* .. External Subroutines ..
*
EXTERNAL BLACS_EXIT, BLACS_GET, BLACS_GRIDINIT,
$ BLACS_PINFO, BLACS_SETUP, DESCINIT, PZSEPTST
* ..
* .. Executable Statements ..
*
IPREPAD = 3
IPOSTPAD = 3
LDA = MAXN
*
* These lines should be replaced by the output from pxSEPdriver
*
*
ISEED( 1 ) = 2312
ISEED( 2 ) = 3709
ISEED( 3 ) = 666
ISEED( 4 ) = 3371
UPLO = 'U'
SUBTESTS = 'Y'
N = 33
NPROW = 2
NPCOL = 2
NB = 4
MATTYPE = 9
* note: the printout often makes a mess of ABSTOL
ABSTOL = 0.1175494351D-37
THRESH = .350000D+01
*
CALL BLACS_PINFO( IAM, NPROCS )
IF( NPROCS.LT.1 ) THEN
*
NPROCS = NPROW*NPCOL
CALL BLACS_SETUP( IAM, NPROCS )
END IF
CALL BLACS_GET( -1, 0, CONTEXT )
CALL BLACS_GRIDINIT( CONTEXT, 'R', NPROW, NPCOL )
*
CALL DESCINIT( DESCA, N, N, NB, NB, 0, 0, CONTEXT, LDA, INFO )
*
CALL PZSEPTST( DESCA, UPLO, N, MATTYPE, SUBTESTS, THRESH, N,
$ ABSTOL, ISEED, A, COPYA, Z, LDA, WIN, WNEW, IFAIL,
$ ICLUSTR, GAP, IPREPAD, IPOSTPAD, WORK,
$ LWORK-IPREPAD-IPOSTPAD, RWORK,
$ LWORK-IPREPAD-IPOSTPAD, IWORK,
$ LIWORK-IPREPAD-IPOSTPAD, 6, INFO )
*
*
*
* Uncomment this line on SUN systems to avoid the useless print out
*
* CALL IEEE_FLAGS( 'clear', 'exception', 'underflow', '')
*
*
*
*
CALL BLACS_EXIT( 0 )
STOP
*
*
*
* End of PZRPTSEPTST
*
END
|
