Diff for /rpl/lapack/lapack/dsycon.f between versions 1.3 and 1.9

version 1.3, 2010/08/06 15:28:48 version 1.9, 2011/07/22 07:38:11
Line 1 Line 1
       SUBROUTINE DSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK,        SUBROUTINE DSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK,
      $                   IWORK, INFO )       $                   IWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.2) --  *  -- LAPACK routine (version 3.3.1) --
 *  -- 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..--
 *     November 2006  *  -- April 2011                                                      --
 *  *
 *     Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH.  *     Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH.
 *  *
Line 61 Line 61
 *  *
 *  WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)  *  WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)
 *  *
 *  IWORK    (workspace) INTEGER array, dimension (N)  *  IWORK   (workspace) INTEGER array, dimension (N)
 *  *
 *  INFO    (output) INTEGER  *  INFO    (output) INTEGER
 *          = 0:  successful exit  *          = 0:  successful exit
Line 148 Line 148
       CALL DLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE )        CALL DLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE )
       IF( KASE.NE.0 ) THEN        IF( KASE.NE.0 ) THEN
 *  *
 *        Multiply by inv(L*D*L') or inv(U*D*U').  *        Multiply by inv(L*D*L**T) or inv(U*D*U**T).
 *  *
          CALL DSYTRS( UPLO, N, 1, A, LDA, IPIV, WORK, N, INFO )           CALL DSYTRS( UPLO, N, 1, A, LDA, IPIV, WORK, N, INFO )
          GO TO 30           GO TO 30

Removed from v.1.3  
changed lines
  Added in v.1.9


CVSweb interface <joel.bertrand@systella.fr>