--- rpl/lapack/lapack/zptts2.f 2011/07/22 07:38:20 1.8 +++ rpl/lapack/lapack/zptts2.f 2011/11/21 20:43:20 1.9 @@ -1,9 +1,122 @@ +*> \brief \b ZPTTS2 +* +* =========== DOCUMENTATION =========== +* +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ +* +*> \htmlonly +*> Download ZPTTS2 + dependencies +*> +*> [TGZ] +*> +*> [ZIP] +*> +*> [TXT] +*> \endhtmlonly +* +* Definition: +* =========== +* +* SUBROUTINE ZPTTS2( IUPLO, N, NRHS, D, E, B, LDB ) +* +* .. Scalar Arguments .. +* INTEGER IUPLO, LDB, N, NRHS +* .. +* .. Array Arguments .. +* DOUBLE PRECISION D( * ) +* COMPLEX*16 B( LDB, * ), E( * ) +* .. +* +* +*> \par Purpose: +* ============= +*> +*> \verbatim +*> +*> ZPTTS2 solves a tridiagonal system of the form +*> A * X = B +*> using the factorization A = U**H *D*U or A = L*D*L**H computed by ZPTTRF. +*> D is a diagonal matrix specified in the vector D, U (or L) is a unit +*> bidiagonal matrix whose superdiagonal (subdiagonal) is specified in +*> the vector E, and X and B are N by NRHS matrices. +*> \endverbatim +* +* Arguments: +* ========== +* +*> \param[in] IUPLO +*> \verbatim +*> IUPLO is INTEGER +*> Specifies the form of the factorization and whether the +*> vector E is the superdiagonal of the upper bidiagonal factor +*> U or the subdiagonal of the lower bidiagonal factor L. +*> = 1: A = U**H *D*U, E is the superdiagonal of U +*> = 0: A = L*D*L**H, E is the subdiagonal of L +*> \endverbatim +*> +*> \param[in] N +*> \verbatim +*> N is INTEGER +*> The order of the tridiagonal matrix A. N >= 0. +*> \endverbatim +*> +*> \param[in] NRHS +*> \verbatim +*> NRHS is INTEGER +*> The number of right hand sides, i.e., the number of columns +*> of the matrix B. NRHS >= 0. +*> \endverbatim +*> +*> \param[in] D +*> \verbatim +*> D is DOUBLE PRECISION array, dimension (N) +*> The n diagonal elements of the diagonal matrix D from the +*> factorization A = U**H *D*U or A = L*D*L**H. +*> \endverbatim +*> +*> \param[in] E +*> \verbatim +*> E is COMPLEX*16 array, dimension (N-1) +*> If IUPLO = 1, the (n-1) superdiagonal elements of the unit +*> bidiagonal factor U from the factorization A = U**H*D*U. +*> If IUPLO = 0, the (n-1) subdiagonal elements of the unit +*> bidiagonal factor L from the factorization A = L*D*L**H. +*> \endverbatim +*> +*> \param[in,out] B +*> \verbatim +*> B is DOUBLE PRECISION array, dimension (LDB,NRHS) +*> On entry, the right hand side vectors B for the system of +*> linear equations. +*> On exit, the solution vectors, X. +*> \endverbatim +*> +*> \param[in] LDB +*> \verbatim +*> LDB is INTEGER +*> The leading dimension of the array B. LDB >= max(1,N). +*> \endverbatim +* +* Authors: +* ======== +* +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. +* +*> \date November 2011 +* +*> \ingroup complex16OTHERcomputational +* +* ===================================================================== SUBROUTINE ZPTTS2( IUPLO, N, NRHS, D, E, B, LDB ) * -* -- LAPACK routine (version 3.3.1) -- +* -- LAPACK computational routine (version 3.4.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* -- April 2011 -- +* November 2011 * * .. Scalar Arguments .. INTEGER IUPLO, LDB, N, NRHS @@ -13,51 +126,6 @@ COMPLEX*16 B( LDB, * ), E( * ) * .. * -* Purpose -* ======= -* -* ZPTTS2 solves a tridiagonal system of the form -* A * X = B -* using the factorization A = U**H *D*U or A = L*D*L**H computed by ZPTTRF. -* D is a diagonal matrix specified in the vector D, U (or L) is a unit -* bidiagonal matrix whose superdiagonal (subdiagonal) is specified in -* the vector E, and X and B are N by NRHS matrices. -* -* Arguments -* ========= -* -* IUPLO (input) INTEGER -* Specifies the form of the factorization and whether the -* vector E is the superdiagonal of the upper bidiagonal factor -* U or the subdiagonal of the lower bidiagonal factor L. -* = 1: A = U**H *D*U, E is the superdiagonal of U -* = 0: A = L*D*L**H, E is the subdiagonal of L -* -* N (input) INTEGER -* The order of the tridiagonal matrix A. N >= 0. -* -* NRHS (input) INTEGER -* The number of right hand sides, i.e., the number of columns -* of the matrix B. NRHS >= 0. -* -* D (input) DOUBLE PRECISION array, dimension (N) -* The n diagonal elements of the diagonal matrix D from the -* factorization A = U**H *D*U or A = L*D*L**H. -* -* E (input) COMPLEX*16 array, dimension (N-1) -* If IUPLO = 1, the (n-1) superdiagonal elements of the unit -* bidiagonal factor U from the factorization A = U**H*D*U. -* If IUPLO = 0, the (n-1) subdiagonal elements of the unit -* bidiagonal factor L from the factorization A = L*D*L**H. -* -* B (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS) -* On entry, the right hand side vectors B for the system of -* linear equations. -* On exit, the solution vectors, X. -* -* LDB (input) INTEGER -* The leading dimension of the array B. LDB >= max(1,N). -* * ===================================================================== * * .. Local Scalars ..