Merge pull request #5910 from martin-frbg/lapack1298

Fix test failures in ?CHKLQ when NB=1 with M=N (Reference-LAPACK PR 1298)

Author: martin-frbg

lapack-netlib/TESTING/LIN/cchklq.f

@@ -394,12 +394,19 @@ SUBROUTINE CCHKLQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL,
                            CALL CLACPY( 'Full', M, N, A, LDA, AF, LDA )
 *
                            SRNAMT = 'CGELS'
-                           CALL CGELS( 'No transpose', M, N, NRHS, AF,
+                            CALL CGELS( 'No transpose', M, N, NRHS, AF,
      $                                 LDA, X, LDA, WORK, LWORK, INFO )
 *
+*                          Re-factorize AF with CGELQF for subsequent
+*                          LQ tests which expect LQ factorization in AF.
+*
+                            CALL CLACPY( 'Full', M, N, A, LDA, AF, LDA )
+                            CALL CGELQF( M, N, AF, LDA, TAU, WORK,
+     $                                   LWORK, INFO )
+*
 *                          Check error code from CGELS.
 *
-                           IF( INFO.NE.0 )
+                            IF( INFO.NE.0 )
      $                        CALL ALAERH( PATH, 'CGELS', INFO, 0, 'N',
      $                                     M, N, NRHS, -1, NB, IMAT,
      $                                     NFAIL, NERRS, NOUT )

lapack-netlib/TESTING/LIN/dchklq.f

@@ -397,12 +397,19 @@ SUBROUTINE DCHKLQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL,
                            CALL DLACPY( 'Full', M, N, A, LDA, AF, LDA )
 *
                            SRNAMT = 'DGELS'
-                           CALL DGELS( 'No transpose', M, N, NRHS, AF,
+                            CALL DGELS( 'No transpose', M, N, NRHS, AF,
      $                                 LDA, X, LDA, WORK, LWORK, INFO )
 *
+*                          Re-factorize AF with DGELQF for subsequent
+*                          LQ tests which expect LQ factorization in AF.
+*
+                            CALL DLACPY( 'Full', M, N, A, LDA, AF, LDA )
+                            CALL DGELQF( M, N, AF, LDA, TAU, WORK,
+     $                                   LWORK, INFO )
+*
 *                          Check error code from DGELS.
 *
-                           IF( INFO.NE.0 )
+                            IF( INFO.NE.0 )
      $                        CALL ALAERH( PATH, 'DGELS', INFO, 0, 'N',
      $                                     M, N, NRHS, -1, NB, IMAT,
      $                                     NFAIL, NERRS, NOUT )

lapack-netlib/TESTING/LIN/schklq.f

@@ -394,12 +394,19 @@ SUBROUTINE SCHKLQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL,
                            CALL SLACPY( 'Full', M, N, A, LDA, AF, LDA )
 *
                            SRNAMT = 'SGELS'
-                           CALL SGELS( 'No transpose', M, N, NRHS, AF,
+                            CALL SGELS( 'No transpose', M, N, NRHS, AF,
      $                                 LDA, X, LDA, WORK, LWORK, INFO )
 *
+*                          Re-factorize AF with SGELQF for subsequent
+*                          LQ tests which expect LQ factorization in AF.
+*
+                            CALL SLACPY( 'Full', M, N, A, LDA, AF, LDA )
+                            CALL SGELQF( M, N, AF, LDA, TAU, WORK,
+     $                                   LWORK, INFO )
+*
 *                          Check error code from SGELS.
 *
-                           IF( INFO.NE.0 )
+                            IF( INFO.NE.0 )
      $                        CALL ALAERH( PATH, 'SGELS', INFO, 0, 'N',
      $                                     M, N, NRHS, -1, NB, IMAT,
      $                                     NFAIL, NERRS, NOUT )

lapack-netlib/TESTING/LIN/zchklq.f

@@ -394,12 +394,19 @@ SUBROUTINE ZCHKLQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL,
                            CALL ZLACPY( 'Full', M, N, A, LDA, AF, LDA )
 *
                            SRNAMT = 'ZGELS'
-                           CALL ZGELS( 'No transpose', M, N, NRHS, AF,
+                            CALL ZGELS( 'No transpose', M, N, NRHS, AF,
      $                                 LDA, X, LDA, WORK, LWORK, INFO )
 *
+*                          Re-factorize AF with ZGELQF for subsequent
+*                          LQ tests which expect LQ factorization in AF.
+*
+                            CALL ZLACPY( 'Full', M, N, A, LDA, AF, LDA )
+                            CALL ZGELQF( M, N, AF, LDA, TAU, WORK,
+     $                                   LWORK, INFO )
+*
 *                          Check error code from ZGELS.
 *
-                           IF( INFO.NE.0 )
+                            IF( INFO.NE.0 )
      $                        CALL ALAERH( PATH, 'ZGELS', INFO, 0, 'N',
      $                                     M, N, NRHS, -1, NB, IMAT,
      $                                     NFAIL, NERRS, NOUT )