[BACK]Return to zlabrd.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/zlabrd.f between versions 1.13 and 1.20

version 1.13, 2012/12/14 14:22:49 version 1.20, 2023/08/07 08:39:28
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 ZLABRD + dependencies   *> Download ZLABRD + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zlabrd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zlabrd.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zlabrd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zlabrd.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zlabrd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zlabrd.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,  *       SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,
 *                          LDY )  *                          LDY )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            LDA, LDX, LDY, M, N, NB  *       INTEGER            LDA, LDX, LDY, M, N, NB
 *       ..  *       ..
Line 29 Line 29
 *       COMPLEX*16         A( LDA, * ), TAUP( * ), TAUQ( * ), X( LDX, * ),  *       COMPLEX*16         A( LDA, * ), TAUP( * ), TAUQ( * ), X( LDX, * ),
 *      $                   Y( LDY, * )  *      $                   Y( LDY, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 111 Line 111
 *>  *>
 *> \param[out] TAUQ  *> \param[out] TAUQ
 *> \verbatim  *> \verbatim
 *>          TAUQ is COMPLEX*16 array dimension (NB)  *>          TAUQ is COMPLEX*16 array, dimension (NB)
 *>          The scalar factors of the elementary reflectors which  *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.  *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim  *> \endverbatim
Line 152 Line 152
 *  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 September 2012  
 *  *
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
Line 212 Line 210
       SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,        SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,
      $                   LDY )       $                   LDY )
 *  *
 *  -- LAPACK auxiliary routine (version 3.4.2) --  *  -- LAPACK auxiliary 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..--
 *     September 2012  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            LDA, LDX, LDY, M, N, NB        INTEGER            LDA, LDX, LDY, M, N, NB
Line 270 Line 267
             ALPHA = A( I, I )              ALPHA = A( I, I )
             CALL ZLARFG( M-I+1, ALPHA, A( MIN( I+1, M ), I ), 1,              CALL ZLARFG( M-I+1, ALPHA, A( MIN( I+1, M ), I ), 1,
      $                   TAUQ( I ) )       $                   TAUQ( I ) )
             D( I ) = ALPHA              D( I ) = DBLE( ALPHA )
             IF( I.LT.N ) THEN              IF( I.LT.N ) THEN
                A( I, I ) = ONE                 A( I, I ) = ONE
 *  *
Line 310 Line 307
                ALPHA = A( I, I+1 )                 ALPHA = A( I, I+1 )
                CALL ZLARFG( N-I, ALPHA, A( I, MIN( I+2, N ) ), LDA,                 CALL ZLARFG( N-I, ALPHA, A( I, MIN( I+2, N ) ), LDA,
      $                      TAUP( I ) )       $                      TAUP( I ) )
                E( I ) = ALPHA                 E( I ) = DBLE( ALPHA )
                A( I, I+1 ) = ONE                 A( I, I+1 ) = ONE
 *  *
 *              Compute X(i+1:m,i)  *              Compute X(i+1:m,i)
Line 354 Line 351
             ALPHA = A( I, I )              ALPHA = A( I, I )
             CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,              CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                   TAUP( I ) )       $                   TAUP( I ) )
             D( I ) = ALPHA              D( I ) = DBLE( ALPHA )
             IF( I.LT.M ) THEN              IF( I.LT.M ) THEN
                A( I, I ) = ONE                 A( I, I ) = ONE
 *  *
Line 388 Line 385
                ALPHA = A( I+1, I )                 ALPHA = A( I+1, I )
                CALL ZLARFG( M-I, ALPHA, A( MIN( I+2, M ), I ), 1,                 CALL ZLARFG( M-I, ALPHA, A( MIN( I+2, M ), I ), 1,
      $                      TAUQ( I ) )       $                      TAUQ( I ) )
                E( I ) = ALPHA                 E( I ) = DBLE( ALPHA )
                A( I+1, I ) = ONE                 A( I+1, I ) = ONE
 *  *
 *              Compute Y(i+1:n,i)  *              Compute Y(i+1:n,i)
Removed from v.1.13  
changed lines
  Added in v.1.20

CVSweb interface <joel.bertrand@systella.fr>