rpl/lapack/lapack/dla_gerfsx_extended.f - diff

Return to dla_gerfsx_extended.f CVS log

Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dla_gerfsx_extended.f between versions 1.7 and 1.17

version 1.7, 2012/08/22 09:48:15	version 1.17, 2023/08/07 08:38:52
Line 1	Line 1
*> \brief \b DLA_GERFSX_EXTENDED	*> \brief \b DLA_GERFSX_EXTENDED improves the computed solution to a system of linear equations for general matrices by performing extra-precise iterative refinement and provides error bounds and backward error estimates for the solution.
*	*
* =========== 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_GERFSX_EXTENDED + dependencies	*> Download DLA_GERFSX_EXTENDED + dependencies
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">
*> [TGZ]</a>	*> [TGZ]</a>
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">
*> [ZIP]</a>	*> [ZIP]</a>
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">	*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dla_gerfsx_extended.f">
*> [TXT]</a>	*> [TXT]</a>
*> \endhtmlonly	*> \endhtmlonly
*	*
* Definition:	* Definition:
* ===========	* ===========
Line 24	Line 24
* ERRS_N, ERRS_C, RES, AYB, DY,	* ERRS_N, ERRS_C, RES, AYB, DY,
* Y_TAIL, RCOND, ITHRESH, RTHRESH,	* Y_TAIL, RCOND, ITHRESH, RTHRESH,
* DZ_UB, IGNORE_CWISE, INFO )	* DZ_UB, IGNORE_CWISE, 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,
* $ TRANS_TYPE, N_NORMS, ITHRESH	* $ TRANS_TYPE, N_NORMS, ITHRESH
Line 38	Line 38
* DOUBLE PRECISION C( * ), AYB( * ), RCOND, BERR_OUT( * ),	* DOUBLE PRECISION C( * ), AYB( * ), RCOND, BERR_OUT( * ),
* $ ERRS_N( NRHS, * ), ERRS_C( NRHS, * )	* $ ERRS_N( NRHS, * ), ERRS_C( NRHS, * )
* ..	* ..
*	*
*	*
*> \par Purpose:	*> \par Purpose:
* =============	* =============
*>	*>
*> \verbatim	*> \verbatim
*>	*>
*>	*>
*> DLA_GERFSX_EXTENDED improves the computed solution to a system of	*> DLA_GERFSX_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 53	Line 53
*> In addition to normwise error bound, the code provides maximum	*> In addition to normwise error bound, the code provides maximum
*> componentwise error bound if possible. See comments for ERRS_N	*> componentwise error bound if possible. See comments for ERRS_N
*> and ERRS_C for details of the error bounds. Note that this	*> and ERRS_C for details of the error bounds. Note that this
*> subroutine is only resonsible for setting the second fields of	*> subroutine is only responsible for setting the second fields of
*> ERRS_N and ERRS_C.	*> ERRS_N and ERRS_C.
*> \endverbatim	*> \endverbatim
*	*
Line 64	Line 64
*> \verbatim	*> \verbatim
*> PREC_TYPE is INTEGER	*> PREC_TYPE is INTEGER
*> Specifies the intermediate precision to be used in refinement.	*> Specifies the intermediate precision to be used in refinement.
*> The value is defined by ILAPREC(P) where P is a CHARACTER and	*> The value is defined by ILAPREC(P) where P is a CHARACTER and P
*> P = 'S': Single	*> = 'S': Single
*> = 'D': Double	*> = 'D': Double
*> = 'I': Indigenous	*> = 'I': Indigenous
*> = 'X', 'E': Extra	*> = 'X' or 'E': Extra
*> \endverbatim	*> \endverbatim
*>	*>
*> \param[in] TRANS_TYPE	*> \param[in] TRANS_TYPE
*> \verbatim	*> \verbatim
*> TRANS_TYPE is INTEGER	*> TRANS_TYPE is INTEGER
*> Specifies the transposition operation on A.	*> Specifies the transposition operation on A.
*> The value is defined by ILATRANS(T) where T is a CHARACTER and	*> The value is defined by ILATRANS(T) where T is a CHARACTER and T
*> T = 'N': No transpose	*> = 'N': No transpose
*> = 'T': Transpose	*> = 'T': Transpose
*> = 'C': Conjugate transpose	*> = 'C': Conjugate transpose
*> \endverbatim	*> \endverbatim
Line 163	Line 163
*>	*>
*> \param[in,out] Y	*> \param[in,out] Y
*> \verbatim	*> \verbatim
*> Y is DOUBLE PRECISION array, dimension	*> Y is DOUBLE PRECISION array, dimension (LDY,NRHS)
*> (LDY,NRHS)
*> On entry, the solution matrix X, as computed by DGETRS.	*> On entry, the solution matrix X, as computed by DGETRS.
*> On exit, the improved solution matrix Y.	*> On exit, the improved solution matrix Y.
*> \endverbatim	*> \endverbatim
Line 196	Line 195
*>	*>
*> \param[in,out] ERRS_N	*> \param[in,out] ERRS_N
*> \verbatim	*> \verbatim
*> ERRS_N is DOUBLE PRECISION array, dimension	*> ERRS_N is DOUBLE PRECISION array, dimension (NRHS, N_ERR_BNDS)
*> (NRHS, N_ERR_BNDS)
*> For each right-hand side, this array contains information about	*> For each right-hand side, this array contains information about
*> various error bounds and condition numbers corresponding to the	*> various error bounds and condition numbers corresponding to the
*> normwise relative error, which is defined as follows:	*> normwise relative error, which is defined as follows:
Line 243	Line 241
*>	*>
*> \param[in,out] ERRS_C	*> \param[in,out] ERRS_C
*> \verbatim	*> \verbatim
*> ERRS_C is DOUBLE PRECISION array, dimension	*> ERRS_C is DOUBLE PRECISION array, dimension (NRHS, N_ERR_BNDS)
*> (NRHS, N_ERR_BNDS)
*> For each right-hand side, this array contains information about	*> For each right-hand side, this array contains information about
*> various error bounds and condition numbers corresponding to the	*> various error bounds and condition numbers corresponding to the
*> componentwise relative error, which is defined as follows:	*> componentwise relative error, which is defined as follows:
Line 259	Line 256
*> information as described below. There currently are up to three	*> information as described below. There currently are up to three
*> pieces of information returned for each right-hand side. If	*> pieces of information returned for each right-hand side. If
*> componentwise accuracy is not requested (PARAMS(3) = 0.0), then	*> componentwise accuracy is not requested (PARAMS(3) = 0.0), then
*> ERRS_C is not accessed. If N_ERR_BNDS .LT. 3, then at most	*> ERRS_C is not accessed. If N_ERR_BNDS < 3, then at most
*> the first (:,N_ERR_BNDS) entries are returned.	*> the first (:,N_ERR_BNDS) entries are returned.
*>	*>
*> The first index in ERRS_C(i,:) corresponds to the ith	*> The first index in ERRS_C(i,:) corresponds to the ith
Line 358	Line 355
*> DZ_UB is DOUBLE PRECISION	*> DZ_UB is DOUBLE PRECISION
*> Determines when to start considering componentwise convergence.	*> Determines when to start considering componentwise convergence.
*> Componentwise convergence is only considered after each component	*> Componentwise convergence is only considered after each component
*> of the solution Y is stable, which we definte as the relative	*> of the solution Y is stable, which we define as the relative
*> change in each component being less than DZ_UB. The default value	*> change in each component being less than DZ_UB. The default value
*> is 0.25, requiring the first bit to be stable. See LAWN 165 for	*> is 0.25, requiring the first bit to be stable. See LAWN 165 for
*> more details.	*> more details.
Line 382	Line 379
* 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 November 2011
*	*
*> \ingroup doubleGEcomputational	*> \ingroup doubleGEcomputational
*	*
Line 399	Line 394
$ Y_TAIL, RCOND, ITHRESH, RTHRESH,	$ Y_TAIL, RCOND, ITHRESH, RTHRESH,
$ DZ_UB, IGNORE_CWISE, INFO )	$ DZ_UB, IGNORE_CWISE, INFO )
*	*
* -- LAPACK computational routine (version 3.4.0) --	* -- LAPACK computational routine --
* -- 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 2011
*	*
* .. 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 666	Line 660
* op(A) = A, AT, or AH depending on TRANS (and type).	* op(A) = A, AT, or AH depending on TRANS (and type).
*	*
CALL DCOPY( N, B( 1, J ), 1, RES, 1 )	CALL DCOPY( N, B( 1, J ), 1, RES, 1 )
CALL DGEMV( TRANS, N, N, -1.0D+0, A, LDA, Y(1,J), 1, 1.0D+0,	CALL DGEMV( TRANS, N, N, -1.0D+0, A, LDA, Y(1,J), 1, 1.0D+0,
$ RES, 1 )	$ RES, 1 )

DO I = 1, N	DO I = 1, N
Line 685	Line 679
END DO	END DO
*	*
RETURN	RETURN
	*
	* End of DLA_GERFSX_EXTENDED
	*
END	END

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.7
changed lines
	Added in v.1.17