LAPACK  3.10.0 LAPACK: Linear Algebra PACKage

## ◆ cchkrq()

 subroutine cchkrq ( logical, dimension( * ) DOTYPE, integer NM, integer, dimension( * ) MVAL, integer NN, integer, dimension( * ) NVAL, integer NNB, integer, dimension( * ) NBVAL, integer, dimension( * ) NXVAL, integer NRHS, real THRESH, logical TSTERR, integer NMAX, complex, dimension( * ) A, complex, dimension( * ) AF, complex, dimension( * ) AQ, complex, dimension( * ) AR, complex, dimension( * ) AC, complex, dimension( * ) B, complex, dimension( * ) X, complex, dimension( * ) XACT, complex, dimension( * ) TAU, complex, dimension( * ) WORK, real, dimension( * ) RWORK, integer, dimension( * ) IWORK, integer NOUT )

CCHKRQ

Purpose:
` CCHKRQ tests CGERQF, CUNGRQ and CUNMRQ.`
Parameters
 [in] DOTYPE ``` DOTYPE is LOGICAL array, dimension (NTYPES) The matrix types to be used for testing. Matrices of type j (for 1 <= j <= NTYPES) are used for testing if DOTYPE(j) = .TRUE.; if DOTYPE(j) = .FALSE., then type j is not used.``` [in] NM ``` NM is INTEGER The number of values of M contained in the vector MVAL.``` [in] MVAL ``` MVAL is INTEGER array, dimension (NM) The values of the matrix row dimension M.``` [in] NN ``` NN is INTEGER The number of values of N contained in the vector NVAL.``` [in] NVAL ``` NVAL is INTEGER array, dimension (NN) The values of the matrix column dimension N.``` [in] NNB ``` NNB is INTEGER The number of values of NB and NX contained in the vectors NBVAL and NXVAL. The blocking parameters are used in pairs (NB,NX).``` [in] NBVAL ``` NBVAL is INTEGER array, dimension (NNB) The values of the blocksize NB.``` [in] NXVAL ``` NXVAL is INTEGER array, dimension (NNB) The values of the crossover point NX.``` [in] NRHS ``` NRHS is INTEGER The number of right hand side vectors to be generated for each linear system.``` [in] THRESH ``` THRESH is REAL The threshold value for the test ratios. A result is included in the output file if RESULT >= THRESH. To have every test ratio printed, use THRESH = 0.``` [in] TSTERR ``` TSTERR is LOGICAL Flag that indicates whether error exits are to be tested.``` [in] NMAX ``` NMAX is INTEGER The maximum value permitted for M or N, used in dimensioning the work arrays.``` [out] A ` A is COMPLEX array, dimension (NMAX*NMAX)` [out] AF ` AF is COMPLEX array, dimension (NMAX*NMAX)` [out] AQ ` AQ is COMPLEX array, dimension (NMAX*NMAX)` [out] AR ` AR is COMPLEX array, dimension (NMAX*NMAX)` [out] AC ` AC is COMPLEX array, dimension (NMAX*NMAX)` [out] B ` B is COMPLEX array, dimension (NMAX*NRHS)` [out] X ` X is COMPLEX array, dimension (NMAX*NRHS)` [out] XACT ` XACT is COMPLEX array, dimension (NMAX*NRHS)` [out] TAU ` TAU is COMPLEX array, dimension (NMAX)` [out] WORK ` WORK is COMPLEX array, dimension (NMAX*NMAX)` [out] RWORK ` RWORK is REAL array, dimension (NMAX)` [out] IWORK ` IWORK is INTEGER array, dimension (NMAX)` [in] NOUT ``` NOUT is INTEGER The unit number for output.```

Definition at line 198 of file cchkrq.f.

201 *
202 * -- LAPACK test routine --
203 * -- LAPACK is a software package provided by Univ. of Tennessee, --
204 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
205 *
206 * .. Scalar Arguments ..
207  LOGICAL TSTERR
208  INTEGER NM, NMAX, NN, NNB, NOUT, NRHS
209  REAL THRESH
210 * ..
211 * .. Array Arguments ..
212  LOGICAL DOTYPE( * )
213  INTEGER IWORK( * ), MVAL( * ), NBVAL( * ), NVAL( * ),
214  \$ NXVAL( * )
215  REAL RWORK( * )
216  COMPLEX A( * ), AC( * ), AF( * ), AQ( * ), AR( * ),
217  \$ B( * ), TAU( * ), WORK( * ), X( * ), XACT( * )
218 * ..
219 *
220 * =====================================================================
221 *
222 * .. Parameters ..
223  INTEGER NTESTS
224  parameter( ntests = 7 )
225  INTEGER NTYPES
226  parameter( ntypes = 8 )
227  REAL ZERO
228  parameter( zero = 0.0e0 )
229 * ..
230 * .. Local Scalars ..
231  CHARACTER DIST, TYPE
232  CHARACTER*3 PATH
233  INTEGER I, IK, IM, IMAT, IN, INB, INFO, K, KL, KU, LDA,
234  \$ LWORK, M, MINMN, MODE, N, NB, NERRS, NFAIL, NK,
235  \$ NRUN, NT, NX
236  REAL ANORM, CNDNUM
237 * ..
238 * .. Local Arrays ..
239  INTEGER ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
240  REAL RESULT( NTESTS )
241 * ..
242 * .. External Subroutines ..
243  EXTERNAL alaerh, alahd, alasum, cerrrq, cgerqs, cget02,
245  \$ crqt03, xlaenv
246 * ..
247 * .. Intrinsic Functions ..
248  INTRINSIC max, min
249 * ..
250 * .. Scalars in Common ..
251  LOGICAL LERR, OK
252  CHARACTER*32 SRNAMT
253  INTEGER INFOT, NUNIT
254 * ..
255 * .. Common blocks ..
256  COMMON / infoc / infot, nunit, ok, lerr
257  COMMON / srnamc / srnamt
258 * ..
259 * .. Data statements ..
260  DATA iseedy / 1988, 1989, 1990, 1991 /
261 * ..
262 * .. Executable Statements ..
263 *
264 * Initialize constants and the random number seed.
265 *
266  path( 1: 1 ) = 'Complex precision'
267  path( 2: 3 ) = 'RQ'
268  nrun = 0
269  nfail = 0
270  nerrs = 0
271  DO 10 i = 1, 4
272  iseed( i ) = iseedy( i )
273  10 CONTINUE
274 *
275 * Test the error exits
276 *
277  IF( tsterr )
278  \$ CALL cerrrq( path, nout )
279  infot = 0
280  CALL xlaenv( 2, 2 )
281 *
282  lda = nmax
283  lwork = nmax*max( nmax, nrhs )
284 *
285 * Do for each value of M in MVAL.
286 *
287  DO 70 im = 1, nm
288  m = mval( im )
289 *
290 * Do for each value of N in NVAL.
291 *
292  DO 60 in = 1, nn
293  n = nval( in )
294  minmn = min( m, n )
295  DO 50 imat = 1, ntypes
296 *
297 * Do the tests only if DOTYPE( IMAT ) is true.
298 *
299  IF( .NOT.dotype( imat ) )
300  \$ GO TO 50
301 *
302 * Set up parameters with CLATB4 and generate a test matrix
303 * with CLATMS.
304 *
305  CALL clatb4( path, imat, m, n, TYPE, KL, KU, ANORM, MODE,
306  \$ CNDNUM, DIST )
307 *
308  srnamt = 'CLATMS'
309  CALL clatms( m, n, dist, iseed, TYPE, RWORK, MODE,
310  \$ CNDNUM, ANORM, KL, KU, 'No packing', A, LDA,
311  \$ WORK, INFO )
312 *
313 * Check error code from CLATMS.
314 *
315  IF( info.NE.0 ) THEN
316  CALL alaerh( path, 'CLATMS', info, 0, ' ', m, n, -1,
317  \$ -1, -1, imat, nfail, nerrs, nout )
318  GO TO 50
319  END IF
320 *
321 * Set some values for K: the first value must be MINMN,
322 * corresponding to the call of CRQT01; other values are
323 * used in the calls of CRQT02, and must not exceed MINMN.
324 *
325  kval( 1 ) = minmn
326  kval( 2 ) = 0
327  kval( 3 ) = 1
328  kval( 4 ) = minmn / 2
329  IF( minmn.EQ.0 ) THEN
330  nk = 1
331  ELSE IF( minmn.EQ.1 ) THEN
332  nk = 2
333  ELSE IF( minmn.LE.3 ) THEN
334  nk = 3
335  ELSE
336  nk = 4
337  END IF
338 *
339 * Do for each value of K in KVAL
340 *
341  DO 40 ik = 1, nk
342  k = kval( ik )
343 *
344 * Do for each pair of values (NB,NX) in NBVAL and NXVAL.
345 *
346  DO 30 inb = 1, nnb
347  nb = nbval( inb )
348  CALL xlaenv( 1, nb )
349  nx = nxval( inb )
350  CALL xlaenv( 3, nx )
351  DO i = 1, ntests
352  result( i ) = zero
353  END DO
354  nt = 2
355  IF( ik.EQ.1 ) THEN
356 *
357 * Test CGERQF
358 *
359  CALL crqt01( m, n, a, af, aq, ar, lda, tau,
360  \$ work, lwork, rwork, result( 1 ) )
361  ELSE IF( m.LE.n ) THEN
362 *
363 * Test CUNGRQ, using factorization
364 * returned by CRQT01
365 *
366  CALL crqt02( m, n, k, a, af, aq, ar, lda, tau,
367  \$ work, lwork, rwork, result( 1 ) )
368  END IF
369  IF( m.GE.k ) THEN
370 *
371 * Test CUNMRQ, using factorization returned
372 * by CRQT01
373 *
374  CALL crqt03( m, n, k, af, ac, ar, aq, lda, tau,
375  \$ work, lwork, rwork, result( 3 ) )
376  nt = nt + 4
377 *
378 * If M>=N and K=N, call CGERQS to solve a system
379 * with NRHS right hand sides and compute the
380 * residual.
381 *
382  IF( k.EQ.m .AND. inb.EQ.1 ) THEN
383 *
384 * Generate a solution and set the right
385 * hand side.
386 *
387  srnamt = 'CLARHS'
388  CALL clarhs( path, 'New', 'Full',
389  \$ 'No transpose', m, n, 0, 0,
390  \$ nrhs, a, lda, xact, lda, b, lda,
391  \$ iseed, info )
392 *
393  CALL clacpy( 'Full', m, nrhs, b, lda,
394  \$ x( n-m+1 ), lda )
395  srnamt = 'CGERQS'
396  CALL cgerqs( m, n, nrhs, af, lda, tau, x,
397  \$ lda, work, lwork, info )
398 *
399 * Check error code from CGERQS.
400 *
401  IF( info.NE.0 )
402  \$ CALL alaerh( path, 'CGERQS', info, 0, ' ',
403  \$ m, n, nrhs, -1, nb, imat,
404  \$ nfail, nerrs, nout )
405 *
406  CALL cget02( 'No transpose', m, n, nrhs, a,
407  \$ lda, x, lda, b, lda, rwork,
408  \$ result( 7 ) )
409  nt = nt + 1
410  END IF
411  END IF
412 *
413 * Print information about the tests that did not
414 * pass the threshold.
415 *
416  DO 20 i = 1, nt
417  IF( result( i ).GE.thresh ) THEN
418  IF( nfail.EQ.0 .AND. nerrs.EQ.0 )
419  \$ CALL alahd( nout, path )
420  WRITE( nout, fmt = 9999 )m, n, k, nb, nx,
421  \$ imat, i, result( i )
422  nfail = nfail + 1
423  END IF
424  20 CONTINUE
425  nrun = nrun + nt
426  30 CONTINUE
427  40 CONTINUE
428  50 CONTINUE
429  60 CONTINUE
430  70 CONTINUE
431 *
432 * Print a summary of the results.
433 *
434  CALL alasum( path, nout, nfail, nrun, nerrs )
435 *
436  9999 FORMAT( ' M=', i5, ', N=', i5, ', K=', i5, ', NB=', i4, ', NX=',
437  \$ i5, ', type ', i2, ', test(', i2, ')=', g12.5 )
438  RETURN
439 *
440 * End of CCHKRQ
441 *
subroutine alasum(TYPE, NOUT, NFAIL, NRUN, NERRS)
ALASUM
Definition: alasum.f:73
subroutine xlaenv(ISPEC, NVALUE)
XLAENV
Definition: xlaenv.f:81
subroutine alahd(IOUNIT, PATH)
ALAHD
Definition: alahd.f:107
subroutine alaerh(PATH, SUBNAM, INFO, INFOE, OPTS, M, N, KL, KU, N5, IMAT, NFAIL, NERRS, NOUT)
ALAERH
Definition: alaerh.f:147
subroutine clarhs(PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B, LDB, ISEED, INFO)
CLARHS
Definition: clarhs.f:208
subroutine cget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
CGET02
Definition: cget02.f:134
subroutine clatb4(PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE, CNDNUM, DIST)
CLATB4
Definition: clatb4.f:121
subroutine crqt03(M, N, K, AF, C, CC, Q, LDA, TAU, WORK, LWORK, RWORK, RESULT)
CRQT03
Definition: crqt03.f:136
subroutine crqt01(M, N, A, AF, Q, R, LDA, TAU, WORK, LWORK, RWORK, RESULT)
CRQT01
Definition: crqt01.f:126
subroutine cgerqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
CGERQS
Definition: cgerqs.f:122
subroutine crqt02(M, N, K, A, AF, Q, R, LDA, TAU, WORK, LWORK, RWORK, RESULT)
CRQT02
Definition: crqt02.f:136
subroutine cerrrq(PATH, NUNIT)
CERRRQ
Definition: cerrrq.f:55
subroutine clatms(M, N, DIST, ISEED, SYM, D, MODE, COND, DMAX, KL, KU, PACK, A, LDA, WORK, INFO)
CLATMS
Definition: clatms.f:332
subroutine clacpy(UPLO, M, N, A, LDA, B, LDB)
CLACPY copies all or part of one two-dimensional array to another.
Definition: clacpy.f:103
Here is the call graph for this function:
Here is the caller graph for this function: