Matrix is the same as in Example 1 for

Function **SELECT** is:

LOGICAL FUNCTION SELECT( WR, WI )

USE LA_PRECISION, ONLY: WPwp

REAL(WP), INTENT(IN) :: WR, WI

INTRINSIC EPSILON, ABS

IF ( (ABS(WR) + ABS(WI) ) 20.0_WP ) THEN

SELECT = .TRUE.

ELSE

SELECT = .FALSE.

END IF

END FUNCTION SELECT

The call:

CALL LA_GEESX( A, WR, WI, SELECT=SELECT, SDIM=SDIM, &

RCONDE=RCONDE, RCONDV=RCONDV )

The eigenvalues of matrix are:

The reciprocal condition number for the average of the eigenvalues is 1.

The reciprocal condition number for the right invariant subspace is 2.21900.

Susan Blackford 2001-08-19