Tall skinny and short wide routines

1 files changed, 225 insertions, 0 deletions

diff --git a/TESTING/LIN/zerrlqtp.f b/TESTING/LIN/zerrlqtp.f
new file mode 100644
index 00000000..25a079ec
--- /dev/null
+++ b/TESTING/LIN/zerrlqtp.f
@@ -0,0 +1,225 @@
+*> \brief \b ZERRLQTP
+*
+*  =========== DOCUMENTATION ===========
+*
+* Online html documentation available at 
+*            http://www.netlib.org/lapack/explore-html/ 
+*
+*  Definition:
+*  ===========
+*
+*       SUBROUTINE ZERRLQTP( PATH, NUNIT )
+* 
+*       .. Scalar Arguments ..
+*       CHARACTER*3        PATH
+*       INTEGER            NUNIT
+*       ..
+*  
+*
+*> \par Purpose:
+*  =============
+*>
+*> \verbatim
+*>
+*> ZERRLQTP tests the error exits for the complex routines
+*> that use the LQT decomposition of a triangular-pentagonal matrix.
+*> \endverbatim
+*
+*  Arguments:
+*  ==========
+*
+*> \param[in] PATH
+*> \verbatim
+*>          PATH is CHARACTER*3
+*>          The LAPACK path name for the routines to be tested.
+*> \endverbatim
+*>
+*> \param[in] NUNIT
+*> \verbatim
+*>          NUNIT is INTEGER
+*>          The unit number for output.
+*> \endverbatim
+*
+*  Authors:
+*  ========
+*
+*> \author Univ. of Tennessee 
+*> \author Univ. of California Berkeley 
+*> \author Univ. of Colorado Denver 
+*> \author NAG Ltd. 
+*
+*> \date November 2011
+*
+*> \ingroup double_lin
+*
+*  =====================================================================
+      SUBROUTINE ZERRLQTP( PATH, NUNIT )
+      IMPLICIT NONE
+*
+*  -- LAPACK test routine (version 3.4.0) --
+*  -- 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 ..
+      CHARACTER*3        PATH
+      INTEGER            NUNIT
+*     ..
+*
+*  =====================================================================
+*
+*     .. Parameters ..
+      INTEGER            NMAX
+      PARAMETER          ( NMAX = 2 )
+*     ..
+*     .. Local Scalars ..
+      INTEGER            I, INFO, J
+*     ..
+*     .. Local Arrays ..
+      COMPLEX*16         A( NMAX, NMAX ), T( NMAX, NMAX ), W( NMAX ),
+     $                   B( NMAX, NMAX ), C( NMAX, NMAX )
+*     ..
+*     .. External Subroutines ..
+      EXTERNAL           ALAESM, CHKXER, ZTPLQT2, ZTPLQT,
+     $                   ZTPMLQT 
+*     ..
+*     .. Scalars in Common ..
+      LOGICAL            LERR, OK
+      CHARACTER*32       SRNAMT
+      INTEGER            INFOT, NOUT
+*     ..
+*     .. Common blocks ..
+      COMMON             / INFOC / INFOT, NOUT, OK, LERR
+      COMMON             / SRNAMC / SRNAMT
+*     ..
+*     .. Intrinsic Functions ..
+      INTRINSIC          DBLE, DCMPLX
+*     ..
+*     .. Executable Statements ..
+*
+      NOUT = NUNIT
+      WRITE( NOUT, FMT = * )
+*
+*     Set the variables to innocuous values.
+*
+      DO J = 1, NMAX
+         DO I = 1, NMAX
+            A( I, J ) = 1.D0 / DCMPLX( DBLE( I+J ), 0.D0 )
+            C( I, J ) = 1.D0 / DCMPLX( DBLE( I+J ), 0.D0 )
+            T( I, J ) = 1.D0 / DCMPLX( DBLE( I+J ), 0.D0 )
+         END DO
+         W( J ) = 0.0
+      END DO
+      OK = .TRUE.
+*
+*     Error exits for TPLQT factorization
+*
+*     ZTPLQT
+*
+      SRNAMT = 'ZTPLQT'
+      INFOT = 1
+      CALL ZTPLQT( -1, 1, 0, 1, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 2
+      CALL ZTPLQT( 1, -1, 0, 1, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 3
+      CALL ZTPLQT( 0, 1, -1, 1, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 3
+      CALL ZTPLQT( 0, 1, 1, 1, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 4
+      CALL ZTPLQT( 0, 1, 0, 0, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 4
+      CALL ZTPLQT( 1, 1, 0, 2, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 6
+      CALL ZTPLQT( 2, 1, 0, 2, A, 1, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 8
+      CALL ZTPLQT( 2, 1, 0, 1, A, 2, B, 1, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 10
+      CALL ZTPLQT( 2, 2, 1, 2, A, 2, B, 2, T, 1, W, INFO )
+      CALL CHKXER( 'ZTPLQT', INFOT, NOUT, LERR, OK )
+*
+*     ZTPLQT2
+*
+      SRNAMT = 'ZTPLQT2'
+      INFOT = 1
+      CALL ZTPLQT2( -1, 0, 0, A, 1, B, 1, T, 1, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+      INFOT = 2
+      CALL ZTPLQT2( 0, -1, 0, A, 1, B, 1, T, 1, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+      INFOT = 3
+      CALL ZTPLQT2( 0, 0, -1, A, 1, B, 1, T, 1, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+      INFOT = 5
+      CALL ZTPLQT2( 2, 2, 0, A, 1, B, 2, T, 2, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+      INFOT = 7
+      CALL ZTPLQT2( 2, 2, 0, A, 2, B, 1, T, 2, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+      INFOT = 9
+      CALL ZTPLQT2( 2, 2, 0, A, 2, B, 2, T, 1, INFO )
+      CALL CHKXER( 'ZTPLQT2', INFOT, NOUT, LERR, OK )
+*
+*     ZTPMLQT
+*
+      SRNAMT = 'ZTPMLQT'
+      INFOT = 1
+      CALL ZTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 2
+      CALL ZTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 3
+      CALL ZTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 4
+      CALL ZTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 5
+      CALL ZTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      INFOT = 6
+      CALL ZTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 7
+      CALL ZTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 9
+      CALL ZTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 11
+      CALL ZTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 13
+      CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+      INFOT = 15
+      CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, 
+     $              W, INFO )
+      CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK )
+*
+*     Print a summary line.
+*
+      CALL ALAESM( PATH, OK, NOUT )
+*
+      RETURN
+*
+*     End of ZERRLQT
+*
+      END