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 /TESTING/LIN/ddrvsyx.f | |
| parent | a4462006fa683e35a6105a0866245ede836387e4 (diff) | |
| download | lapack-9205713fbc07fa5bcca7d17e74394430762d9aad.tar.gz lapack-9205713fbc07fa5bcca7d17e74394430762d9aad.tar.bz2 lapack-9205713fbc07fa5bcca7d17e74394430762d9aad.zip | |
[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 'TESTING/LIN/ddrvsyx.f')
| -rw-r--r-- | TESTING/LIN/ddrvsyx.f | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/TESTING/LIN/ddrvsyx.f b/TESTING/LIN/ddrvsyx.f index 11b1116f..fca0f372 100644 --- a/TESTING/LIN/ddrvsyx.f +++ b/TESTING/LIN/ddrvsyx.f @@ -112,7 +112,7 @@ * .. External Subroutines .. EXTERNAL ALADHD, ALAERH, ALASVM, DERRVX, DGET04, DLACPY, $ DLARHS, DLASET, DLATB4, DLATMS, DPOT02, DPOT05, - $ DSYSV, DSYSVX, DSYT01, DSYTRF, DSYTRI, XLAENV, + $ DSYSV, DSYSVX, DSYT01, DSYTRF, DSYTRI2, XLAENV, $ DSYSVXX * .. * .. Scalars in Common .. @@ -301,8 +301,9 @@ * Compute inv(A) and take its norm. * CALL DLACPY( UPLO, N, N, AFAC, LDA, AINV, LDA ) - CALL DSYTRI( UPLO, N, AINV, LDA, IWORK, WORK, - $ INFO ) + LWORK = (N+NB+1)*(NB+3) + CALL DSYTRI2( UPLO, N, AINV, LDA, IWORK, WORK, + $ LWORK, INFO ) AINVNM = DLANSY( '1', UPLO, N, AINV, LDA, RWORK ) * * Compute the 1-norm condition number of A. |