--- rpl/lapack/lapack/dlasd1.f 2010/12/21 13:53:32 1.8 +++ rpl/lapack/lapack/dlasd1.f 2011/07/22 07:38:07 1.9 @@ -1,10 +1,10 @@ SUBROUTINE DLASD1( NL, NR, SQRE, D, ALPHA, BETA, U, LDU, VT, LDVT, $ IDXQ, IWORK, WORK, INFO ) * -* -- LAPACK auxiliary routine (version 3.2.2) -- +* -- LAPACK auxiliary routine (version 3.3.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* June 2010 +* -- April 2011 -- * * .. Scalar Arguments .. INTEGER INFO, LDU, LDVT, NL, NR, SQRE @@ -26,13 +26,13 @@ * * DLASD1 computes the SVD as follows: * -* ( D1(in) 0 0 0 ) -* B = U(in) * ( Z1' a Z2' b ) * VT(in) -* ( 0 0 D2(in) 0 ) +* ( D1(in) 0 0 0 ) +* B = U(in) * ( Z1**T a Z2**T b ) * VT(in) +* ( 0 0 D2(in) 0 ) * * = U(out) * ( D(out) 0) * VT(out) * -* where Z' = (Z1' a Z2' b) = u' VT', and u is a vector of dimension M +* where Z**T = (Z1**T a Z2**T b) = u**T VT**T, and u is a vector of dimension M * with ALPHA and BETA in the NL+1 and NL+2 th entries and zeros * elsewhere; and the entry b is empty if SQRE = 0. * @@ -98,10 +98,10 @@ * * VT (input/output) DOUBLE PRECISION array, dimension(LDVT,M) * where M = N + SQRE. -* On entry VT(1:NL+1, 1:NL+1)' contains the right singular -* vectors of the upper block; VT(NL+2:M, NL+2:M)' contains +* On entry VT(1:NL+1, 1:NL+1)**T contains the right singular +* vectors of the upper block; VT(NL+2:M, NL+2:M)**T contains * the right singular vectors of the lower block. On exit -* VT' contains the right singular vectors of the +* VT**T contains the right singular vectors of the * bidiagonal matrix. * * LDVT (input) INTEGER