From 9205713fbc07fa5bcca7d17e74394430762d9aad Mon Sep 17 00:00:00 2001 From: julie Date: Tue, 2 Nov 2010 18:53:38 +0000 Subject: [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 --- SRC/csysv.f | 12 ++---------- 1 file changed, 2 insertions(+), 10 deletions(-) (limited to 'SRC/csysv.f') 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,17 +160,9 @@ 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 * -- cgit v1.2.3