[Gretl-devel] Re: different eigenvalue error handling for non-numeric input
On 25/03/2025 15:46, Sven Schreiber wrote:
Hi,
consider this with gretl 2025a:
<console>
? matrix hello = {1/0} # to get invalid input
? = eigensym(hello)
Not a Number geschah bei Berechnung
? = eigen(hello)
1.#INF
? = svd(hello)
nan
</console>
So we have everything: a hard-stopping error, some IEEE code (?), and a
gretl code. Plus, with a 6x6 matrix input I also saw that eigen() spat
out a "Data error".
Ouch
I think it would be good if this could be a little harmonized, no?
Is it even possible to do without breaking backwards compatibility? I we
don't care about that, I'd favour stopping with an error rather than
returning NAs (FWIW, stopping with an error is what both Octave and R do).
-------------------------------------------------------
Riccardo (Jack) Lucchetti
Dipartimento di Scienze Economiche e Sociali (DiSES)
Università Politecnica delle Marche
(formerly known as Università di Ancona)
r.lucchetti(a)univpm.it
http://www2.econ.univpm.it/servizi/hpp/lucchetti
-------------------------------------------------------