version 1.13, 2016/08/27 15:34:38
|
version 1.19, 2023/08/07 08:39:06
|
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 DSPOSV + dependencies |
*> Download DSPOSV + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsposv.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsposv.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsposv.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsposv.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsposv.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsposv.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* SUBROUTINE DSPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK, |
* SUBROUTINE DSPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK, |
* SWORK, ITER, INFO ) |
* SWORK, ITER, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* CHARACTER UPLO |
* CHARACTER UPLO |
* INTEGER INFO, ITER, LDA, LDB, LDX, N, NRHS |
* INTEGER INFO, ITER, LDA, LDB, LDX, N, NRHS |
Line 30
|
Line 30
|
* DOUBLE PRECISION A( LDA, * ), B( LDB, * ), WORK( N, * ), |
* DOUBLE PRECISION A( LDA, * ), B( LDB, * ), WORK( N, * ), |
* $ X( LDX, * ) |
* $ X( LDX, * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 106
|
Line 106
|
*> triangular part of the matrix A, and the strictly upper |
*> triangular part of the matrix A, and the strictly upper |
*> triangular part of A is not referenced. |
*> triangular part of A is not referenced. |
*> On exit, if iterative refinement has been successfully used |
*> On exit, if iterative refinement has been successfully used |
*> (INFO.EQ.0 and ITER.GE.0, see description below), then A is |
*> (INFO = 0 and ITER >= 0, see description below), then A is |
*> unchanged, if double precision factorization has been used |
*> unchanged, if double precision factorization has been used |
*> (INFO.EQ.0 and ITER.LT.0, see description below), then the |
*> (INFO = 0 and ITER < 0, see description below), then the |
*> array A contains the factor U or L from the Cholesky |
*> array A contains the factor U or L from the Cholesky |
*> factorization A = U**T*U or A = L*L**T. |
*> factorization A = U**T*U or A = L*L**T. |
*> \endverbatim |
*> \endverbatim |
Line 186
|
Line 186
|
* 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 |
|
* |
* |
*> \ingroup doublePOsolve |
*> \ingroup doublePOsolve |
* |
* |
Line 199
|
Line 197
|
SUBROUTINE DSPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK, |
SUBROUTINE DSPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK, |
$ SWORK, ITER, INFO ) |
$ SWORK, ITER, INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.6.1) -- |
* -- LAPACK driver 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..-- |
* June 2016 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 235
|
Line 232
|
* |
* |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL DAXPY, DSYMM, DLACPY, DLAT2S, DLAG2S, SLAG2D, |
EXTERNAL DAXPY, DSYMM, DLACPY, DLAT2S, DLAG2S, SLAG2D, |
$ SPOTRF, SPOTRS, XERBLA |
$ SPOTRF, SPOTRS, DPOTRF, DPOTRS, XERBLA |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
INTEGER IDAMAX |
INTEGER IDAMAX |
Line 413
|
Line 410
|
30 CONTINUE |
30 CONTINUE |
* |
* |
* If we are at this place of the code, this is because we have |
* If we are at this place of the code, this is because we have |
* performed ITER=ITERMAX iterations and never satisified the |
* performed ITER=ITERMAX iterations and never satisfied the |
* stopping criterion, set up the ITER flag accordingly and follow |
* stopping criterion, set up the ITER flag accordingly and follow |
* up on double precision routine. |
* up on double precision routine. |
* |
* |
Line 434
|
Line 431
|
* |
* |
RETURN |
RETURN |
* |
* |
* End of DSPOSV. |
* End of DSPOSV |
* |
* |
END |
END |