1:24 p.m.
Dear all,
I've got a function which involves matrix inversion and is frequently
called in loop. From time to time matrix inversion, however, fails. Thus
I tried to use the catch command similar to this tiny example:
<hansl>
function matrix matI(matrix mat)
return inv(mat)
end function
mat = zeros(2,2)
scalar err = 0
catch mI = matI(mat)
err = $error
if err==0
printf "WORKS \n"
else
printf "FAILS \n"
endif
</hansl>
In this example-script everything works fine. However, in my larger
script involving many different things before function "matI" is called,
gretl stops giving me the error message:
<gretl-output>
Unmatched "endif"
</gretl-output>
and the corresponding terminal output (using latest CVS on ubuntu) is:
<terminal>
gretl_invert_symmetric_matrix: dpotrf failed with info = 2 (n = 2)
dsycon: rcond = 5.47239e-18
gretl_clipboard_get: info = 0 (TARGET_UTF8_STRING)
gretl_clipboard_get: info = 0 (TARGET_UTF8_STRING)
</terminal>
Unfortunately I am not able to reproduce this error in a smaller example
script.
Artur
1:41 p.m.
Dear Artur,
Does your code involve something like the following lines?
if erro
quit
endif
Best,
Henrique
1:47 p.m.
No, actually not, Henrique :-/
Artur
Am 06.04.2016 um 15:41 schrieb Henrique Andrade:
Dear Artur,
Does your code involve something like the following lines?
if erro
quit
endif
Best,
Henrique
_______________________________________________
Gretl-devel mailing list
Gretl-devel(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-devel
6:04 p.m.
On Wed, 6 Apr 2016, Artur T. wrote:
Dear all,
I've got a function which involves matrix inversion and is frequently
called in loop. From time to time matrix inversion, however, fails. Thus
I tried to use the catch command similar to this tiny example:
<hansl>
function matrix matI(matrix mat)
return inv(mat)
end function
mat = zeros(2,2)
scalar err = 0
catch mI = matI(mat)
err = $error
if err==0
printf "WORKS \n"
else
printf "FAILS \n"
endif
</hansl>
In this example-script everything works fine. However, in my larger
script involving many different things before function "matI" is called,
gretl stops giving me the error message:
<gretl-output>
Unmatched "endif"
</gretl-output>
Executing the script above with current gretl one sees, inter alia:
<output>
? catch mI = matI(mat)
In regard to function matI:
Warning: "catch" should not be used on calls to user-defined functions
</output>
This should be taken seriously. For "catch" to work reliably it must
be applied directly to a built-in command or function. On return from
a user-defined function it's too late to ensure integrity of gretl's
state in case of error, so (although you may be OK in simple cases)
problems such as you mention are to be expected.
You should move the "catch" into the function where inv() is called.
Here's one possible approach:
<hansl>
function matrix try_invert (matrix m, scalar *err)
catch matrix ret = inv(m)
err = $error
if err
ret = {}
endif
return ret
end function
mat = zeros(2,2)
scalar err = 0
mI = try_invert(mat, &err)
if !err
printf "WORKS\n"
else
printf "FAILS\n"
endif
</hansl>
And here's another:
<hansl>
function scalar try_invert (matrix m, matrix *minv)
catch minv = inv(m)
return $error
end function
mat = zeros(2,2)
matrix minv
err = try_invert(mat, &minv)
if !err
printf "WORKS\n"
else
printf "FAILS\n"
endif
</hansl>
Allin
11:12 a.m.
That was very helpful, Allin! Thank you for this.
This problem is solved now.
Best,
Artur
Am 07.04.2016 um 20:04 schrieb Allin Cottrell:
<output>
? catch mI = matI(mat)
In regard to function matI:
Warning: "catch" should not be used on calls to user-defined functions
</output>
This should be taken seriously. For "catch" to work reliably it must be
applied directly to a built-in command or function. On return from a
user-defined function it's too late to ensure integrity of gretl's state
in case of error, so (although you may be OK in simple cases) problems
such as you mention are to be expected.
You should move the "catch" into the function where inv() is called.
Here's one possible approach:
<hansl>
function matrix try_invert (matrix m, scalar *err)
catch matrix ret = inv(m)
err = $error
if err
ret = {}
endif
return ret
end function
mat = zeros(2,2)
scalar err = 0
mI = try_invert(mat, &err)
if !err
printf "WORKS\n"
else
printf "FAILS\n"
endif
</hansl>
And here's another:
<hansl>
function scalar try_invert (matrix m, matrix *minv)
catch minv = inv(m)
return $error
end function
mat = zeros(2,2)
matrix minv
err = try_invert(mat, &minv)
if !err
printf "WORKS\n"
else
printf "FAILS\n"
endif
</hansl>
Allin
_______________________________________________
Gretl-devel mailing list
Gretl-devel(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-devel
3661
days inactive
3663
days old
4 comments
3 participants
participants (3)
-
Allin Cottrell -
Artur T. -
Henrique Andrade