diff options
| author | julie <julielangou@users.noreply.github.com> | 2010-11-02 18:53:38 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2010-11-02 18:53:38 +0000 |
| commit | 9205713fbc07fa5bcca7d17e74394430762d9aad (patch) | |
| tree | 8ed1d54491b09b3d768f14003c34b03f96b08d99 /SRC/csysv.f | |
| parent | a4462006fa683e35a6105a0866245ede836387e4 (diff) | |
[xSYTRS/xSYSV] Hide the call to syconv inside trs2 to avoid changing the SYTRS interface.
Update the testing accordingly
[DSYTRI2] Comit dsytri2 to get some feedback
Update the testing accordingly
DSYTRI2 is the Level 3 blas Version of DSYTRI
The actual routine that does the work is DSYTRI2X (name can be changed)
DSYTRI2 is just a wrapper to allow to hide the 2D Workspace required by the routine.
The interface had to be changed to integrate the possibility of doing a workspace query.
DSYTRI2x implementation will be documented in a LAWN. This algorithm was inspired by the following paper:
"Families of Algorithms Related to the Inversion of a Symmetric Positive Definite Matrix"
PAOLO BIENTINESI Duke University and BRIAN GUNTER Delft University of Technology and ROBERT A. VAN DE GEIJN The University of Texas at Austin
Diffstat (limited to 'SRC/csysv.f')
| -rw-r--r-- | SRC/csysv.f | 12 |
1 files changed, 2 insertions, 10 deletions
diff --git a/SRC/csysv.f b/SRC/csysv.f index 0654af7d..1c97a79d 100644 --- a/SRC/csysv.f +++ b/SRC/csysv.f @@ -105,7 +105,7 @@ * * .. Local Scalars .. LOGICAL LQUERY - INTEGER IINFO, LWKOPT, NB + INTEGER LWKOPT, NB * .. * .. External Functions .. LOGICAL LSAME @@ -113,7 +113,7 @@ EXTERNAL ILAENV, LSAME * .. * .. External Subroutines .. - EXTERNAL CSYCONV, CSYTRF, CSYTRS2, XERBLA + EXTERNAL CSYTRF, CSYTRS2, XERBLA * .. * .. Intrinsic Functions .. INTRINSIC MAX @@ -160,18 +160,10 @@ CALL CSYTRF( UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO ) IF( INFO.EQ.0 ) THEN * -* Convert A -* - CALL CSYCONV( UPLO, 'C', N, A, LDA, IPIV, WORK, IINFO ) -* * Solve the system A*X = B, overwriting B with X. * CALL CSYTRS2( UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, INFO ) * -* Revert A -* - CALL CSYCONV( UPLO, 'R', N, A, LDA, IPIV, WORK, IINFO ) -* END IF * WORK( 1 ) = LWKOPT |