Diff for /rpl/lapack/lapack/zlarft.f between versions 1.11 and 1.17

version 1.11, 2012/07/31 11:06:39 version 1.17, 2016/08/27 15:27:14
Line 1 Line 1
 *> \brief \b ZLARFT  *> \brief \b ZLARFT forms the triangular factor T of a block reflector H = I - vtvH
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
Line 130 Line 130
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date April 2012  *> \date June 2016
 *  *
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
Line 163 Line 163
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT )        SUBROUTINE ZLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT )
 *  *
 *  -- LAPACK auxiliary routine (version 3.4.1) --  *  -- LAPACK auxiliary routine (version 3.6.1) --
 *  -- 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  *     June 2016
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          DIRECT, STOREV        CHARACTER          DIRECT, STOREV
Line 187 Line 187
       INTEGER            I, J, PREVLASTV, LASTV        INTEGER            I, J, PREVLASTV, LASTV
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           ZGEMV, ZLACGV, ZTRMV        EXTERNAL           ZGEMV, ZLACGV, ZTRMV, ZGEMM
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            LSAME        LOGICAL            LSAME
Line 303 Line 303
 *  *
                      CALL ZGEMM( 'N', 'C', K-I, 1, N-K+I-J, -TAU( I ),                       CALL ZGEMM( 'N', 'C', K-I, 1, N-K+I-J, -TAU( I ),
      $                           V( I+1, J ), LDV, V( I, J ), LDV,       $                           V( I+1, J ), LDV, V( I, J ), LDV,
      $                           ONE, T( I+1, I ), LDT )                            $                           ONE, T( I+1, I ), LDT )
                   END IF                    END IF
 *  *
 *                 T(i+1:k,i) := T(i+1:k,i+1:k) * T(i+1:k,i)  *                 T(i+1:k,i) := T(i+1:k,i+1:k) * T(i+1:k,i)

Removed from v.1.11  
changed lines
  Added in v.1.17


CVSweb interface <joel.bertrand@systella.fr>