| Ilya Zakharevich on Sun, 1 Oct 2000 23:00:03 -0400 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| eigen() strikes again |
With a week-old version:
? mateigen([1,2;3,4])
*** precision too low in eigen.
? \p 38
realprecision = 38 significant digits
? mateigen([1,2;3,4])
*** precision too low in eigen.
? \p80
realprecision = 86 significant digits (80 digits displayed)
? mateigen([1,2;3,4])
%1 =
[-1.45742710...
Looking at the code for the Gauss algorithm (in 2.0.15, which would
just silently give a wrong answer), I would expect that the correct
message should be not "precision too low", but "some random collection
of coincidences did not happen"... Right?
Ilya