58
59
60
61
62
63
64 CHARACTER*3 PATH
65 INTEGER NUNIT
66
67
68
69
70
71 INTEGER NMAX, LW
72 parameter( nmax = 4, lw = 3*nmax )
73
74
75 CHARACTER EQ
76 CHARACTER*2 C2
77 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
78 REAL ANRM, CCOND, RCOND, BERR
79
80
81 INTEGER IP( NMAX ), IW( NMAX )
82 REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
83 $ C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
84 $ W( LW ), X( NMAX ), ERR_BNDS_N( NMAX, 3 ),
85 $ ERR_BNDS_C( NMAX, 3 ), PARAMS( 1 )
86
87
88 LOGICAL LSAMEN
90
91
96
97
98 LOGICAL LERR, OK
99 CHARACTER*32 SRNAMT
100 INTEGER INFOT, NOUT
101
102
103 COMMON / infoc / infot, nout, ok, lerr
104 COMMON / srnamc / srnamt
105
106
107 INTRINSIC real
108
109
110
111 nout = nunit
112 WRITE( nout, fmt = * )
113 c2 = path( 2: 3 )
114
115
116
117 DO 20 j = 1, nmax
118 DO 10 i = 1, nmax
119 a( i, j ) = 1. / real( i+j )
120 af( i, j ) = 1. / real( i+j )
121 10 CONTINUE
122 b( j ) = 0.
123 r1( j ) = 0.
124 r2( j ) = 0.
125 w( j ) = 0.
126 x( j ) = 0.
127 c( j ) = 0.
128 r( j ) = 0.
129 ip( j ) = j
130 iw( j ) = j
131 20 CONTINUE
132 ok = .true.
133
134 IF(
lsamen( 2, c2,
'GE' ) )
THEN
135
136
137
138
139
140
141 srnamt = 'SGETRF'
142 infot = 1
143 CALL sgetrf( -1, 0, a, 1, ip, info )
144 CALL chkxer(
'SGETRF', infot, nout, lerr, ok )
145 infot = 2
146 CALL sgetrf( 0, -1, a, 1, ip, info )
147 CALL chkxer(
'SGETRF', infot, nout, lerr, ok )
148 infot = 4
149 CALL sgetrf( 2, 1, a, 1, ip, info )
150 CALL chkxer(
'SGETRF', infot, nout, lerr, ok )
151
152
153
154 srnamt = 'SGETF2'
155 infot = 1
156 CALL sgetf2( -1, 0, a, 1, ip, info )
157 CALL chkxer(
'SGETF2', infot, nout, lerr, ok )
158 infot = 2
159 CALL sgetf2( 0, -1, a, 1, ip, info )
160 CALL chkxer(
'SGETF2', infot, nout, lerr, ok )
161 infot = 4
162 CALL sgetf2( 2, 1, a, 1, ip, info )
163 CALL chkxer(
'SGETF2', infot, nout, lerr, ok )
164
165
166
167 srnamt = 'SGETRI'
168 infot = 1
169 CALL sgetri( -1, a, 1, ip, w, lw, info )
170 CALL chkxer(
'SGETRI', infot, nout, lerr, ok )
171 infot = 3
172 CALL sgetri( 2, a, 1, ip, w, lw, info )
173 CALL chkxer(
'SGETRI', infot, nout, lerr, ok )
174
175
176
177 srnamt = 'SGETRS'
178 infot = 1
179 CALL sgetrs(
'/', 0, 0, a, 1, ip, b, 1, info )
180 CALL chkxer(
'SGETRS', infot, nout, lerr, ok )
181 infot = 2
182 CALL sgetrs(
'N', -1, 0, a, 1, ip, b, 1, info )
183 CALL chkxer(
'SGETRS', infot, nout, lerr, ok )
184 infot = 3
185 CALL sgetrs(
'N', 0, -1, a, 1, ip, b, 1, info )
186 CALL chkxer(
'SGETRS', infot, nout, lerr, ok )
187 infot = 5
188 CALL sgetrs(
'N', 2, 1, a, 1, ip, b, 2, info )
189 CALL chkxer(
'SGETRS', infot, nout, lerr, ok )
190 infot = 8
191 CALL sgetrs(
'N', 2, 1, a, 2, ip, b, 1, info )
192 CALL chkxer(
'SGETRS', infot, nout, lerr, ok )
193
194
195
196 srnamt = 'SGERFS'
197 infot = 1
198 CALL sgerfs(
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1, r1, r2, w,
199 $ iw, info )
200 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
201 infot = 2
202 CALL sgerfs(
'N', -1, 0, a, 1, af, 1, ip, b, 1, x, 1, r1, r2,
203 $ w, iw, info )
204 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
205 infot = 3
206 CALL sgerfs(
'N', 0, -1, a, 1, af, 1, ip, b, 1, x, 1, r1, r2,
207 $ w, iw, info )
208 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
209 infot = 5
210 CALL sgerfs(
'N', 2, 1, a, 1, af, 2, ip, b, 2, x, 2, r1, r2, w,
211 $ iw, info )
212 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
213 infot = 7
214 CALL sgerfs(
'N', 2, 1, a, 2, af, 1, ip, b, 2, x, 2, r1, r2, w,
215 $ iw, info )
216 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
217 infot = 10
218 CALL sgerfs(
'N', 2, 1, a, 2, af, 2, ip, b, 1, x, 2, r1, r2, w,
219 $ iw, info )
220 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
221 infot = 12
222 CALL sgerfs(
'N', 2, 1, a, 2, af, 2, ip, b, 2, x, 1, r1, r2, w,
223 $ iw, info )
224 CALL chkxer(
'SGERFS', infot, nout, lerr, ok )
225
226
227
228 n_err_bnds = 3
229 nparams = 0
230 srnamt = 'SGERFSX'
231 infot = 1
232 CALL sgerfsx(
'/', eq, 0, 0, a, 1, af, 1, ip, r, c, b, 1, x,
233 $ 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
234 $ nparams, params, w, iw, info )
235 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
236 infot = 2
237 eq = '/'
238 CALL sgerfsx(
'N', eq, 2, 1, a, 1, af, 2, ip, r, c, b, 2, x,
239 $ 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
240 $ nparams, params, w, iw, info )
241 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
242 infot = 3
243 eq = 'R'
244 CALL sgerfsx(
'N', eq, -1, 0, a, 1, af, 1, ip, r, c, b, 1, x,
245 $ 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
246 $ nparams, params, w, iw, info )
247 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
248 infot = 4
249 CALL sgerfsx(
'N', eq, 0, -1, a, 1, af, 1, ip, r, c, b, 1, x,
250 $ 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
251 $ nparams, params, w, iw, info )
252 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
253 infot = 6
254 CALL sgerfsx(
'N', eq, 2, 1, a, 1, af, 2, ip, r, c, b, 2, x,
255 $ 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
256 $ nparams, params, w, iw, info )
257 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
258 infot = 8
259 CALL sgerfsx(
'N', eq, 2, 1, a, 2, af, 1, ip, r, c, b, 2, x,
260 $ 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
261 $ nparams, params, w, iw, info )
262 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
263 infot = 13
264 eq = 'C'
265 CALL sgerfsx(
'N', eq, 2, 1, a, 2, af, 2, ip, r, c, b, 1, x,
266 $ 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
267 $ nparams, params, w, iw, info )
268 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
269 infot = 15
270 CALL sgerfsx(
'N', eq, 2, 1, a, 2, af, 2, ip, r, c, b, 2, x,
271 $ 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
272 $ nparams, params, w, iw, info )
273 CALL chkxer(
'SGERFSX', infot, nout, lerr, ok )
274
275
276
277 srnamt = 'SGECON'
278 infot = 1
279 CALL sgecon(
'/', 0, a, 1, anrm, rcond, w, iw, info )
280 CALL chkxer(
'SGECON', infot, nout, lerr, ok )
281 infot = 2
282 CALL sgecon(
'1', -1, a, 1, anrm, rcond, w, iw, info )
283 CALL chkxer(
'SGECON', infot, nout, lerr, ok )
284 infot = 4
285 CALL sgecon(
'1', 2, a, 1, anrm, rcond, w, iw, info )
286 CALL chkxer(
'SGECON', infot, nout, lerr, ok )
287
288
289
290 srnamt = 'SGEEQU'
291 infot = 1
292 CALL sgeequ( -1, 0, a, 1, r1, r2, rcond, ccond, anrm, info )
293 CALL chkxer(
'SGEEQU', infot, nout, lerr, ok )
294 infot = 2
295 CALL sgeequ( 0, -1, a, 1, r1, r2, rcond, ccond, anrm, info )
296 CALL chkxer(
'SGEEQU', infot, nout, lerr, ok )
297 infot = 4
298 CALL sgeequ( 2, 2, a, 1, r1, r2, rcond, ccond, anrm, info )
299 CALL chkxer(
'SGEEQU', infot, nout, lerr, ok )
300
301
302
303 srnamt = 'SGEEQUB'
304 infot = 1
305 CALL sgeequb( -1, 0, a, 1, r1, r2, rcond, ccond, anrm, info )
306 CALL chkxer(
'SGEEQUB', infot, nout, lerr, ok )
307 infot = 2
308 CALL sgeequb( 0, -1, a, 1, r1, r2, rcond, ccond, anrm, info )
309 CALL chkxer(
'SGEEQUB', infot, nout, lerr, ok )
310 infot = 4
311 CALL sgeequb( 2, 2, a, 1, r1, r2, rcond, ccond, anrm, info )
312 CALL chkxer(
'SGEEQUB', infot, nout, lerr, ok )
313
314 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
315
316
317
318
319
320
321 srnamt = 'SGBTRF'
322 infot = 1
323 CALL sgbtrf( -1, 0, 0, 0, a, 1, ip, info )
324 CALL chkxer(
'SGBTRF', infot, nout, lerr, ok )
325 infot = 2
326 CALL sgbtrf( 0, -1, 0, 0, a, 1, ip, info )
327 CALL chkxer(
'SGBTRF', infot, nout, lerr, ok )
328 infot = 3
329 CALL sgbtrf( 1, 1, -1, 0, a, 1, ip, info )
330 CALL chkxer(
'SGBTRF', infot, nout, lerr, ok )
331 infot = 4
332 CALL sgbtrf( 1, 1, 0, -1, a, 1, ip, info )
333 CALL chkxer(
'SGBTRF', infot, nout, lerr, ok )
334 infot = 6
335 CALL sgbtrf( 2, 2, 1, 1, a, 3, ip, info )
336 CALL chkxer(
'SGBTRF', infot, nout, lerr, ok )
337
338
339
340 srnamt = 'SGBTF2'
341 infot = 1
342 CALL sgbtf2( -1, 0, 0, 0, a, 1, ip, info )
343 CALL chkxer(
'SGBTF2', infot, nout, lerr, ok )
344 infot = 2
345 CALL sgbtf2( 0, -1, 0, 0, a, 1, ip, info )
346 CALL chkxer(
'SGBTF2', infot, nout, lerr, ok )
347 infot = 3
348 CALL sgbtf2( 1, 1, -1, 0, a, 1, ip, info )
349 CALL chkxer(
'SGBTF2', infot, nout, lerr, ok )
350 infot = 4
351 CALL sgbtf2( 1, 1, 0, -1, a, 1, ip, info )
352 CALL chkxer(
'SGBTF2', infot, nout, lerr, ok )
353 infot = 6
354 CALL sgbtf2( 2, 2, 1, 1, a, 3, ip, info )
355 CALL chkxer(
'SGBTF2', infot, nout, lerr, ok )
356
357
358
359 srnamt = 'SGBTRS'
360 infot = 1
361 CALL sgbtrs(
'/', 0, 0, 0, 1, a, 1, ip, b, 1, info )
362 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
363 infot = 2
364 CALL sgbtrs(
'N', -1, 0, 0, 1, a, 1, ip, b, 1, info )
365 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
366 infot = 3
367 CALL sgbtrs(
'N', 1, -1, 0, 1, a, 1, ip, b, 1, info )
368 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
369 infot = 4
370 CALL sgbtrs(
'N', 1, 0, -1, 1, a, 1, ip, b, 1, info )
371 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
372 infot = 5
373 CALL sgbtrs(
'N', 1, 0, 0, -1, a, 1, ip, b, 1, info )
374 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
375 infot = 7
376 CALL sgbtrs(
'N', 2, 1, 1, 1, a, 3, ip, b, 2, info )
377 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
378 infot = 10
379 CALL sgbtrs(
'N', 2, 0, 0, 1, a, 1, ip, b, 1, info )
380 CALL chkxer(
'SGBTRS', infot, nout, lerr, ok )
381
382
383
384 srnamt = 'SGBRFS'
385 infot = 1
386 CALL sgbrfs(
'/', 0, 0, 0, 0, a, 1, af, 1, ip, b, 1, x, 1, r1,
387 $ r2, w, iw, info )
388 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
389 infot = 2
390 CALL sgbrfs(
'N', -1, 0, 0, 0, a, 1, af, 1, ip, b, 1, x, 1, r1,
391 $ r2, w, iw, info )
392 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
393 infot = 3
394 CALL sgbrfs(
'N', 1, -1, 0, 0, a, 1, af, 1, ip, b, 1, x, 1, r1,
395 $ r2, w, iw, info )
396 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
397 infot = 4
398 CALL sgbrfs(
'N', 1, 0, -1, 0, a, 1, af, 1, ip, b, 1, x, 1, r1,
399 $ r2, w, iw, info )
400 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
401 infot = 5
402 CALL sgbrfs(
'N', 1, 0, 0, -1, a, 1, af, 1, ip, b, 1, x, 1, r1,
403 $ r2, w, iw, info )
404 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
405 infot = 7
406 CALL sgbrfs(
'N', 2, 1, 1, 1, a, 2, af, 4, ip, b, 2, x, 2, r1,
407 $ r2, w, iw, info )
408 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
409 infot = 9
410 CALL sgbrfs(
'N', 2, 1, 1, 1, a, 3, af, 3, ip, b, 2, x, 2, r1,
411 $ r2, w, iw, info )
412 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
413 infot = 12
414 CALL sgbrfs(
'N', 2, 0, 0, 1, a, 1, af, 1, ip, b, 1, x, 2, r1,
415 $ r2, w, iw, info )
416 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
417 infot = 14
418 CALL sgbrfs(
'N', 2, 0, 0, 1, a, 1, af, 1, ip, b, 2, x, 1, r1,
419 $ r2, w, iw, info )
420 CALL chkxer(
'SGBRFS', infot, nout, lerr, ok )
421
422
423
424 n_err_bnds = 3
425 nparams = 0
426 srnamt = 'SGBRFSX'
427 infot = 1
428 CALL sgbrfsx(
'/', eq, 0, 0, 0, 0, a, 1, af, 1, ip, r, c, b, 1,
429 $ x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
430 $ nparams, params, w, iw, info )
431 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
432 infot = 2
433 eq = '/'
434 CALL sgbrfsx(
'N', eq, 2, 1, 1, 1, a, 1, af, 2, ip, r, c, b, 2,
435 $ x, 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
436 $ nparams, params, w, iw, info )
437 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
438 infot = 3
439 eq = 'R'
440 CALL sgbrfsx(
'N', eq, -1, 1, 1, 0, a, 1, af, 1, ip, r, c, b,
441 $ 1, x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
442 $ nparams, params, w, iw, info )
443 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
444 infot = 4
445 eq = 'R'
446 CALL sgbrfsx(
'N', eq, 2, -1, 1, 1, a, 3, af, 4, ip, r, c, b,
447 $ 1, x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
448 $ nparams, params, w, iw, info )
449 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
450 infot = 5
451 eq = 'R'
452 CALL sgbrfsx(
'N', eq, 2, 1, -1, 1, a, 3, af, 4, ip, r, c, b,
453 $ 1, x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
454 $ nparams, params, w, iw, info )
455 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
456 infot = 6
457 CALL sgbrfsx(
'N', eq, 0, 0, 0, -1, a, 1, af, 1, ip, r, c, b,
458 $ 1, x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
459 $ nparams, params, w, iw, info )
460 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
461 infot = 8
462 CALL sgbrfsx(
'N', eq, 2, 1, 1, 1, a, 1, af, 2, ip, r, c, b,
463 $ 2, x, 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
464 $ nparams, params, w, iw, info )
465 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
466 infot = 10
467 CALL sgbrfsx(
'N', eq, 2, 1, 1, 1, a, 3, af, 3, ip, r, c, b, 2,
468 $ x, 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
469 $ nparams, params, w, iw, info )
470 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
471 infot = 13
472 eq = 'C'
473 CALL sgbrfsx(
'N', eq, 2, 1, 1, 1, a, 3, af, 5, ip, r, c, b,
474 $ 1, x, 2, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
475 $ nparams, params, w, iw, info )
476 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
477 infot = 15
478 CALL sgbrfsx(
'N', eq, 2, 1, 1, 1, a, 3, af, 5, ip, r, c, b, 2,
479 $ x, 1, rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c,
480 $ nparams, params, w, iw, info )
481 CALL chkxer(
'SGBRFSX', infot, nout, lerr, ok )
482
483
484
485 srnamt = 'SGBCON'
486 infot = 1
487 CALL sgbcon(
'/', 0, 0, 0, a, 1, ip, anrm, rcond, w, iw, info )
488 CALL chkxer(
'SGBCON', infot, nout, lerr, ok )
489 infot = 2
490 CALL sgbcon(
'1', -1, 0, 0, a, 1, ip, anrm, rcond, w, iw,
491 $ info )
492 CALL chkxer(
'SGBCON', infot, nout, lerr, ok )
493 infot = 3
494 CALL sgbcon(
'1', 1, -1, 0, a, 1, ip, anrm, rcond, w, iw,
495 $ info )
496 CALL chkxer(
'SGBCON', infot, nout, lerr, ok )
497 infot = 4
498 CALL sgbcon(
'1', 1, 0, -1, a, 1, ip, anrm, rcond, w, iw,
499 $ info )
500 CALL chkxer(
'SGBCON', infot, nout, lerr, ok )
501 infot = 6
502 CALL sgbcon(
'1', 2, 1, 1, a, 3, ip, anrm, rcond, w, iw, info )
503 CALL chkxer(
'SGBCON', infot, nout, lerr, ok )
504
505
506
507 srnamt = 'SGBEQU'
508 infot = 1
509 CALL sgbequ( -1, 0, 0, 0, a, 1, r1, r2, rcond, ccond, anrm,
510 $ info )
511 CALL chkxer(
'SGBEQU', infot, nout, lerr, ok )
512 infot = 2
513 CALL sgbequ( 0, -1, 0, 0, a, 1, r1, r2, rcond, ccond, anrm,
514 $ info )
515 CALL chkxer(
'SGBEQU', infot, nout, lerr, ok )
516 infot = 3
517 CALL sgbequ( 1, 1, -1, 0, a, 1, r1, r2, rcond, ccond, anrm,
518 $ info )
519 CALL chkxer(
'SGBEQU', infot, nout, lerr, ok )
520 infot = 4
521 CALL sgbequ( 1, 1, 0, -1, a, 1, r1, r2, rcond, ccond, anrm,
522 $ info )
523 CALL chkxer(
'SGBEQU', infot, nout, lerr, ok )
524 infot = 6
525 CALL sgbequ( 2, 2, 1, 1, a, 2, r1, r2, rcond, ccond, anrm,
526 $ info )
527 CALL chkxer(
'SGBEQU', infot, nout, lerr, ok )
528
529
530
531 srnamt = 'SGBEQUB'
532 infot = 1
533 CALL sgbequb( -1, 0, 0, 0, a, 1, r1, r2, rcond, ccond, anrm,
534 $ info )
535 CALL chkxer(
'SGBEQUB', infot, nout, lerr, ok )
536 infot = 2
537 CALL sgbequb( 0, -1, 0, 0, a, 1, r1, r2, rcond, ccond, anrm,
538 $ info )
539 CALL chkxer(
'SGBEQUB', infot, nout, lerr, ok )
540 infot = 3
541 CALL sgbequb( 1, 1, -1, 0, a, 1, r1, r2, rcond, ccond, anrm,
542 $ info )
543 CALL chkxer(
'SGBEQUB', infot, nout, lerr, ok )
544 infot = 4
545 CALL sgbequb( 1, 1, 0, -1, a, 1, r1, r2, rcond, ccond, anrm,
546 $ info )
547 CALL chkxer(
'SGBEQUB', infot, nout, lerr, ok )
548 infot = 6
549 CALL sgbequb( 2, 2, 1, 1, a, 2, r1, r2, rcond, ccond, anrm,
550 $ info )
551 CALL chkxer(
'SGBEQUB', infot, nout, lerr, ok )
552 END IF
553
554
555
556 CALL alaesm( path, ok, nout )
557
558 RETURN
559
560
561
subroutine alaesm(path, ok, nout)
ALAESM
subroutine chkxer(srnamt, infot, nout, lerr, ok)
subroutine sgbcon(norm, n, kl, ku, ab, ldab, ipiv, anorm, rcond, work, iwork, info)
SGBCON
subroutine sgbequ(m, n, kl, ku, ab, ldab, r, c, rowcnd, colcnd, amax, info)
SGBEQU
subroutine sgbequb(m, n, kl, ku, ab, ldab, r, c, rowcnd, colcnd, amax, info)
SGBEQUB
subroutine sgbrfs(trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, b, ldb, x, ldx, ferr, berr, work, iwork, info)
SGBRFS
subroutine sgbrfsx(trans, equed, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, r, c, b, ldb, x, ldx, rcond, berr, n_err_bnds, err_bnds_norm, err_bnds_comp, nparams, params, work, iwork, info)
SGBRFSX
subroutine sgbtf2(m, n, kl, ku, ab, ldab, ipiv, info)
SGBTF2 computes the LU factorization of a general band matrix using the unblocked version of the algo...
subroutine sgbtrf(m, n, kl, ku, ab, ldab, ipiv, info)
SGBTRF
subroutine sgbtrs(trans, n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
SGBTRS
subroutine sgecon(norm, n, a, lda, anorm, rcond, work, iwork, info)
SGECON
subroutine sgeequ(m, n, a, lda, r, c, rowcnd, colcnd, amax, info)
SGEEQU
subroutine sgeequb(m, n, a, lda, r, c, rowcnd, colcnd, amax, info)
SGEEQUB
subroutine sgerfs(trans, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, ferr, berr, work, iwork, info)
SGERFS
subroutine sgerfsx(trans, equed, n, nrhs, a, lda, af, ldaf, ipiv, r, c, b, ldb, x, ldx, rcond, berr, n_err_bnds, err_bnds_norm, err_bnds_comp, nparams, params, work, iwork, info)
SGERFSX
subroutine sgetf2(m, n, a, lda, ipiv, info)
SGETF2 computes the LU factorization of a general m-by-n matrix using partial pivoting with row inter...
subroutine sgetrf(m, n, a, lda, ipiv, info)
SGETRF
subroutine sgetri(n, a, lda, ipiv, work, lwork, info)
SGETRI
subroutine sgetrs(trans, n, nrhs, a, lda, ipiv, b, ldb, info)
SGETRS
logical function lsamen(n, ca, cb)
LSAMEN