Diff for /rpl/lapack/lapack/zungbr.f between versions 1.18 and 1.19

version 1.18, 2018/05/29 07:18:41 version 1.19, 2023/08/07 08:39:43
Line 150 Line 150
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date April 2012  
 *  
 *> \ingroup complex16GBcomputational  *> \ingroup complex16GBcomputational
 *  *
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.7.0) --  *  -- LAPACK computational routine --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     April 2012  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          VECT        CHARACTER          VECT
Line 222 Line 219
                CALL ZUNGQR( M, N, K, A, LDA, TAU, WORK, -1, IINFO )                 CALL ZUNGQR( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE              ELSE
                IF( M.GT.1 ) THEN                 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, LDA, TAU, WORK, -1,
      $                         -1, IINFO )       $                         IINFO )
                END IF                 END IF
             END IF              END IF
          ELSE           ELSE
Line 231 Line 228
                CALL ZUNGLQ( M, N, K, A, LDA, TAU, WORK, -1, IINFO )                 CALL ZUNGLQ( M, N, K, A, LDA, TAU, WORK, -1, IINFO )
             ELSE              ELSE
                IF( N.GT.1 ) THEN                 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, LDA, TAU, WORK, -1,
      $                         -1, IINFO )       $                         IINFO )
                END IF                 END IF
             END IF              END IF
          END IF           END IF
          LWKOPT = WORK( 1 )           LWKOPT = INT( DBLE( WORK( 1 ) ) )
          LWKOPT = MAX (LWKOPT, MN)           LWKOPT = MAX (LWKOPT, MN)
       END IF        END IF
 *  *

Removed from v.1.18  
changed lines
  Added in v.1.19


CVSweb interface <joel.bertrand@systella.fr>