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

Diff for /rpl/lapack/lapack/zgebal.f between versions 1.10 and 1.19

version 1.10, 2011/11/21 22:19:44 version 1.19, 2018/05/29 06:55:22
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 ZGEBAL + dependencies   *> Download ZGEBAL + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgebal.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgebal.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgebal.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgebal.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgebal.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgebal.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )  *       SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          JOB  *       CHARACTER          JOB
 *       INTEGER            IHI, ILO, INFO, LDA, N  *       INTEGER            IHI, ILO, INFO, LDA, N
Line 28 Line 28
 *       DOUBLE PRECISION   SCALE( * )  *       DOUBLE PRECISION   SCALE( * )
 *       COMPLEX*16         A( LDA, * )  *       COMPLEX*16         A( LDA, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 83 Line 83
 *>  *>
 *> \param[out] ILO  *> \param[out] ILO
 *> \verbatim  *> \verbatim
   *>          ILO is INTEGER
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[out] IHI  *> \param[out] IHI
 *> \verbatim  *> \verbatim
   *>          IHI is INTEGER
 *>          ILO and IHI are set to INTEGER such that on exit  *>          ILO and IHI are set to INTEGER such that on exit
 *>          A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N.  *>          A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N.
 *>          If JOB = 'N' or 'S', ILO = 1 and IHI = N.  *>          If JOB = 'N' or 'S', ILO = 1 and IHI = N.
Line 116 Line 118
 *  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 2011  *> \date June 2017
 *  *
 *> \ingroup complex16GEcomputational  *> \ingroup complex16GEcomputational
 *  *
Line 160 Line 162
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )        SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.0) --  *  -- LAPACK computational routine (version 3.7.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..--
 *     November 2011  *     June 2017
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOB        CHARACTER          JOB
Line 189 Line 191
       INTEGER            I, ICA, IEXC, IRA, J, K, L, M        INTEGER            I, ICA, IEXC, IRA, J, K, L, M
       DOUBLE PRECISION   C, CA, F, G, R, RA, S, SFMAX1, SFMAX2, SFMIN1,        DOUBLE PRECISION   C, CA, F, G, R, RA, S, SFMAX1, SFMAX2, SFMIN1,
      $                   SFMIN2       $                   SFMIN2
       COMPLEX*16         CDUM  
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       LOGICAL            DISNAN, LSAME        LOGICAL            DISNAN, LSAME
       INTEGER            IZAMAX        INTEGER            IZAMAX
       DOUBLE PRECISION   DLAMCH        DOUBLE PRECISION   DLAMCH, DZNRM2
       EXTERNAL           DISNAN, LSAME, IZAMAX, DLAMCH        EXTERNAL           DISNAN, LSAME, IZAMAX, DLAMCH, DZNRM2
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZDSCAL, ZSWAP        EXTERNAL           XERBLA, ZDSCAL, ZSWAP
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, DIMAG, MAX, MIN        INTRINSIC          ABS, DBLE, DIMAG, MAX, MIN
 *     ..  
 *     .. Statement Functions ..  
       DOUBLE PRECISION   CABS1  
 *     ..  
 *     .. Statement Function definitions ..  
       CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( DIMAG( CDUM ) )  
 *     ..  
 *     .. Executable Statements ..  
 *  *
 *     Test the input parameters  *     Test the input parameters
 *  *
Line 324 Line 317
       NOCONV = .FALSE.        NOCONV = .FALSE.
 *  *
       DO 200 I = K, L        DO 200 I = K, L
          C = ZERO  
          R = ZERO  
 *  *
          DO 150 J = K, L           C = DZNRM2( L-K+1, A( K, I ), 1 )
             IF( J.EQ.I )           R = DZNRM2( L-K+1, A( I, K ), LDA )
      $         GO TO 150  
             C = C + CABS1( A( J, I ) )  
             R = R + CABS1( A( I, J ) )  
   150    CONTINUE  
          ICA = IZAMAX( L, A( 1, I ), 1 )           ICA = IZAMAX( L, A( 1, I ), 1 )
          CA = ABS( A( ICA, I ) )           CA = ABS( A( ICA, I ) )
          IRA = IZAMAX( N-K+1, A( I, K ), LDA )           IRA = IZAMAX( N-K+1, A( I, K ), LDA )
Removed from v.1.10  
changed lines
  Added in v.1.19

CVSweb interface <joel.bertrand@systella.fr>