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
|
PROGRAM PSRPTGSEPTST
*
* -- ScaLAPACK routine (version 1.7) --
* University of Tennessee, Knoxville, Oak Ridge National Laboratory,
* and University of California, Berkeley.
* August 14, 2001
*
* Repeat generalized parallel symmetric 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, IBTYPE, INFO, IPOSTPAD, IPREPAD,
$ LDA, MATTYPE, N, NB, NPCOL, NPROCS, NPROW
REAL ABSTOL, THRESH
* ..
* .. Local Arrays ..
INTEGER DESCA( DLEN_ ), ICLUSTR( MAXN ), IFAIL( MAXN ),
$ ISEED( 4 ), IWORK( LIWORK )
REAL A( MAXN*MAXN ), B( MAXN, MAXN ),
$ COPYA( MAXN*MAXN ), COPYB( MAXN, MAXN ),
$ GAP( MAXN ), WIN( MAXN ), WNEW( MAXN ),
$ WORK( LWORK ), Z( MAXN*MAXN )
* ..
*
*
* .. External Subroutines ..
*
EXTERNAL BLACS_EXIT, BLACS_GET, BLACS_GRIDINIT,
$ BLACS_PINFO, BLACS_SETUP, DESCINIT, PSGSEPTST
* ..
* .. Executable Statements ..
*
IPREPAD = 3
IPOSTPAD = 3
LDA = MAXN
*
* These lines should be replaced by the output from pxGSEPdriver
*
*
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
IBTYPE = 1
* note: the printout often makes a mess of ABSTOL
ABSTOL = 0.1175494351E-37
THRESH = .350000E+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 PSGSEPTST( DESCA, UPLO, N, MATTYPE, IBTYPE, SUBTESTS, THRESH,
$ N, ABSTOL, ISEED, A, COPYA, B, COPYB, Z, LDA, WIN,
$ WNEW, IFAIL, ICLUSTR, GAP, IPREPAD, IPOSTPAD,
$ WORK, 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 PSRPTGSEPTST
*
END
|
