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
|
LOGICAL FUNCTION SLCTSX( AR, AI, BETA )
*
* -- LAPACK test routine (version 3.0) --
* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
* Courant Institute, Argonne National Lab, and Rice University
* June 30, 1999
*
* .. Scalar Arguments ..
REAL AI, AR, BETA
* ..
*
* Purpose
* =======
*
* This function is used to determine what eigenvalues will be
* selected. If this is part of the test driver SDRGSX, do not
* change the code UNLESS you are testing input examples and not
* using the built-in examples.
*
* Arguments
* =========
*
* AR (input) REAL
* The numerator of the real part of a complex eigenvalue
* (AR/BETA) + i*(AI/BETA).
*
* AI (input) REAL
* The numerator of the imaginary part of a complex eigenvalue
* (AR/BETA) + i*(AI).
*
* BETA (input) REAL
* The denominator part of a complex eigenvalue
* (AR/BETA) + i*(AI/BETA).
*
* =====================================================================
*
* .. Scalars in Common ..
LOGICAL FS
INTEGER I, M, MPLUSN, N
* ..
* .. Common blocks ..
COMMON / MN / M, N, MPLUSN, I, FS
* ..
* .. Save statement ..
SAVE
* ..
* .. Executable Statements ..
*
IF( FS ) THEN
I = I + 1
IF( I.LE.M ) THEN
SLCTSX = .FALSE.
ELSE
SLCTSX = .TRUE.
END IF
IF( I.EQ.MPLUSN ) THEN
FS = .FALSE.
I = 0
END IF
ELSE
I = I + 1
IF( I.LE.N ) THEN
SLCTSX = .TRUE.
ELSE
SLCTSX = .FALSE.
END IF
IF( I.EQ.MPLUSN ) THEN
FS = .TRUE.
I = 0
END IF
END IF
*
* IF( AR/BETA.GT.0.0 )THEN
* SLCTSX = .TRUE.
* ELSE
* SLCTSX = .FALSE.
* END IF
*
RETURN
*
* End of SLCTSX
*
END
|
