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