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