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

Diff for /rpl/lapack/lapack/zgetrf2.f between versions 1.1 and 1.7

version 1.1, 2015/11/26 11:44:22 version 1.7, 2023/08/07 08:39:20
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/
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       RECURSIVE SUBROUTINE ZGETRF2( M, N, A, LDA, IPIV, INFO )  *       RECURSIVE SUBROUTINE ZGETRF2( M, N, A, LDA, IPIV, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            INFO, LDA, M, N  *       INTEGER            INFO, LDA, M, N
 *       ..  *       ..
Line 17 Line 17
 *       INTEGER            IPIV( * )  *       INTEGER            IPIV( * )
 *       COMPLEX*16         A( LDA, * )  *       COMPLEX*16         A( LDA, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 35 Line 35
 *>  *>
 *> This is the recursive version of the algorithm. It divides  *> This is the recursive version of the algorithm. It divides
 *> the matrix into four submatrices:  *> the matrix into four submatrices:
 *>              *>
 *>        [  A11 | A12  ]  where A11 is n1 by n1 and A22 is n2 by n2  *>        [  A11 | A12  ]  where A11 is n1 by n1 and A22 is n2 by n2
 *>    A = [ -----|----- ]  with n1 = min(m,n)  *>    A = [ -----|----- ]  with n1 = min(m,n)/2
 *>        [  A21 | A22  ]       n2 = n-n1  *>        [  A21 | A22  ]       n2 = n-n1
 *>              *>
 *>                                       [ A11 ]  *>                                       [ A11 ]
 *> The subroutine calls itself to factor [ --- ],  *> The subroutine calls itself to factor [ --- ],
 *>                                       [ A12 ]  *>                                       [ A12 ]
Line 101 Line 101
 *  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 November 2015  
 *  *
 *> \ingroup complex16GEcomputational  *> \ingroup complex16GEcomputational
 *  *
 *  =====================================================================  *  =====================================================================
       RECURSIVE SUBROUTINE ZGETRF2( M, N, A, LDA, IPIV, INFO )        RECURSIVE SUBROUTINE ZGETRF2( M, N, A, LDA, IPIV, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.6.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..--
 *     November 2015  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, M, N        INTEGER            INFO, LDA, M, N
Line 144 Line 141
       EXTERNAL           DLAMCH, IZAMAX        EXTERNAL           DLAMCH, IZAMAX
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           ZGEMM, ZSCAL, ZLASWP, ZTRSM, ZERBLA        EXTERNAL           ZGEMM, ZSCAL, ZLASWP, ZTRSM, XERBLA
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN        INTRINSIC          MAX, MIN
Line 241 Line 238
 *  *
 *        Solve A12  *        Solve A12
 *  *
          CALL ZTRSM( 'L', 'L', 'N', 'U', N1, N2, ONE, A, LDA,            CALL ZTRSM( 'L', 'L', 'N', 'U', N1, N2, ONE, A, LDA,
      $               A( 1, N1+1 ), LDA )       $               A( 1, N1+1 ), LDA )
 *  *
 *        Update A22  *        Update A22
 *  *
          CALL ZGEMM( 'N', 'N', M-N1, N2, N1, -ONE, A( N1+1, 1 ), LDA,            CALL ZGEMM( 'N', 'N', M-N1, N2, N1, -ONE, A( N1+1, 1 ), LDA,
      $               A( 1, N1+1 ), LDA, ONE, A( N1+1, N1+1 ), LDA )       $               A( 1, N1+1 ), LDA, ONE, A( N1+1, N1+1 ), LDA )
 *  *
 *        Factor A22  *        Factor A22
Removed from v.1.1  
changed lines
  Added in v.1.7

CVSweb interface <joel.bertrand@systella.fr>