To whom it may concern,
As it is posted on the netlib.org webpage, the source for dsyevr.f
reads:
*
=====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0, TWO = 0.0D+0 )
* ..
* .. Local Scalars ..
Clearly, TWO is set to have the value of ZERO above (0.0D+0).
Later, in line 442, the parameter TWO is used:
*
IF (ABSTOL .LE. TWO*N*EPS) THEN
TRYRAC = .TRUE.
ELSE
TRYRAC = .FALSE.
END IF
Why is TWO set to be ZERO? Is this a bug? Is this some intentional
obscurity in the code?
I would greatly appreciate some clarification!
Thanks!
Misha
--
Mikhail Wolfson
Graduate Researcher
Laboratory for Computational Immunology
Department of Chemistry
Massachusetts Institute of Technology
77 Massachusetts Avenue
Room E-19 534
Cambridge, MA 02139
(617) 252-1742
wolfsonm@Domain.Removed
http://web.mit.edu/wolfsonm/
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
http://lists.cs.utk.edu/private/lapack/attachments/20070525/566be666/attachment.htm
|