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

version 1.15, 2016/08/27 15:35:13 version 1.19, 2023/08/07 08:39:43
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *> \htmlonly  *> \htmlonly
 *> Download ZUNGBR + dependencies   *> Download ZUNGBR + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zungbr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zungbr.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zungbr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zungbr.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zungbr.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zungbr.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )  *       SUBROUTINE ZUNGBR( VECT, M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          VECT  *       CHARACTER          VECT
 *       INTEGER            INFO, K, LDA, LWORK, M, N  *       INTEGER            INFO, K, LDA, LWORK, M, N
Line 27 Line 27
 *       .. Array Arguments ..  *       .. Array Arguments ..
 *       COMPLEX*16         A( LDA, * ), TAU( * ), WORK( * )  *       COMPLEX*16         A( LDA, * ), TAU( * ), WORK( * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 145 Line 145
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \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.4.1) --  *  -- 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 183 Line 180
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            LSAME        LOGICAL            LSAME
       INTEGER            ILAENV        EXTERNAL           LSAME
       EXTERNAL           LSAME, ILAENV  
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZUNGLQ, ZUNGQR        EXTERNAL           XERBLA, ZUNGLQ, ZUNGQR
Line 223 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 232 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.15  
changed lines
  Added in v.1.19


CVSweb interface <joel.bertrand@systella.fr>