Re: [Matlisp-users] Re: Problem with function geev

[Paul L. I. <aea...@gm...>]

My email probably seemed somewhat strange coming after yours, since we
came to similar conclusions.  My first email was delayed by more than
a day(for moderation purposes?) and you verified what I was asking for
before it could arrive :)

As I think we have both found, just about any change in the code
before that loop causes dlamch to behave correctly.  I chose to set
B=3D2 because it's a trivial improvement to the LAPACK code(it
eliminates an unnecessary loop iteration).

I have witnessed the error on two computers, one with a celeron cpu,
and the other with a pentium III.  In both cases I am using cmucl-19b
on FreeBSD 5.4.  My f77 on both is 3.4.2.

This bug is annoying.  I hope this helps.

Paul