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

Diff for /rpl/lapack/lapack/zgelss.f between versions 1.14 and 1.15

version 1.14, 2016/08/27 15:34:45 version 1.15, 2017/06/17 10:54:09
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 ZGELSS + dependencies   *> Download ZGELSS + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgelss.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgelss.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgelss.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgelss.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgelss.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgelss.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE ZGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,  *       SUBROUTINE ZGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,
 *                          WORK, LWORK, RWORK, INFO )  *                          WORK, LWORK, RWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS, RANK  *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS, RANK
 *       DOUBLE PRECISION   RCOND  *       DOUBLE PRECISION   RCOND
Line 29 Line 29
 *       DOUBLE PRECISION   RWORK( * ), S( * )  *       DOUBLE PRECISION   RWORK( * ), S( * )
 *       COMPLEX*16         A( LDA, * ), B( LDB, * ), WORK( * )  *       COMPLEX*16         A( LDA, * ), B( LDB, * ), WORK( * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 165 Line 165
 *  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 June 2016  *> \date June 2016
 *  *
Line 178 Line 178
       SUBROUTINE ZGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,        SUBROUTINE ZGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,
      $                   WORK, LWORK, RWORK, INFO )       $                   WORK, LWORK, RWORK, INFO )
 *  *
 *  -- LAPACK driver routine (version 3.6.1) --  *  -- LAPACK driver routine (version 3.7.0) --
 *  -- 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..--
 *     June 2016  *     June 2016
Line 296 Line 296
                CALL ZUNGBR( 'P', N, N, N, A, LDA, DUM(1),                 CALL ZUNGBR( 'P', N, N, N, A, LDA, DUM(1),
      $                   DUM(1), -1, INFO )       $                   DUM(1), -1, INFO )
                LWORK_ZUNGBR=DUM(1)                 LWORK_ZUNGBR=DUM(1)
 *              Compute total workspace needed   *              Compute total workspace needed
                MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZGEBRD )                 MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZGEBRD )
                MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZUNMBR )                 MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZUNMBR )
                MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZUNGBR )                 MAXWRK = MAX( MAXWRK, 2*N + LWORK_ZUNGBR )
Line 319 Line 319
      $                         DUM(1), -1, INFO )       $                         DUM(1), -1, INFO )
                   LWORK_ZGEBRD=DUM(1)                    LWORK_ZGEBRD=DUM(1)
 *                 Compute space needed for ZUNMBR  *                 Compute space needed for ZUNMBR
                   CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, N, A, LDA,                     CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, N, A, LDA,
      $                DUM(1), B, LDB, DUM(1), -1, INFO )       $                DUM(1), B, LDB, DUM(1), -1, INFO )
                   LWORK_ZUNMBR=DUM(1)                    LWORK_ZUNMBR=DUM(1)
 *                 Compute space needed for ZUNGBR  *                 Compute space needed for ZUNGBR
Line 330 Line 330
                   CALL ZUNMLQ( 'L', 'C', N, NRHS, M, A, LDA, DUM(1),                    CALL ZUNMLQ( 'L', 'C', N, NRHS, M, A, LDA, DUM(1),
      $                 B, LDB, DUM(1), -1, INFO )       $                 B, LDB, DUM(1), -1, INFO )
                   LWORK_ZUNMLQ=DUM(1)                    LWORK_ZUNMLQ=DUM(1)
 *                 Compute total workspace needed   *                 Compute total workspace needed
                   MAXWRK = M + LWORK_ZGELQF                    MAXWRK = M + LWORK_ZGELQF
                   MAXWRK = MAX( MAXWRK, 3*M + M*M + LWORK_ZGEBRD )                    MAXWRK = MAX( MAXWRK, 3*M + M*M + LWORK_ZGEBRD )
                   MAXWRK = MAX( MAXWRK, 3*M + M*M + LWORK_ZUNMBR )                    MAXWRK = MAX( MAXWRK, 3*M + M*M + LWORK_ZUNMBR )
Line 350 Line 350
      $                         DUM(1), -1, INFO )       $                         DUM(1), -1, INFO )
                   LWORK_ZGEBRD=DUM(1)                    LWORK_ZGEBRD=DUM(1)
 *                 Compute space needed for ZUNMBR  *                 Compute space needed for ZUNMBR
                   CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, M, A, LDA,                     CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, M, A, LDA,
      $                DUM(1), B, LDB, DUM(1), -1, INFO )       $                DUM(1), B, LDB, DUM(1), -1, INFO )
                   LWORK_ZUNMBR=DUM(1)                    LWORK_ZUNMBR=DUM(1)
 *                 Compute space needed for ZUNGBR  *                 Compute space needed for ZUNGBR
Removed from v.1.14  
changed lines
  Added in v.1.15

CVSweb interface <joel.bertrand@systella.fr>