1*> \brief \b DERRQL
2*
3*  =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6*            http://www.netlib.org/lapack/explore-html/
7*
8*  Definition:
9*  ===========
10*
11*       SUBROUTINE DERRQL( PATH, NUNIT )
12*
13*       .. Scalar Arguments ..
14*       CHARACTER*3        PATH
15*       INTEGER            NUNIT
16*       ..
17*
18*
19*> \par Purpose:
20*  =============
21*>
22*> \verbatim
23*>
24*> DERRQL tests the error exits for the DOUBLE PRECISION routines
25*> that use the QL decomposition of a general matrix.
26*> \endverbatim
27*
28*  Arguments:
29*  ==========
30*
31*> \param[in] PATH
32*> \verbatim
33*>          PATH is CHARACTER*3
34*>          The LAPACK path name for the routines to be tested.
35*> \endverbatim
36*>
37*> \param[in] NUNIT
38*> \verbatim
39*>          NUNIT is INTEGER
40*>          The unit number for output.
41*> \endverbatim
42*
43*  Authors:
44*  ========
45*
46*> \author Univ. of Tennessee
47*> \author Univ. of California Berkeley
48*> \author Univ. of Colorado Denver
49*> \author NAG Ltd.
50*
51*> \ingroup double_lin
52*
53*  =====================================================================
54      SUBROUTINE DERRQL( PATH, NUNIT )
55*
56*  -- LAPACK test routine --
57*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
58*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
59*
60*     .. Scalar Arguments ..
61      CHARACTER*3        PATH
62      INTEGER            NUNIT
63*     ..
64*
65*  =====================================================================
66*
67*     .. Parameters ..
68      INTEGER            NMAX
69      PARAMETER          ( NMAX = 2 )
70*     ..
71*     .. Local Scalars ..
72      INTEGER            I, INFO, J
73*     ..
74*     .. Local Arrays ..
75      DOUBLE PRECISION   A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
76     $                   W( NMAX ), X( NMAX )
77*     ..
78*     .. External Subroutines ..
79      EXTERNAL           ALAESM, CHKXER, DGEQL2, DGEQLF, DGEQLS, DORG2L,
80     $                   DORGQL, DORM2L, DORMQL
81*     ..
82*     .. Scalars in Common ..
83      LOGICAL            LERR, OK
84      CHARACTER*32       SRNAMT
85      INTEGER            INFOT, NOUT
86*     ..
87*     .. Common blocks ..
88      COMMON             / INFOC / INFOT, NOUT, OK, LERR
89      COMMON             / SRNAMC / SRNAMT
90*     ..
91*     .. Intrinsic Functions ..
92      INTRINSIC          DBLE
93*     ..
94*     .. Executable Statements ..
95*
96      NOUT = NUNIT
97      WRITE( NOUT, FMT = * )
98*
99*     Set the variables to innocuous values.
100*
101      DO 20 J = 1, NMAX
102         DO 10 I = 1, NMAX
103            A( I, J ) = 1.D0 / DBLE( I+J )
104            AF( I, J ) = 1.D0 / DBLE( I+J )
105   10    CONTINUE
106         B( J ) = 0.D0
107         W( J ) = 0.D0
108         X( J ) = 0.D0
109   20 CONTINUE
110      OK = .TRUE.
111*
112*     Error exits for QL factorization
113*
114*     DGEQLF
115*
116      SRNAMT = 'DGEQLF'
117      INFOT = 1
118      CALL DGEQLF( -1, 0, A, 1, B, W, 1, INFO )
119      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
120      INFOT = 2
121      CALL DGEQLF( 0, -1, A, 1, B, W, 1, INFO )
122      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
123      INFOT = 4
124      CALL DGEQLF( 2, 1, A, 1, B, W, 1, INFO )
125      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
126      INFOT = 7
127      CALL DGEQLF( 1, 2, A, 1, B, W, 1, INFO )
128      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
129*
130*     DGEQL2
131*
132      SRNAMT = 'DGEQL2'
133      INFOT = 1
134      CALL DGEQL2( -1, 0, A, 1, B, W, INFO )
135      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
136      INFOT = 2
137      CALL DGEQL2( 0, -1, A, 1, B, W, INFO )
138      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
139      INFOT = 4
140      CALL DGEQL2( 2, 1, A, 1, B, W, INFO )
141      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
142*
143*     DGEQLS
144*
145      SRNAMT = 'DGEQLS'
146      INFOT = 1
147      CALL DGEQLS( -1, 0, 0, A, 1, X, B, 1, W, 1, INFO )
148      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
149      INFOT = 2
150      CALL DGEQLS( 0, -1, 0, A, 1, X, B, 1, W, 1, INFO )
151      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
152      INFOT = 2
153      CALL DGEQLS( 1, 2, 0, A, 1, X, B, 1, W, 1, INFO )
154      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
155      INFOT = 3
156      CALL DGEQLS( 0, 0, -1, A, 1, X, B, 1, W, 1, INFO )
157      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
158      INFOT = 5
159      CALL DGEQLS( 2, 1, 0, A, 1, X, B, 2, W, 1, INFO )
160      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
161      INFOT = 8
162      CALL DGEQLS( 2, 1, 0, A, 2, X, B, 1, W, 1, INFO )
163      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
164      INFOT = 10
165      CALL DGEQLS( 1, 1, 2, A, 1, X, B, 1, W, 1, INFO )
166      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
167*
168*     DORGQL
169*
170      SRNAMT = 'DORGQL'
171      INFOT = 1
172      CALL DORGQL( -1, 0, 0, A, 1, X, W, 1, INFO )
173      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
174      INFOT = 2
175      CALL DORGQL( 0, -1, 0, A, 1, X, W, 1, INFO )
176      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
177      INFOT = 2
178      CALL DORGQL( 1, 2, 0, A, 1, X, W, 2, INFO )
179      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
180      INFOT = 3
181      CALL DORGQL( 0, 0, -1, A, 1, X, W, 1, INFO )
182      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
183      INFOT = 3
184      CALL DORGQL( 1, 1, 2, A, 1, X, W, 1, INFO )
185      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
186      INFOT = 5
187      CALL DORGQL( 2, 1, 0, A, 1, X, W, 1, INFO )
188      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
189      INFOT = 8
190      CALL DORGQL( 2, 2, 0, A, 2, X, W, 1, INFO )
191      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
192*
193*     DORG2L
194*
195      SRNAMT = 'DORG2L'
196      INFOT = 1
197      CALL DORG2L( -1, 0, 0, A, 1, X, W, INFO )
198      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
199      INFOT = 2
200      CALL DORG2L( 0, -1, 0, A, 1, X, W, INFO )
201      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
202      INFOT = 2
203      CALL DORG2L( 1, 2, 0, A, 1, X, W, INFO )
204      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
205      INFOT = 3
206      CALL DORG2L( 0, 0, -1, A, 1, X, W, INFO )
207      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
208      INFOT = 3
209      CALL DORG2L( 2, 1, 2, A, 2, X, W, INFO )
210      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
211      INFOT = 5
212      CALL DORG2L( 2, 1, 0, A, 1, X, W, INFO )
213      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
214*
215*     DORMQL
216*
217      SRNAMT = 'DORMQL'
218      INFOT = 1
219      CALL DORMQL( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
220      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
221      INFOT = 2
222      CALL DORMQL( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
223      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
224      INFOT = 3
225      CALL DORMQL( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
226      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
227      INFOT = 4
228      CALL DORMQL( 'L', 'N', 0, -1, 0, A, 1, X, AF, 1, W, 1, INFO )
229      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
230      INFOT = 5
231      CALL DORMQL( 'L', 'N', 0, 0, -1, A, 1, X, AF, 1, W, 1, INFO )
232      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
233      INFOT = 5
234      CALL DORMQL( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, 1, INFO )
235      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
236      INFOT = 5
237      CALL DORMQL( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, 1, INFO )
238      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
239      INFOT = 7
240      CALL DORMQL( 'L', 'N', 2, 1, 0, A, 1, X, AF, 2, W, 1, INFO )
241      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
242      INFOT = 7
243      CALL DORMQL( 'R', 'N', 1, 2, 0, A, 1, X, AF, 1, W, 1, INFO )
244      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
245      INFOT = 10
246      CALL DORMQL( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, 1, INFO )
247      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
248      INFOT = 12
249      CALL DORMQL( 'L', 'N', 1, 2, 0, A, 1, X, AF, 1, W, 1, INFO )
250      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
251      INFOT = 12
252      CALL DORMQL( 'R', 'N', 2, 1, 0, A, 1, X, AF, 2, W, 1, INFO )
253      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
254*
255*     DORM2L
256*
257      SRNAMT = 'DORM2L'
258      INFOT = 1
259      CALL DORM2L( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
260      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
261      INFOT = 2
262      CALL DORM2L( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
263      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
264      INFOT = 3
265      CALL DORM2L( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, INFO )
266      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
267      INFOT = 4
268      CALL DORM2L( 'L', 'N', 0, -1, 0, A, 1, X, AF, 1, W, INFO )
269      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
270      INFOT = 5
271      CALL DORM2L( 'L', 'N', 0, 0, -1, A, 1, X, AF, 1, W, INFO )
272      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
273      INFOT = 5
274      CALL DORM2L( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, INFO )
275      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
276      INFOT = 5
277      CALL DORM2L( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, INFO )
278      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
279      INFOT = 7
280      CALL DORM2L( 'L', 'N', 2, 1, 0, A, 1, X, AF, 2, W, INFO )
281      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
282      INFOT = 7
283      CALL DORM2L( 'R', 'N', 1, 2, 0, A, 1, X, AF, 1, W, INFO )
284      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
285      INFOT = 10
286      CALL DORM2L( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, INFO )
287      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
288*
289*     Print a summary line.
290*
291      CALL ALAESM( PATH, OK, NOUT )
292*
293      RETURN
294*
295*     End of DERRQL
296*
297      END
298