diff options
Diffstat (limited to 'TESTING/EIG/cckgqr.f')
| -rw-r--r-- | TESTING/EIG/cckgqr.f | 304 |
1 files changed, 214 insertions, 90 deletions
diff --git a/TESTING/EIG/cckgqr.f b/TESTING/EIG/cckgqr.f index c439a98f..83b8e272 100644 --- a/TESTING/EIG/cckgqr.f +++ b/TESTING/EIG/cckgqr.f @@ -1,10 +1,222 @@ +*> \brief \b CCKGQR +* +* =========== DOCUMENTATION =========== +* +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ +* +* Definition +* ========== +* +* SUBROUTINE CCKGQR( NM, MVAL, NP, PVAL, NN, NVAL, NMATS, ISEED, +* THRESH, NMAX, A, AF, AQ, AR, TAUA, B, BF, BZ, +* BT, BWK, TAUB, WORK, RWORK, NIN, NOUT, INFO ) +* +* .. Scalar Arguments .. +* INTEGER INFO, NIN, NM, NMATS, NMAX, NN, NOUT, NP +* REAL THRESH +* .. +* .. Array Arguments .. +* INTEGER ISEED( 4 ), MVAL( * ), NVAL( * ), PVAL( * ) +* REAL RWORK( * ) +* COMPLEX A( * ), AF( * ), AQ( * ), AR( * ), B( * ), +* $ BF( * ), BT( * ), BWK( * ), BZ( * ), TAUA( * ), +* $ TAUB( * ), WORK( * ) +* .. +* +* Purpose +* ======= +* +*>\details \b Purpose: +*>\verbatim +*> +*> CCKGQR tests +*> CGGQRF: GQR factorization for N-by-M matrix A and N-by-P matrix B, +*> CGGRQF: GRQ factorization for M-by-N matrix A and P-by-N matrix B. +*> +*>\endverbatim +* +* Arguments +* ========= +* +*> \param[in] NM +*> \verbatim +*> NM is INTEGER +*> The number of values of M contained in the vector MVAL. +*> \endverbatim +*> +*> \param[in] MVAL +*> \verbatim +*> MVAL is INTEGER array, dimension (NM) +*> The values of the matrix row(column) dimension M. +*> \endverbatim +*> +*> \param[in] NP +*> \verbatim +*> NP is INTEGER +*> The number of values of P contained in the vector PVAL. +*> \endverbatim +*> +*> \param[in] PVAL +*> \verbatim +*> PVAL is INTEGER array, dimension (NP) +*> The values of the matrix row(column) dimension P. +*> \endverbatim +*> +*> \param[in] NN +*> \verbatim +*> NN is INTEGER +*> The number of values of N contained in the vector NVAL. +*> \endverbatim +*> +*> \param[in] NVAL +*> \verbatim +*> NVAL is INTEGER array, dimension (NN) +*> The values of the matrix column(row) dimension N. +*> \endverbatim +*> +*> \param[in] NMATS +*> \verbatim +*> NMATS is INTEGER +*> The number of matrix types to be tested for each combination +*> of matrix dimensions. If NMATS >= NTYPES (the maximum +*> number of matrix types), then all the different types are +*> generated for testing. If NMATS < NTYPES, another input line +*> is read to get the numbers of the matrix types to be used. +*> \endverbatim +*> +*> \param[in,out] ISEED +*> \verbatim +*> ISEED is INTEGER array, dimension (4) +*> On entry, the seed of the random number generator. The array +*> elements should be between 0 and 4095, otherwise they will be +*> reduced mod 4096, and ISEED(4) must be odd. +*> On exit, the next seed in the random number sequence after +*> all the test matrices have been generated. +*> \endverbatim +*> +*> \param[in] THRESH +*> \verbatim +*> 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. +*> \endverbatim +*> +*> \param[in] NMAX +*> \verbatim +*> NMAX is INTEGER +*> The maximum value permitted for M or N, used in dimensioning +*> the work arrays. +*> \endverbatim +*> +*> \param[out] A +*> \verbatim +*> A is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] AF +*> \verbatim +*> AF is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] AQ +*> \verbatim +*> AQ is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] AR +*> \verbatim +*> AR is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] TAUA +*> \verbatim +*> TAUA is COMPLEX array, dimension (NMAX) +*> \endverbatim +*> +*> \param[out] B +*> \verbatim +*> B is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] BF +*> \verbatim +*> BF is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] BZ +*> \verbatim +*> BZ is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] BT +*> \verbatim +*> BT is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] BWK +*> \verbatim +*> BWK is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] TAUB +*> \verbatim +*> TAUB is COMPLEX array, dimension (NMAX) +*> \endverbatim +*> +*> \param[out] WORK +*> \verbatim +*> WORK is COMPLEX array, dimension (NMAX*NMAX) +*> \endverbatim +*> +*> \param[out] RWORK +*> \verbatim +*> RWORK is REAL array, dimension (NMAX) +*> \endverbatim +*> +*> \param[in] NIN +*> \verbatim +*> NIN is INTEGER +*> The unit number for input. +*> \endverbatim +*> +*> \param[in] NOUT +*> \verbatim +*> NOUT is INTEGER +*> The unit number for output. +*> \endverbatim +*> +*> \param[out] INFO +*> \verbatim +*> INFO is INTEGER +*> = 0 : successful exit +*> > 0 : If CLATMS returns an error code, the absolute value +*> of it is returned. +*> \endverbatim +*> +* +* Authors +* ======= +* +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. +* +*> \date November 2011 +* +*> \ingroup complex_eig +* +* ===================================================================== SUBROUTINE CCKGQR( NM, MVAL, NP, PVAL, NN, NVAL, NMATS, ISEED, $ THRESH, NMAX, A, AF, AQ, AR, TAUA, B, BF, BZ, $ BT, BWK, TAUB, WORK, RWORK, NIN, NOUT, INFO ) * * -- LAPACK test routine (version 3.1) -- -* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. -* November 2006 +* -- LAPACK is a software package provided by Univ. of Tennessee, -- +* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- +* November 2011 * * .. Scalar Arguments .. INTEGER INFO, NIN, NM, NMATS, NMAX, NN, NOUT, NP @@ -18,94 +230,6 @@ $ TAUB( * ), WORK( * ) * .. * -* Purpose -* ======= -* -* CCKGQR tests -* CGGQRF: GQR factorization for N-by-M matrix A and N-by-P matrix B, -* CGGRQF: GRQ factorization for M-by-N matrix A and P-by-N matrix B. -* -* Arguments -* ========= -* -* 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(column) dimension M. -* -* NP (input) INTEGER -* The number of values of P contained in the vector PVAL. -* -* PVAL (input) INTEGER array, dimension (NP) -* The values of the matrix row(column) dimension P. -* -* 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(row) dimension N. -* -* NMATS (input) INTEGER -* The number of matrix types to be tested for each combination -* of matrix dimensions. If NMATS >= NTYPES (the maximum -* number of matrix types), then all the different types are -* generated for testing. If NMATS < NTYPES, another input line -* is read to get the numbers of the matrix types to be used. -* -* ISEED (input/output) INTEGER array, dimension (4) -* On entry, the seed of the random number generator. The array -* elements should be between 0 and 4095, otherwise they will be -* reduced mod 4096, and ISEED(4) must be odd. -* On exit, the next seed in the random number sequence after -* all the test matrices have been generated. -* -* 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. -* -* NMAX (input) INTEGER -* The maximum value permitted for M or N, used in dimensioning -* the work arrays. -* -* A (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* AF (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* AQ (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* AR (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* TAUA (workspace) COMPLEX array, dimension (NMAX) -* -* B (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* BF (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* BZ (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* BT (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* BWK (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* TAUB (workspace) COMPLEX array, dimension (NMAX) -* -* WORK (workspace) COMPLEX array, dimension (NMAX*NMAX) -* -* RWORK (workspace) REAL array, dimension (NMAX) -* -* NIN (input) INTEGER -* The unit number for input. -* -* NOUT (input) INTEGER -* The unit number for output. -* -* INFO (output) INTEGER -* = 0 : successful exit -* > 0 : If CLATMS returns an error code, the absolute value -* of it is returned. -* * ===================================================================== * * .. Parameters .. |