fix #467

Fix issue ``Wrong declaration of CDUM in c/zgesdd #467`` from @vzecca by
passing A(1,1) instead of A(2,2) on lines 225 and 234 of SRC/[c/z]ungbr.f

Author: langou

SRC/cungbr.f

@@ -222,7 +222,7 @@ SUBROUTINE CUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
                CALL CUNGQR( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE
                IF( M.GT.1 ) THEN
-                  CALL CUNGQR( M-1, M-1, M-1, A( 2, 2 ), LDA, TAU, WORK,
+                  CALL CUNGQR( M-1, M-1, M-1, A( 1, 1 ), LDA, TAU, WORK,
      $                         -1, IINFO )
                END IF
             END IF
@@ -231,7 +231,7 @@ SUBROUTINE CUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
                CALL CUNGLQ( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE
                IF( N.GT.1 ) THEN
-                  CALL CUNGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,
+                  CALL CUNGLQ( N-1, N-1, N-1, A( 1, 1 ), LDA, TAU, WORK,
      $                         -1, IINFO )
                END IF
             END IF

SRC/zungbr.f

@@ -222,7 +222,7 @@ SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
                CALL ZUNGQR( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE
                IF( M.GT.1 ) THEN
-                  CALL ZUNGQR( M-1, M-1, M-1, A( 2, 2 ), LDA, TAU, WORK,
+                  CALL ZUNGQR( M-1, M-1, M-1, A( 1, 1 ), LDA, TAU, WORK,
      $                         -1, IINFO )
                END IF
             END IF
@@ -231,7 +231,7 @@ SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
                CALL ZUNGLQ( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE
                IF( N.GT.1 ) THEN
-                  CALL ZUNGLQ( N-1, N-1, N-1, A( 2, 2 ), LDA, TAU, WORK,
+                  CALL ZUNGLQ( N-1, N-1, N-1, A( 1, 1 ), LDA, TAU, WORK,
      $                         -1, IINFO )
                END IF
             END IF