Hi,

I am posting this to both Octave and Octave-Forge as I think this 
problem encountered while running a function from the 'optim' package
actually exposes some bug in Octave itself. 
If others can confirm this behaviour 'll try to better identify the bug 
and report it to the bug tracker.

I'm using the latest developement version of Octave and I am
trying to use leasqr from the optim package to fit some data 
but I experience some strange behaviour:

------------------------------------
F = @(x, p) p .* exp (- p .* x);
x = linspace (0, 1, 10);
p0 = .1; p = .2;
[f, p1] =  leasqr (x, F(x, p), p0, F)

>> F = @(x, p) p .* exp (- p .* x);
>> x = linspace (0, 1, 10);
>> p0 = .1; p = .2;
>> [f, p1] =  leasqr (x, F(x, p), p0, F)
error: `m' undefined near line 81 column 16
error: evaluating argument list element number 1
error: called from:
error:   /Users/carlo/octave/optim-1.0.16/dfdp.m at line 81, column 7
error:    at line -1, column -1
error:   /Users/carlo/octave/optim-1.0.16/__lm_svd__.m at line 176, column 9
error:   /Users/carlo/octave/optim-1.0.16/leasqr.m at line 583, column 26
error:   ?unknown? at line -1, column -1
>> 
------------------------------------

which is strange for two reasons:

1) there is no line 81 in dfdp.m, I rather believe this refers to line 81 in 
__dfdp__.m:

  prt = zeros (m, n); % initialise Jacobian to Zero

2) m is defined on line 31 of __dfdp__.m and should not be undefined on line 81


but even more strange, after running leasqr, I get

------------------------------------
>> x
error: `x' undefined near line 33 column 1
>> p
error: `p' undefined near line 33 column 1
>> F
error: `F' undefined near line 33 column 1
>> whos
>> 
------------------------------------

So all variables defined before running leasqr have been cleared!
can anyone confirm this behaviour? Any hints how to debug it?

c.






------------------------------------------------------------------------------
All of the data generated in your IT infrastructure is seriously valuable.
Why? It contains a definitive record of application performance, security 
threats, fraudulent activity, and more. Splunk takes this data and makes 
sense of it. IT sense. And common sense.
http://p.sf.net/sfu/splunk-d2d-c2
_______________________________________________
Octave-dev mailing list
Octave-dev@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/octave-dev

Reply via email to