[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.9 and 1.15

version 1.9, 2011/11/21 22:19: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 November 2011  *> \date June 2016
 *  *
 *> \ingroup complex16GEsolve  *> \ingroup complex16GEsolve
 *  *
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.4.0) --  *  -- 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..--
 *     November 2011  *     June 2016
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS, RANK        INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS, RANK
Line 285 Line 285
 *              Path 1 - overdetermined or exactly determined  *              Path 1 - overdetermined or exactly determined
 *  *
 *              Compute space needed for ZGEBRD  *              Compute space needed for ZGEBRD
                CALL ZGEBRD( MM, N, A, LDA, S, DUM(1), DUM(1),                 CALL ZGEBRD( MM, N, A, LDA, S, S, DUM(1), DUM(1), DUM(1),
      $                      DUM(1), DUM(1), -1, INFO )       $                      -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', MM, NRHS, N, A, LDA, DUM(1),                 CALL ZUNMBR( 'Q', 'L', 'C', MM, NRHS, N, A, LDA, DUM(1),
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 315 Line 315
      $                -1, INFO )       $                -1, INFO )
                   LWORK_ZGELQF=DUM(1)                    LWORK_ZGELQF=DUM(1)
 *                 Compute space needed for ZGEBRD  *                 Compute space needed for ZGEBRD
                   CALL ZGEBRD( M, M, A, LDA, S, DUM(1), DUM(1),                    CALL ZGEBRD( M, M, A, LDA, S, S, DUM(1), DUM(1),
      $                      DUM(1), 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 346 Line 346
 *                 Path 2 - underdetermined  *                 Path 2 - underdetermined
 *  *
 *                 Compute space needed for ZGEBRD  *                 Compute space needed for ZGEBRD
                   CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1),                    CALL ZGEBRD( M, N, A, LDA, S, S, DUM(1), DUM(1),
      $                      DUM(1), 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.9  
changed lines
  Added in v.1.15

CVSweb interface <joel.bertrand@systella.fr>