Browse Files

DDRVLS


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

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

Purpose

DDRVLS tests the least squares driver routines DGELS, DGELSS, DGELSX,
DGELSY and DGELSD.

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 orthogonal 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 orthogonal 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.
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.
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.
THRESH
(input) DOUBLE PRECISION
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) DOUBLE PRECISION 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) DOUBLE PRECISION array, dimension (MMAX*NMAX)
B
(workspace) DOUBLE PRECISION 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) DOUBLE PRECISION array, dimension (MMAX*NSMAX)
C
(workspace) DOUBLE PRECISION array, dimension (MMAX*NSMAX)
S
(workspace) DOUBLE PRECISION array, dimension
            (min(MMAX,NMAX))
COPYS
(workspace) DOUBLE PRECISION array, dimension
            (min(MMAX,NMAX))
WORK
(workspace) DOUBLE PRECISION array,
            dimension (MMAX*NMAX + 4*NMAX + MMAX).
IWORK
(workspace) INTEGER array, dimension (15*NMAX)
NOUT
(input) INTEGER
The unit number for output.

Call Graph

Caller Graph