LAPACK Archives

[Lapack] xlanhf div by zero

Hi Ake,
I did not forget you on that one.
It seems that enabling no-ieee, a lot of testing are generating exception.
I am using the -ffpe-trap=zero option of gfortran.
I will get back to you after talking to the team to see what we do about it.
Regards,
Julie


On Oct 22, 2011, at 1:22 PM, Ake Sandgren wrote:

SRC/clanhf.f
line 975
S = S + ( AA / SCALE )**2
SCALE is 0 when N=1
Not sure if this is expected but that S (inf) is then used in the
VALUE = SCALE*SQRT( S )
both of these should cause floating point exception on non-IEEE systems

_______________________________________________
Lapack mailing list
Lapack@Domain.Removed
http://lists.eecs.utk.edu/mailman/listinfo/lapack


<Prev in Thread] Current Thread [Next in Thread>


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or