Browse Files

CDRVLS


      SUBROUTINE CDRVLS( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL, NNB,
     $                   NBVAL, NXVAL, THRESH, TSTERR, A, COPYA, B,
     $                   COPYB, C, S, COPYS, WORK, RWORK, IWORK,
     $                   NOUT )

   Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
   January 2007

Purpose

CDRVLS tests the least squares driver routines CGELS, CGELSX, CGELSS,
CGELSY and CGELSD.

Arguments

DOTYPE
(input) 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.
The matrix of type j is generated as follows:
j=1: A = U*D*V where U and V are random unitary matrices
     and D has random entries (> 0.1) taken from a uniform
     distribution (0,1). A is full rank.
j=2: The same of 1, but A is scaled up.
j=3: The same of 1, but A is scaled down.
j=4: A = U*D*V where U and V are random unitary matrices
     and D has 3*min(M,N)/4 random entries (> 0.1) taken
     from a uniform distribution (0,1) and the remaining
     entries set to 0. A is rank-deficient.
j=5: The same of 4, but A is scaled up.
j=6: The same of 5, but A is scaled down.
NM
(input) INTEGER
The number of values of M contained in the vector MVAL.
MVAL
(input) INTEGER array, dimension (NM)
The values of the matrix row dimension M.
NN
(input) INTEGER
The number of values of N contained in the vector NVAL.
NVAL
(input) INTEGER array, dimension (NN)
The values of the matrix column dimension N.
NNB
(input) 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).
NBVAL
(input) INTEGER array, dimension (NNB)
The values of the blocksize NB.
NXVAL
(input) INTEGER array, dimension (NNB)
The values of the crossover point NX.
NNS
(input) INTEGER
The number of values of NRHS contained in the vector NSVAL.
NSVAL
(input) INTEGER array, dimension (NNS)
The values of the number of right hand sides NRHS.
THRESH
(input) 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.
TSTERR
(input) LOGICAL
Flag that indicates whether error exits are to be tested.
A
(workspace) COMPLEX array, dimension (MMAX*NMAX)
where MMAX is the maximum value of M in MVAL and NMAX is the
maximum value of N in NVAL.
COPYA
(workspace) COMPLEX array, dimension (MMAX*NMAX)
B
(workspace) COMPLEX array, dimension (MMAX*NSMAX)
where MMAX is the maximum value of M in MVAL and NSMAX is the
maximum value of NRHS in NSVAL.
COPYB
(workspace) COMPLEX array, dimension (MMAX*NSMAX)
C
(workspace) COMPLEX array, dimension (MMAX*NSMAX)
S
(workspace) REAL array, dimension
            (min(MMAX,NMAX))
COPYS
(workspace) REAL array, dimension
            (min(MMAX,NMAX))
WORK
(workspace) COMPLEX array, dimension
            (MMAX*NMAX + 4*NMAX + MMAX).
RWORK
(workspace) REAL array, dimension (5*NMAX-1)
IWORK
(workspace) INTEGER array, dimension (15*NMAX)
NOUT
(input) INTEGER
The unit number for output.

Call Graph

Caller Graph