Diff for /rpl/lapack/lapack/zgeqrf.f between versions 1.18 and 1.19

version 1.18, 2020/05/21 21:46:04 version 1.19, 2023/08/07 08:39:19
Line 95 Line 95
 *> \param[in] LWORK  *> \param[in] LWORK
 *> \verbatim  *> \verbatim
 *>          LWORK is INTEGER  *>          LWORK is INTEGER
 *>          The dimension of the array WORK.  LWORK >= max(1,N).  *>          The dimension of the array WORK.
   *>          LWORK >= 1, if MIN(M,N) = 0, and LWORK >= N, otherwise.
 *>          For optimum performance LWORK >= N*NB, where NB is  *>          For optimum performance LWORK >= N*NB, where NB is
 *>          the optimal blocksize.  *>          the optimal blocksize.
 *>  *>
Line 120 Line 121
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date November 2019  
 *  
 *> \ingroup complex16GEcomputational  *> \ingroup complex16GEcomputational
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 145 Line 144
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZGEQRF( M, N, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE ZGEQRF( M, N, A, LDA, TAU, WORK, LWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.9.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 2019  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, LWORK, M, N        INTEGER            INFO, LDA, LWORK, M, N
Line 178 Line 176
 *  *
 *     Test the input arguments  *     Test the input arguments
 *  *
         K = MIN( M, N )
       INFO = 0        INFO = 0
       NB = ILAENV( 1, 'ZGEQRF', ' ', M, N, -1, -1 )        NB = ILAENV( 1, 'ZGEQRF', ' ', M, N, -1, -1 )
       LWKOPT = N*NB  
       WORK( 1 ) = LWKOPT  
       LQUERY = ( LWORK.EQ.-1 )        LQUERY = ( LWORK.EQ.-1 )
       IF( M.LT.0 ) THEN        IF( M.LT.0 ) THEN
          INFO = -1           INFO = -1
Line 189 Line 186
          INFO = -2           INFO = -2
       ELSE IF( LDA.LT.MAX( 1, M ) ) THEN        ELSE IF( LDA.LT.MAX( 1, M ) ) THEN
          INFO = -4           INFO = -4
       ELSE IF( LWORK.LT.MAX( 1, N ) .AND. .NOT.LQUERY ) THEN        ELSE IF( .NOT.LQUERY ) THEN
          INFO = -7           IF( LWORK.LE.0 .OR. ( M.GT.0 .AND. LWORK.LT.MAX( 1, N ) ) )
        $      INFO = -7
       END IF        END IF
       IF( INFO.NE.0 ) THEN        IF( INFO.NE.0 ) THEN
          CALL XERBLA( 'ZGEQRF', -INFO )           CALL XERBLA( 'ZGEQRF', -INFO )
          RETURN           RETURN
       ELSE IF( LQUERY ) THEN        ELSE IF( LQUERY ) THEN
            IF( K.EQ.0 ) THEN
               LWKOPT = 1
            ELSE
               LWKOPT = N*NB
            END IF
            WORK( 1 ) = LWKOPT
          RETURN           RETURN
       END IF        END IF
 *  *
 *     Quick return if possible  *     Quick return if possible
 *  *
       K = MIN( M, N )  
       IF( K.EQ.0 ) THEN        IF( K.EQ.0 ) THEN
          WORK( 1 ) = 1           WORK( 1 ) = 1
          RETURN           RETURN

Removed from v.1.18  
changed lines
  Added in v.1.19


CVSweb interface <joel.bertrand@systella.fr>