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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
---
:name: slasq3
:md5sum: 2f410ac97a48fb4a07e823012fc4b2d8
:category: :subroutine
:arguments:
- i0:
:type: integer
:intent: input
- n0:
:type: integer
:intent: input/output
- z:
:type: real
:intent: input
:dims:
- 4*n0
- pp:
:type: integer
:intent: input/output
- dmin:
:type: real
:intent: output
- sigma:
:type: real
:intent: output
- desig:
:type: real
:intent: input/output
- qmax:
:type: real
:intent: input
- nfail:
:type: integer
:intent: output
- iter:
:type: integer
:intent: output
- ndiv:
:type: integer
:intent: output
- ieee:
:type: logical
:intent: input
- ttype:
:type: integer
:intent: input/output
- dmin1:
:type: real
:intent: input/output
- dmin2:
:type: real
:intent: input/output
- dn:
:type: real
:intent: input/output
- dn1:
:type: real
:intent: input/output
- dn2:
:type: real
:intent: input/output
- g:
:type: real
:intent: input/output
- tau:
:type: real
:intent: input/output
:substitutions: {}
:fortran_help: " SUBROUTINE SLASQ3( I0, N0, Z, PP, DMIN, SIGMA, DESIG, QMAX, NFAIL, ITER, NDIV, IEEE, TTYPE, DMIN1, DMIN2, DN, DN1, DN2, G, TAU )\n\n\
* Purpose\n\
* =======\n\
*\n\
* SLASQ3 checks for deflation, computes a shift (TAU) and calls dqds.\n\
* In case of failure it changes shifts, and tries again until output\n\
* is positive.\n\
*\n\n\
* Arguments\n\
* =========\n\
*\n\
* I0 (input) INTEGER\n\
* First index.\n\
*\n\
* N0 (input/output) INTEGER\n\
* Last index.\n\
*\n\
* Z (input) REAL array, dimension ( 4*N )\n\
* Z holds the qd array.\n\
*\n\
* PP (input/output) INTEGER\n\
* PP=0 for ping, PP=1 for pong.\n\
* PP=2 indicates that flipping was applied to the Z array \n\
* and that the initial tests for deflation should not be \n\
* performed.\n\
*\n\
* DMIN (output) REAL\n\
* Minimum value of d.\n\
*\n\
* SIGMA (output) REAL\n\
* Sum of shifts used in current segment.\n\
*\n\
* DESIG (input/output) REAL\n\
* Lower order part of SIGMA\n\
*\n\
* QMAX (input) REAL\n\
* Maximum value of q.\n\
*\n\
* NFAIL (output) INTEGER\n\
* Number of times shift was too big.\n\
*\n\
* ITER (output) INTEGER\n\
* Number of iterations.\n\
*\n\
* NDIV (output) INTEGER\n\
* Number of divisions.\n\
*\n\
* IEEE (input) LOGICAL\n\
* Flag for IEEE or non IEEE arithmetic (passed to SLASQ5).\n\
*\n\
* TTYPE (input/output) INTEGER\n\
* Shift type.\n\
*\n\
* DMIN1 (input/output) REAL\n\
*\n\
* DMIN2 (input/output) REAL\n\
*\n\
* DN (input/output) REAL\n\
*\n\
* DN1 (input/output) REAL\n\
*\n\
* DN2 (input/output) REAL\n\
*\n\
* G (input/output) REAL\n\
*\n\
* TAU (input/output) REAL\n\
*\n\
* These are passed as arguments in order to save their values\n\
* between calls to SLASQ3.\n\
*\n\n\
* =====================================================================\n\
*\n"
|
