1 DOUBLE PRECISION FUNCTION DGET06( RCOND, RCONDC )
2 *
3 * -- LAPACK test routine (version 3.1) --
4 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
5 * November 2006
6 *
7 * .. Scalar Arguments ..
8 DOUBLE PRECISION RCOND, RCONDC
9 * ..
10 *
11 * Purpose
12 * =======
13 *
14 * DGET06 computes a test ratio to compare two values for RCOND.
15 *
16 * Arguments
17 * ==========
18 *
19 * RCOND (input) DOUBLE PRECISION
20 * The estimate of the reciprocal of the condition number of A,
21 * as computed by DGECON.
22 *
23 * RCONDC (input) DOUBLE PRECISION
24 * The reciprocal of the condition number of A, computed as
25 * ( 1/norm(A) ) / norm(inv(A)).
26 *
27 * =====================================================================
28 *
29 * .. Parameters ..
30 DOUBLE PRECISION ZERO, ONE
31 PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
32 * ..
33 * .. Local Scalars ..
34 DOUBLE PRECISION EPS, RAT
35 * ..
36 * .. External Functions ..
37 DOUBLE PRECISION DLAMCH
38 EXTERNAL DLAMCH
39 * ..
40 * .. Intrinsic Functions ..
41 INTRINSIC MAX, MIN
42 * ..
43 * .. Executable Statements ..
44 *
45 EPS = DLAMCH( 'Epsilon' )
46 IF( RCOND.GT.ZERO ) THEN
47 IF( RCONDC.GT.ZERO ) THEN
48 RAT = MAX( RCOND, RCONDC ) / MIN( RCOND, RCONDC ) -
49 $ ( ONE-EPS )
50 ELSE
51 RAT = RCOND / EPS
52 END IF
53 ELSE
54 IF( RCONDC.GT.ZERO ) THEN
55 RAT = RCONDC / EPS
56 ELSE
57 RAT = ZERO
58 END IF
59 END IF
60 DGET06 = RAT
61 RETURN
62 *
63 * End of DGET06
64 *
65 END
2 *
3 * -- LAPACK test routine (version 3.1) --
4 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
5 * November 2006
6 *
7 * .. Scalar Arguments ..
8 DOUBLE PRECISION RCOND, RCONDC
9 * ..
10 *
11 * Purpose
12 * =======
13 *
14 * DGET06 computes a test ratio to compare two values for RCOND.
15 *
16 * Arguments
17 * ==========
18 *
19 * RCOND (input) DOUBLE PRECISION
20 * The estimate of the reciprocal of the condition number of A,
21 * as computed by DGECON.
22 *
23 * RCONDC (input) DOUBLE PRECISION
24 * The reciprocal of the condition number of A, computed as
25 * ( 1/norm(A) ) / norm(inv(A)).
26 *
27 * =====================================================================
28 *
29 * .. Parameters ..
30 DOUBLE PRECISION ZERO, ONE
31 PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
32 * ..
33 * .. Local Scalars ..
34 DOUBLE PRECISION EPS, RAT
35 * ..
36 * .. External Functions ..
37 DOUBLE PRECISION DLAMCH
38 EXTERNAL DLAMCH
39 * ..
40 * .. Intrinsic Functions ..
41 INTRINSIC MAX, MIN
42 * ..
43 * .. Executable Statements ..
44 *
45 EPS = DLAMCH( 'Epsilon' )
46 IF( RCOND.GT.ZERO ) THEN
47 IF( RCONDC.GT.ZERO ) THEN
48 RAT = MAX( RCOND, RCONDC ) / MIN( RCOND, RCONDC ) -
49 $ ( ONE-EPS )
50 ELSE
51 RAT = RCOND / EPS
52 END IF
53 ELSE
54 IF( RCONDC.GT.ZERO ) THEN
55 RAT = RCONDC / EPS
56 ELSE
57 RAT = ZERO
58 END IF
59 END IF
60 DGET06 = RAT
61 RETURN
62 *
63 * End of DGET06
64 *
65 END