Re: [Gretl-users] eigensym() issue?
On Fri, 11 Dec 2015, Allin Cottrell wrote:
>> I am not sure if it is my computer or a "bug" ...
or just me (i have a
>> feeling i am doing something terribly wrong..)
>
> At first sight your code looks perfectly normal, and in any case a program
> crash is never really the user's fault.
Absolutely. And it's most unlikely that your computer is at fault.
However, I ran this code on 64-bit Linux under valgrind: no crash, and no
complaints from valgrind about memory management.
Same here. If the problem was really the call to eigensym() the following
(simpler) code would probably trigger the same error. Could you please
check?
<hansl>
cT = 600
loop i=128..144 --quiet
matrix A = mnormal(cT, i)
A = (A .- meanc(A))/cT
matrix AA = A'A
matrix Ve0 = {} #ascending eigenvectors
matrix Va0 = eigensym( AA , &Ve0 )
endloop
</hansl>
-------------------------------------------------------
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
-------------------------------------------------------