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*> \date December 2016
52*
53*> \ingroup double_lin
54*
55*  =====================================================================
56      SUBROUTINE DERRQL( PATH, NUNIT )
57*
58*  -- LAPACK test routine (version 3.7.0) --
59*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
60*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
61*     December 2016
62*
63*     .. Scalar Arguments ..
64      CHARACTER*3        PATH
65      INTEGER            NUNIT
66*     ..
67*
68*  =====================================================================
69*
70*     .. Parameters ..
71      INTEGER            NMAX
72      PARAMETER          ( NMAX = 2 )
73*     ..
74*     .. Local Scalars ..
75      INTEGER            I, INFO, J
76*     ..
77*     .. Local Arrays ..
78      DOUBLE PRECISION   A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
79     $                   W( NMAX ), X( NMAX )
80*     ..
81*     .. External Subroutines ..
82      EXTERNAL           ALAESM, CHKXER, DGEQL2, DGEQLF, DGEQLS, DORG2L,
83     $                   DORGQL, DORM2L, DORMQL
84*     ..
85*     .. Scalars in Common ..
86      LOGICAL            LERR, OK
87      CHARACTER*32       SRNAMT
88      INTEGER            INFOT, NOUT
89*     ..
90*     .. Common blocks ..
91      COMMON             / INFOC / INFOT, NOUT, OK, LERR
92      COMMON             / SRNAMC / SRNAMT
93*     ..
94*     .. Intrinsic Functions ..
95      INTRINSIC          DBLE
96*     ..
97*     .. Executable Statements ..
98*
99      NOUT = NUNIT
100      WRITE( NOUT, FMT = * )
101*
102*     Set the variables to innocuous values.
103*
104      DO 20 J = 1, NMAX
105         DO 10 I = 1, NMAX
106            A( I, J ) = 1.D0 / DBLE( I+J )
107            AF( I, J ) = 1.D0 / DBLE( I+J )
108   10    CONTINUE
109         B( J ) = 0.D0
110         W( J ) = 0.D0
111         X( J ) = 0.D0
112   20 CONTINUE
113      OK = .TRUE.
114*
115*     Error exits for QL factorization
116*
117*     DGEQLF
118*
119      SRNAMT = 'DGEQLF'
120      INFOT = 1
121      CALL DGEQLF( -1, 0, A, 1, B, W, 1, INFO )
122      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
123      INFOT = 2
124      CALL DGEQLF( 0, -1, A, 1, B, W, 1, INFO )
125      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
126      INFOT = 4
127      CALL DGEQLF( 2, 1, A, 1, B, W, 1, INFO )
128      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
129      INFOT = 7
130      CALL DGEQLF( 1, 2, A, 1, B, W, 1, INFO )
131      CALL CHKXER( 'DGEQLF', INFOT, NOUT, LERR, OK )
132*
133*     DGEQL2
134*
135      SRNAMT = 'DGEQL2'
136      INFOT = 1
137      CALL DGEQL2( -1, 0, A, 1, B, W, INFO )
138      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
139      INFOT = 2
140      CALL DGEQL2( 0, -1, A, 1, B, W, INFO )
141      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
142      INFOT = 4
143      CALL DGEQL2( 2, 1, A, 1, B, W, INFO )
144      CALL CHKXER( 'DGEQL2', INFOT, NOUT, LERR, OK )
145*
146*     DGEQLS
147*
148      SRNAMT = 'DGEQLS'
149      INFOT = 1
150      CALL DGEQLS( -1, 0, 0, A, 1, X, B, 1, W, 1, INFO )
151      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
152      INFOT = 2
153      CALL DGEQLS( 0, -1, 0, A, 1, X, B, 1, W, 1, INFO )
154      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
155      INFOT = 2
156      CALL DGEQLS( 1, 2, 0, A, 1, X, B, 1, W, 1, INFO )
157      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
158      INFOT = 3
159      CALL DGEQLS( 0, 0, -1, A, 1, X, B, 1, W, 1, INFO )
160      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
161      INFOT = 5
162      CALL DGEQLS( 2, 1, 0, A, 1, X, B, 2, W, 1, INFO )
163      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
164      INFOT = 8
165      CALL DGEQLS( 2, 1, 0, A, 2, X, B, 1, W, 1, INFO )
166      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
167      INFOT = 10
168      CALL DGEQLS( 1, 1, 2, A, 1, X, B, 1, W, 1, INFO )
169      CALL CHKXER( 'DGEQLS', INFOT, NOUT, LERR, OK )
170*
171*     DORGQL
172*
173      SRNAMT = 'DORGQL'
174      INFOT = 1
175      CALL DORGQL( -1, 0, 0, A, 1, X, W, 1, INFO )
176      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
177      INFOT = 2
178      CALL DORGQL( 0, -1, 0, A, 1, X, W, 1, INFO )
179      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
180      INFOT = 2
181      CALL DORGQL( 1, 2, 0, A, 1, X, W, 2, INFO )
182      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
183      INFOT = 3
184      CALL DORGQL( 0, 0, -1, A, 1, X, W, 1, INFO )
185      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
186      INFOT = 3
187      CALL DORGQL( 1, 1, 2, A, 1, X, W, 1, INFO )
188      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
189      INFOT = 5
190      CALL DORGQL( 2, 1, 0, A, 1, X, W, 1, INFO )
191      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
192      INFOT = 8
193      CALL DORGQL( 2, 2, 0, A, 2, X, W, 1, INFO )
194      CALL CHKXER( 'DORGQL', INFOT, NOUT, LERR, OK )
195*
196*     DORG2L
197*
198      SRNAMT = 'DORG2L'
199      INFOT = 1
200      CALL DORG2L( -1, 0, 0, A, 1, X, W, INFO )
201      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
202      INFOT = 2
203      CALL DORG2L( 0, -1, 0, A, 1, X, W, INFO )
204      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
205      INFOT = 2
206      CALL DORG2L( 1, 2, 0, A, 1, X, W, INFO )
207      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
208      INFOT = 3
209      CALL DORG2L( 0, 0, -1, A, 1, X, W, INFO )
210      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
211      INFOT = 3
212      CALL DORG2L( 2, 1, 2, A, 2, X, W, INFO )
213      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
214      INFOT = 5
215      CALL DORG2L( 2, 1, 0, A, 1, X, W, INFO )
216      CALL CHKXER( 'DORG2L', INFOT, NOUT, LERR, OK )
217*
218*     DORMQL
219*
220      SRNAMT = 'DORMQL'
221      INFOT = 1
222      CALL DORMQL( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
223      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
224      INFOT = 2
225      CALL DORMQL( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
226      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
227      INFOT = 3
228      CALL DORMQL( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
229      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
230      INFOT = 4
231      CALL DORMQL( 'L', 'N', 0, -1, 0, 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, 0, -1, A, 1, X, AF, 1, W, 1, INFO )
235      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
236      INFOT = 5
237      CALL DORMQL( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, 1, INFO )
238      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
239      INFOT = 5
240      CALL DORMQL( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, 1, INFO )
241      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
242      INFOT = 7
243      CALL DORMQL( 'L', 'N', 2, 1, 0, A, 1, X, AF, 2, W, 1, INFO )
244      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
245      INFOT = 7
246      CALL DORMQL( 'R', 'N', 1, 2, 0, A, 1, X, AF, 1, W, 1, INFO )
247      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
248      INFOT = 10
249      CALL DORMQL( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, 1, INFO )
250      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
251      INFOT = 12
252      CALL DORMQL( 'L', 'N', 1, 2, 0, A, 1, X, AF, 1, W, 1, INFO )
253      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
254      INFOT = 12
255      CALL DORMQL( 'R', 'N', 2, 1, 0, A, 1, X, AF, 2, W, 1, INFO )
256      CALL CHKXER( 'DORMQL', INFOT, NOUT, LERR, OK )
257*
258*     DORM2L
259*
260      SRNAMT = 'DORM2L'
261      INFOT = 1
262      CALL DORM2L( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
263      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
264      INFOT = 2
265      CALL DORM2L( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
266      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
267      INFOT = 3
268      CALL DORM2L( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, INFO )
269      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
270      INFOT = 4
271      CALL DORM2L( 'L', 'N', 0, -1, 0, A, 1, X, AF, 1, W, INFO )
272      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
273      INFOT = 5
274      CALL DORM2L( 'L', 'N', 0, 0, -1, A, 1, X, AF, 1, W, INFO )
275      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
276      INFOT = 5
277      CALL DORM2L( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, INFO )
278      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
279      INFOT = 5
280      CALL DORM2L( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, INFO )
281      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
282      INFOT = 7
283      CALL DORM2L( 'L', 'N', 2, 1, 0, A, 1, X, AF, 2, W, INFO )
284      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
285      INFOT = 7
286      CALL DORM2L( 'R', 'N', 1, 2, 0, A, 1, X, AF, 1, W, INFO )
287      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
288      INFOT = 10
289      CALL DORM2L( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, INFO )
290      CALL CHKXER( 'DORM2L', INFOT, NOUT, LERR, OK )
291*
292*     Print a summary line.
293*
294      CALL ALAESM( PATH, OK, NOUT )
295*
296      RETURN
297*
298*     End of DERRQL
299*
300      END
301