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

Diff for /rpl/lapack/lapack/dorcsd.f between versions 1.10 and 1.16

version 1.10, 2014/01/27 09:28:24 version 1.16, 2023/08/07 08:39:01
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 DORCSD + dependencies   *> Download DORCSD + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dorcsd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dorcsd.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dorcsd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dorcsd.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dorcsd.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dorcsd.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 23 Line 23
 *                                    LDX12, X21, LDX21, X22, LDX22, THETA,  *                                    LDX12, X21, LDX21, X22, LDX22, THETA,
 *                                    U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,  *                                    U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,
 *                                    LDV2T, WORK, LWORK, IWORK, INFO )  *                                    LDV2T, WORK, LWORK, IWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS  *       CHARACTER          JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS
 *       INTEGER            INFO, LDU1, LDU2, LDV1T, LDV2T, LDX11, LDX12,  *       INTEGER            INFO, LDU1, LDU2, LDV1T, LDV2T, LDX11, LDX12,
Line 37 Line 37
 *      $                   X12( LDX12, * ), X21( LDX21, * ), X22( LDX22,  *      $                   X12( LDX12, * ), X21( LDX21, * ), X22( LDX22,
 *      $                   * )  *      $                   * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 186 Line 186
 *>  *>
 *> \param[out] U1  *> \param[out] U1
 *> \verbatim  *> \verbatim
 *>          U1 is DOUBLE PRECISION array, dimension (P)  *>          U1 is DOUBLE PRECISION array, dimension (LDU1,P)
 *>          If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.  *>          If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 199 Line 199
 *>  *>
 *> \param[out] U2  *> \param[out] U2
 *> \verbatim  *> \verbatim
 *>          U2 is DOUBLE PRECISION array, dimension (M-P)  *>          U2 is DOUBLE PRECISION array, dimension (LDU2,M-P)
 *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal  *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal
 *>          matrix U2.  *>          matrix U2.
 *> \endverbatim  *> \endverbatim
Line 213 Line 213
 *>  *>
 *> \param[out] V1T  *> \param[out] V1T
 *> \verbatim  *> \verbatim
 *>          V1T is DOUBLE PRECISION array, dimension (Q)  *>          V1T is DOUBLE PRECISION array, dimension (LDV1T,Q)
 *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal  *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal
 *>          matrix V1**T.  *>          matrix V1**T.
 *> \endverbatim  *> \endverbatim
Line 227 Line 227
 *>  *>
 *> \param[out] V2T  *> \param[out] V2T
 *> \verbatim  *> \verbatim
 *>          V2T is DOUBLE PRECISION array, dimension (M-Q)  *>          V2T is DOUBLE PRECISION array, dimension (LDV2T,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal  *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
 *>          matrix V2**T.  *>          matrix V2**T.
 *> \endverbatim  *> \endverbatim
Line 284 Line 284
 *  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 2013  
 *  *
 *> \ingroup doubleOTHERcomputational  *> \ingroup doubleOTHERcomputational
 *  *
Line 300 Line 298
      $                             U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,       $                             U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,
      $                             LDV2T, WORK, LWORK, IWORK, INFO )       $                             LDV2T, WORK, LWORK, IWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.5.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 2013  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS        CHARACTER          JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS
Line 340 Line 337
      $                   WANTV1T, WANTV2T       $                   WANTV1T, WANTV2T
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DBBCSD, DLACPY, DLAPMR, DLAPMT, DLASCL, DLASET,        EXTERNAL           DBBCSD, DLACPY, DLAPMR, DLAPMT,
      $                   DORBDB, DORGLQ, DORGQR, XERBLA       $                   DORBDB, DORGLQ, DORGQR, XERBLA
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
Line 464 Line 461
          IB22D = IB21E + MAX( 1, Q - 1 )           IB22D = IB21E + MAX( 1, Q - 1 )
          IB22E = IB22D + MAX( 1, Q )           IB22E = IB22D + MAX( 1, Q )
          IBBCSD = IB22E + MAX( 1, Q - 1 )           IBBCSD = IB22E + MAX( 1, Q - 1 )
          CALL DBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,            CALL DBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
      $                THETA, THETA, U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,       $                THETA, THETA, U1, LDU1, U2, LDU2, V1T, LDV1T, V2T,
      $                LDV2T, U1, U1, U1, U1, U1, U1, U1, U1, WORK, -1,       $                LDV2T, U1, U1, U1, U1, U1, U1, U1, U1, WORK, -1,
      $                CHILDINFO )       $                CHILDINFO )
Line 579 Line 576
 *     Permute rows and columns to place identity submatrices in top-  *     Permute rows and columns to place identity submatrices in top-
 *     left corner of (1,1)-block and/or bottom-right corner of (1,2)-  *     left corner of (1,1)-block and/or bottom-right corner of (1,2)-
 *     block and/or bottom-right corner of (2,1)-block and/or top-left  *     block and/or bottom-right corner of (2,1)-block and/or top-left
 *     corner of (2,2)-block   *     corner of (2,2)-block
 *  *
       IF( Q .GT. 0 .AND. WANTU2 ) THEN        IF( Q .GT. 0 .AND. WANTU2 ) THEN
          DO I = 1, Q           DO I = 1, Q
Removed from v.1.10  
changed lines
  Added in v.1.16

CVSweb interface <joel.bertrand@systella.fr>