The LAPACK forum has moved to https://github.com/Reference-LAPACK/lapack/discussions.

Error in ILAENV for lapack 3.1.1

Open discussion regarding features, bugs, issues, vendors, etc.

Error in ILAENV for lapack 3.1.1

Postby whaley » Thu Jun 26, 2008 3:21 pm

It appears to me there is an error in ILAENV where ISPEC=10 and 11 are actually
swapped. The present code says:
*
* ISPEC = 10: ieee NaN arithmetic can be trusted not to trap
*
* ILAENV = 0
ILAENV = 1
IF( ILAENV.EQ.1 ) THEN
ILAENV = IEEECK( 0, 0.0, 1.0 )
END IF
RETURN
*
150 CONTINUE
*
* ISPEC = 11: infinity arithmetic can be trusted not to trap
*
* ILAENV = 0
ILAENV = 1
IF( ILAENV.EQ.1 ) THEN
ILAENV = IEEECK( 1, 0.0, 1.0 )
END IF
RETURN

But if you scope IEEECK, it says:
* ISPEC (input) INTEGER
* Specifies whether to test just for inifinity arithmetic
* or whether to test for infinity and NaN arithmetic.
* = 0: Verify infinity arithmetic only.
* = 1: Verify infinity and NaN arithmetic.


I.e., ISPEC=11 is testing NaNs, whereas ISPEC=10 is testing infinity only,
the opposite of what they are supposed to do.
whaley
 
Posts: 3
Joined: Thu Apr 28, 2005 12:23 pm
Location: University of Texas at San Antonio

Return to User Discussion

Who is online

Users browsing this forum: No registered users and 8 guests