Re: [Gretl-devel] numerical_hessian

Hello everyone, Sorry for the long delay. As I explained in my first post, I does not have access directly to the code and I need to ask help by email to have a way to reproduce the problem. After applying your fix, ie, replace "for (k=0; k<r-m; k++) {" by "for (k=0; k<r-m-1; k++) {", purify does not find anymore a memory corruption. Nevertheless, there is a new problem at the line 377 "*err = gretl_invert_symmetric_matrix(H);" in gretl_bfgs.c gretl_invert_symmetric_matrix: dpotrf failed with info = 4 (n = 9) To reproduce the problem, you can compile testarma.c and run with the following parameter "a.out EntreesArmasJDD58_precision20.txt 21228" We try this test on an old linux (not reproduce with 1.9.1), sun Solaris ( reproduce with 1.9.1). It is possible there is a bug inside gretl_invert_symmetric_matrix ? It is possible the first fix is not correct ? Is there any requirements on the arma input dataset to work correctly ? Thanks for your help. Best regards, Olivier Gaça 2011/5/5 Olivier Ga <ogaca42(a)gmail.com&gt;