1*> \brief \b DERRLQ
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 DERRLQ( 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*> DERRLQ tests the error exits for the DOUBLE PRECISION routines
25*> that use the LQ 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 DERRLQ( 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, DGELQ2, DGELQF, DGELQS, DORGL2,
80     $                   DORGLQ, DORML2, DORMLQ
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 LQ factorization
113*
114*     DGELQF
115*
116      SRNAMT = 'DGELQF'
117      INFOT = 1
118      CALL DGELQF( -1, 0, A, 1, B, W, 1, INFO )
119      CALL CHKXER( 'DGELQF', INFOT, NOUT, LERR, OK )
120      INFOT = 2
121      CALL DGELQF( 0, -1, A, 1, B, W, 1, INFO )
122      CALL CHKXER( 'DGELQF', INFOT, NOUT, LERR, OK )
123      INFOT = 4
124      CALL DGELQF( 2, 1, A, 1, B, W, 2, INFO )
125      CALL CHKXER( 'DGELQF', INFOT, NOUT, LERR, OK )
126      INFOT = 7
127      CALL DGELQF( 2, 1, A, 2, B, W, 1, INFO )
128      CALL CHKXER( 'DGELQF', INFOT, NOUT, LERR, OK )
129*
130*     DGELQ2
131*
132      SRNAMT = 'DGELQ2'
133      INFOT = 1
134      CALL DGELQ2( -1, 0, A, 1, B, W, INFO )
135      CALL CHKXER( 'DGELQ2', INFOT, NOUT, LERR, OK )
136      INFOT = 2
137      CALL DGELQ2( 0, -1, A, 1, B, W, INFO )
138      CALL CHKXER( 'DGELQ2', INFOT, NOUT, LERR, OK )
139      INFOT = 4
140      CALL DGELQ2( 2, 1, A, 1, B, W, INFO )
141      CALL CHKXER( 'DGELQ2', INFOT, NOUT, LERR, OK )
142*
143*     DGELQS
144*
145      SRNAMT = 'DGELQS'
146      INFOT = 1
147      CALL DGELQS( -1, 0, 0, A, 1, X, B, 1, W, 1, INFO )
148      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
149      INFOT = 2
150      CALL DGELQS( 0, -1, 0, A, 1, X, B, 1, W, 1, INFO )
151      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
152      INFOT = 2
153      CALL DGELQS( 2, 1, 0, A, 2, X, B, 1, W, 1, INFO )
154      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
155      INFOT = 3
156      CALL DGELQS( 0, 0, -1, A, 1, X, B, 1, W, 1, INFO )
157      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
158      INFOT = 5
159      CALL DGELQS( 2, 2, 0, A, 1, X, B, 2, W, 1, INFO )
160      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
161      INFOT = 8
162      CALL DGELQS( 1, 2, 0, A, 1, X, B, 1, W, 1, INFO )
163      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
164      INFOT = 10
165      CALL DGELQS( 1, 1, 2, A, 1, X, B, 1, W, 1, INFO )
166      CALL CHKXER( 'DGELQS', INFOT, NOUT, LERR, OK )
167*
168*     DORGLQ
169*
170      SRNAMT = 'DORGLQ'
171      INFOT = 1
172      CALL DORGLQ( -1, 0, 0, A, 1, X, W, 1, INFO )
173      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
174      INFOT = 2
175      CALL DORGLQ( 0, -1, 0, A, 1, X, W, 1, INFO )
176      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
177      INFOT = 2
178      CALL DORGLQ( 2, 1, 0, A, 2, X, W, 2, INFO )
179      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
180      INFOT = 3
181      CALL DORGLQ( 0, 0, -1, A, 1, X, W, 1, INFO )
182      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
183      INFOT = 3
184      CALL DORGLQ( 1, 1, 2, A, 1, X, W, 1, INFO )
185      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
186      INFOT = 5
187      CALL DORGLQ( 2, 2, 0, A, 1, X, W, 2, INFO )
188      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
189      INFOT = 8
190      CALL DORGLQ( 2, 2, 0, A, 2, X, W, 1, INFO )
191      CALL CHKXER( 'DORGLQ', INFOT, NOUT, LERR, OK )
192*
193*     DORGL2
194*
195      SRNAMT = 'DORGL2'
196      INFOT = 1
197      CALL DORGL2( -1, 0, 0, A, 1, X, W, INFO )
198      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
199      INFOT = 2
200      CALL DORGL2( 0, -1, 0, A, 1, X, W, INFO )
201      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
202      INFOT = 2
203      CALL DORGL2( 2, 1, 0, A, 2, X, W, INFO )
204      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
205      INFOT = 3
206      CALL DORGL2( 0, 0, -1, A, 1, X, W, INFO )
207      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
208      INFOT = 3
209      CALL DORGL2( 1, 1, 2, A, 1, X, W, INFO )
210      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
211      INFOT = 5
212      CALL DORGL2( 2, 2, 0, A, 1, X, W, INFO )
213      CALL CHKXER( 'DORGL2', INFOT, NOUT, LERR, OK )
214*
215*     DORMLQ
216*
217      SRNAMT = 'DORMLQ'
218      INFOT = 1
219      CALL DORMLQ( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
220      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
221      INFOT = 2
222      CALL DORMLQ( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
223      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
224      INFOT = 3
225      CALL DORMLQ( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, 1, INFO )
226      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
227      INFOT = 4
228      CALL DORMLQ( 'L', 'N', 0, -1, 0, A, 1, X, AF, 1, W, 1, INFO )
229      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
230      INFOT = 5
231      CALL DORMLQ( 'L', 'N', 0, 0, -1, A, 1, X, AF, 1, W, 1, INFO )
232      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
233      INFOT = 5
234      CALL DORMLQ( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, 1, INFO )
235      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
236      INFOT = 5
237      CALL DORMLQ( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, 1, INFO )
238      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
239      INFOT = 7
240      CALL DORMLQ( 'L', 'N', 2, 0, 2, A, 1, X, AF, 2, W, 1, INFO )
241      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
242      INFOT = 7
243      CALL DORMLQ( 'R', 'N', 0, 2, 2, A, 1, X, AF, 1, W, 1, INFO )
244      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
245      INFOT = 10
246      CALL DORMLQ( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, 1, INFO )
247      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
248      INFOT = 12
249      CALL DORMLQ( 'L', 'N', 1, 2, 0, A, 1, X, AF, 1, W, 1, INFO )
250      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
251      INFOT = 12
252      CALL DORMLQ( 'R', 'N', 2, 1, 0, A, 1, X, AF, 2, W, 1, INFO )
253      CALL CHKXER( 'DORMLQ', INFOT, NOUT, LERR, OK )
254*
255*     DORML2
256*
257      SRNAMT = 'DORML2'
258      INFOT = 1
259      CALL DORML2( '/', 'N', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
260      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
261      INFOT = 2
262      CALL DORML2( 'L', '/', 0, 0, 0, A, 1, X, AF, 1, W, INFO )
263      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
264      INFOT = 3
265      CALL DORML2( 'L', 'N', -1, 0, 0, A, 1, X, AF, 1, W, INFO )
266      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
267      INFOT = 4
268      CALL DORML2( 'L', 'N', 0, -1, 0, A, 1, X, AF, 1, W, INFO )
269      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
270      INFOT = 5
271      CALL DORML2( 'L', 'N', 0, 0, -1, A, 1, X, AF, 1, W, INFO )
272      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
273      INFOT = 5
274      CALL DORML2( 'L', 'N', 0, 1, 1, A, 1, X, AF, 1, W, INFO )
275      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
276      INFOT = 5
277      CALL DORML2( 'R', 'N', 1, 0, 1, A, 1, X, AF, 1, W, INFO )
278      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
279      INFOT = 7
280      CALL DORML2( 'L', 'N', 2, 1, 2, A, 1, X, AF, 2, W, INFO )
281      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
282      INFOT = 7
283      CALL DORML2( 'R', 'N', 1, 2, 2, A, 1, X, AF, 1, W, INFO )
284      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
285      INFOT = 10
286      CALL DORML2( 'L', 'N', 2, 1, 0, A, 2, X, AF, 1, W, INFO )
287      CALL CHKXER( 'DORML2', INFOT, NOUT, LERR, OK )
288*
289*     Print a summary line.
290*
291      CALL ALAESM( PATH, OK, NOUT )
292*
293      RETURN
294*
295*     End of DERRLQ
296*
297      END
298