[BACK]Return to dla_syrfsx_extended.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dla_syrfsx_extended.f between versions 1.12 and 1.13

version 1.12, 2016/08/27 15:34:25 version 1.13, 2017/06/17 10:53:52
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 DLA_SYRFSX_EXTENDED + dependencies   *> Download DLA_SYRFSX_EXTENDED + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dla_syrfsx_extended.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 25 Line 25
 *                                       AYB, DY, Y_TAIL, RCOND, ITHRESH,  *                                       AYB, DY, Y_TAIL, RCOND, ITHRESH,
 *                                       RTHRESH, DZ_UB, IGNORE_CWISE,  *                                       RTHRESH, DZ_UB, IGNORE_CWISE,
 *                                       INFO )  *                                       INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE,  *       INTEGER            INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE,
 *      $                   N_NORMS, ITHRESH  *      $                   N_NORMS, ITHRESH
Line 41 Line 41
 *      $                   ERR_BNDS_NORM( NRHS, * ),  *      $                   ERR_BNDS_NORM( NRHS, * ),
 *      $                   ERR_BNDS_COMP( NRHS, * )  *      $                   ERR_BNDS_COMP( NRHS, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
 *>  *>
 *> \verbatim  *> \verbatim
 *>  *>
 *>   *>
 *> DLA_SYRFSX_EXTENDED improves the computed solution to a system of  *> DLA_SYRFSX_EXTENDED improves the computed solution to a system of
 *> linear equations by performing extra-precise iterative refinement  *> linear equations by performing extra-precise iterative refinement
 *> and provides error bounds and backward error estimates for the solution.  *> and provides error bounds and backward error estimates for the solution.
Line 381 Line 381
 *  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 September 2012  *> \date December 2016
 *  *
 *> \ingroup doubleSYcomputational  *> \ingroup doubleSYcomputational
 *  *
Line 399 Line 399
      $                                RTHRESH, DZ_UB, IGNORE_CWISE,       $                                RTHRESH, DZ_UB, IGNORE_CWISE,
      $                                INFO )       $                                INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.2) --  *  -- LAPACK computational routine (version 3.7.0) --
 *  -- 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..--
 *     September 2012  *     December 2016
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE,        INTEGER            INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE,
Line 545 Line 545
                CALL BLAS_DSYMV2_X(UPLO2, N, -1.0D+0, A, LDA,                 CALL BLAS_DSYMV2_X(UPLO2, N, -1.0D+0, A, LDA,
      $              Y(1, J), Y_TAIL, 1, 1.0D+0, RES, 1, PREC_TYPE)       $              Y(1, J), Y_TAIL, 1, 1.0D+0, RES, 1, PREC_TYPE)
             END IF              END IF
               
 !         XXX: RES is no longer needed.  !         XXX: RES is no longer needed.
             CALL DCOPY( N, RES, 1, DY, 1 )              CALL DCOPY( N, RES, 1, DY, 1 )
             CALL DSYTRS( UPLO, N, 1, AF, LDAF, IPIV, DY, N, INFO )              CALL DSYTRS( UPLO, N, 1, AF, LDAF, IPIV, DY, N, INFO )
Line 557 Line 557
             NORMDX = 0.0D+0              NORMDX = 0.0D+0
             DZ_Z = 0.0D+0              DZ_Z = 0.0D+0
             YMIN = HUGEVAL              YMIN = HUGEVAL
               
             DO I = 1, N              DO I = 1, N
                YK = ABS( Y( I, J ) )                 YK = ABS( Y( I, J ) )
                DYK = ABS( DY( I ) )                 DYK = ABS( DY( I ) )
                  
                IF ( YK .NE. 0.0D+0 ) THEN                 IF ( YK .NE. 0.0D+0 ) THEN
                   DZ_Z = MAX( DZ_Z, DYK / YK )                    DZ_Z = MAX( DZ_Z, DYK / YK )
                ELSE IF ( DYK .NE. 0.0D+0 ) THEN                 ELSE IF ( DYK .NE. 0.0D+0 ) THEN
Line 660 Line 660
             ELSE              ELSE
                CALL DLA_WWADDW( N, Y(1,J), Y_TAIL, DY )                 CALL DLA_WWADDW( N, Y(1,J), Y_TAIL, DY )
             END IF              END IF
               
          END DO           END DO
 *        Target of "IF (Z_STOP .AND. X_STOP)".  Sun's f77 won't EXIT.  *        Target of "IF (Z_STOP .AND. X_STOP)".  Sun's f77 won't EXIT.
  666     CONTINUE   666     CONTINUE
Line 689 Line 689
 *        Compute residual RES = B_s - op(A_s) * Y,  *        Compute residual RES = B_s - op(A_s) * Y,
 *            op(A) = A, A**T, or A**H depending on TRANS (and type).  *            op(A) = A, A**T, or A**H depending on TRANS (and type).
          CALL DCOPY( N, B( 1, J ), 1, RES, 1 )           CALL DCOPY( N, B( 1, J ), 1, RES, 1 )
          CALL DSYMV( UPLO, N, -1.0D+0, A, LDA, Y(1,J), 1, 1.0D+0, RES,            CALL DSYMV( UPLO, N, -1.0D+0, A, LDA, Y(1,J), 1, 1.0D+0, RES,
      $     1 )       $     1 )
            
          DO I = 1, N           DO I = 1, N
             AYB( I ) = ABS( B( I, J ) )              AYB( I ) = ABS( B( I, J ) )
          END DO           END DO
Line 700 Line 700
 *  *
          CALL DLA_SYAMV( UPLO2, N, 1.0D+0,           CALL DLA_SYAMV( UPLO2, N, 1.0D+0,
      $        A, LDA, Y(1, J), 1, 1.0D+0, AYB, 1 )       $        A, LDA, Y(1, J), 1, 1.0D+0, AYB, 1 )
            
          CALL DLA_LIN_BERR( N, N, 1, RES, AYB, BERR_OUT( J ) )           CALL DLA_LIN_BERR( N, N, 1, RES, AYB, BERR_OUT( J ) )
 *  *
 *     End of loop for each RHS.  *     End of loop for each RHS.
Removed from v.1.12  
changed lines
  Added in v.1.13

CVSweb interface <joel.bertrand@systella.fr>