xGESVDQ tests of the error exits

Adding the first set of tests for xGESVDQ

Author: julielangou

TESTING/EIG/cerred.f

@@ -36,6 +36,8 @@
 *>       CGEJSV   compute SVD of an M-by-N matrix A where M >= N
 *>       CGESVDX  compute SVD of an M-by-N matrix A(by bisection
 *>                and inverse iteration)
+*>       CGESVDQ  compute SVD of an M-by-N matrix A(with a 
+*>                QR-Preconditioned )
 *> \endverbatim
 *
 *  Arguments:
@@ -101,7 +103,7 @@ SUBROUTINE CERRED( PATH, NUNIT )
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV,
-     $                   CGESDD, CGESVD
+     $                   CGESDD, CGESVD, CGESVDX, CGESVDQ
 *     ..
 *     .. External Functions ..
       LOGICAL            LSAMEN, CSLECT
@@ -495,6 +497,61 @@ SUBROUTINE CERRED( PATH, NUNIT )
          ELSE
             WRITE( NOUT, FMT = 9998 )
          END IF
+*
+*        Test CGESVDQ
+*
+         SRNAMT = 'CGESVDQ'
+         INFOT = 1
+         CALL CGESVDQ( 'X', 'P', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL CGESVDQ( 'A', 'X', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL CGESVDQ( 'A', 'P', 'X', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL CGESVDQ( 'A', 'P', 'T', 'X', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'X', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 6
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', -1, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 7
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', 0, 1, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 9
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 0, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 12
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 -1, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 14
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, -1, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 17
+         CALL CGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, 1, NS, IW, -5, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'CGESVDQ', INFOT, NOUT, LERR, OK )
+         NT = 11
+         IF( OK ) THEN
+            WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
+     $           NT
+         ELSE
+            WRITE( NOUT, FMT = 9998 )
+         END IF
       END IF
 *
 *     Print a summary line.

TESTING/EIG/derred.f

@@ -36,6 +36,8 @@
 *>       DGEJSV   compute SVD of an M-by-N matrix A where M >= N
 *>       DGESVDX  compute SVD of an M-by-N matrix A(by bisection
 *>                and inverse iteration)
+*>       DGESVDQ  compute SVD of an M-by-N matrix A(with a 
+*>                QR-Preconditioned )
 *> \endverbatim
 *
 *  Arguments:
@@ -100,7 +102,7 @@ SUBROUTINE DERRED( PATH, NUNIT )
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           CHKXER, DGEES, DGEESX, DGEEV, DGEEVX, DGEJSV,
-     $                   DGESDD, DGESVD
+     $                   DGESDD, DGESVD, DGESVDX, DGESVQ
 *     ..
 *     .. External Functions ..
       LOGICAL            DSLECT, LSAMEN
@@ -486,6 +488,61 @@ SUBROUTINE DERRED( PATH, NUNIT )
          ELSE
             WRITE( NOUT, FMT = 9998 )
          END IF
+*
+*        Test DGESVDQ
+*
+         SRNAMT = 'DGESVDQ'
+         INFOT = 1
+         CALL DGESVDQ( 'X', 'P', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL DGESVDQ( 'A', 'X', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL DGESVDQ( 'A', 'P', 'X', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL DGESVDQ( 'A', 'P', 'T', 'X', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'X', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 6
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', -1, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 7
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', 0, 1, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 9
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 0, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 12
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 -1, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 14
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, -1, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 17
+         CALL DGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, 1, NS, IW, -5, W, 1, W, 1, INFO )
+         CALL CHKXER( 'DGESVDQ', INFOT, NOUT, LERR, OK )
+         NT = 11
+         IF( OK ) THEN
+            WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
+     $           NT
+         ELSE
+            WRITE( NOUT, FMT = 9998 )
+         END IF
       END IF
 *
 *     Print a summary line.

TESTING/EIG/serred.f

@@ -36,6 +36,8 @@
 *>       SGEJSV   compute SVD of an M-by-N matrix A where M >= N
 *>       SGESVDX  compute SVD of an M-by-N matrix A(by bisection
 *>                and inverse iteration)
+*>       SGESVDQ  compute SVD of an M-by-N matrix A(with a 
+*>                QR-Preconditioned )
 *> \endverbatim
 *
 *  Arguments:
@@ -100,7 +102,7 @@ SUBROUTINE SERRED( PATH, NUNIT )
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           CHKXER, SGEES, SGEESX, SGEEV, SGEEVX, SGEJSV,
-     $                   SGESDD, SGESVD
+     $                   SGESDD, SGESVD, SGESVDX, SGESVDQ
 *     ..
 *     .. External Functions ..
       LOGICAL            SSLECT, LSAMEN
@@ -486,6 +488,61 @@ SUBROUTINE SERRED( PATH, NUNIT )
          ELSE
             WRITE( NOUT, FMT = 9998 )
          END IF
+*
+*        Test SGESVDQ
+*
+         SRNAMT = 'SGESVDQ'
+         INFOT = 1
+         CALL SGESVDQ( 'X', 'P', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL SGESVDQ( 'A', 'X', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL SGESVDQ( 'A', 'P', 'X', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL SGESVDQ( 'A', 'P', 'T', 'X', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'X', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 6
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', -1, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 7
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', 0, 1, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 9
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 0, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 12
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 -1, VT, 0, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 14
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, -1, NS, IW, 1, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 17
+         CALL SGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, 1, NS, IW, -5, W, 1, W, 1, INFO )
+         CALL CHKXER( 'SGESVDQ', INFOT, NOUT, LERR, OK )
+         NT = 11
+         IF( OK ) THEN
+            WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
+     $           NT
+         ELSE
+            WRITE( NOUT, FMT = 9998 )
+         END IF
       END IF
 *
 *     Print a summary line.

TESTING/EIG/zerred.f

@@ -36,6 +36,8 @@
 *>       ZGEJSV   compute SVD of an M-by-N matrix A where M >= N
 *>       ZGESVDX  compute SVD of an M-by-N matrix A(by bisection
 *>                and inverse iteration)
+*>       ZGESVDQ  compute SVD of an M-by-N matrix A(with a 
+*>                QR-Preconditioned )
 *> \endverbatim
 *
 *  Arguments:
@@ -101,7 +103,7 @@ SUBROUTINE ZERRED( PATH, NUNIT )
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           CHKXER, ZGEES, ZGEESX, ZGEEV, ZGEEVX, ZGESVJ,
-     $                   ZGESDD, ZGESVD
+     $                   ZGESDD, ZGESVD, ZGESVDX, ZGESVQ
 *     ..
 *     .. External Functions ..
       LOGICAL            LSAMEN, ZSLECT
@@ -495,6 +497,61 @@ SUBROUTINE ZERRED( PATH, NUNIT )
          ELSE
             WRITE( NOUT, FMT = 9998 )
          END IF
+*
+*        Test ZGESVDQ
+*
+         SRNAMT = 'ZGESVDQ'
+         INFOT = 1
+         CALL ZGESVDQ( 'X', 'P', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL ZGESVDQ( 'A', 'X', 'T', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL ZGESVDQ( 'A', 'P', 'X', 'A', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL ZGESVDQ( 'A', 'P', 'T', 'X', 'A', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'X', 0, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 6
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', -1, 0, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 7
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', 0, 1, A, 1, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 9
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 0, S, U,
+     $                 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 12
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 -1, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 14
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, -1, NS, IW, 1, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         INFOT = 17
+         CALL ZGESVDQ( 'A', 'P', 'T', 'A', 'A', 1, 1, A, 1, S, U,
+     $                 1, VT, 1, NS, IW, -5, W, 1, RW, 1, INFO )
+         CALL CHKXER( 'ZGESVDQ', INFOT, NOUT, LERR, OK )
+         NT = 11
+         IF( OK ) THEN
+            WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
+     $           NT
+         ELSE
+            WRITE( NOUT, FMT = 9998 )
+         END IF
       END IF
 *
 *     Print a summary line.