[Octave-bug-tracker] [bug #60620] [octave forge] (optim) 'expfit' function from optim-1.6.1 package terminates with segmentation fault

[anonymous]

Follow-up Comment #9, bug #60620 (project octave):

I've added a number of checkpoints to source code of 'expfit.m'. The piece of
code of interest:

"
fprintf(stderr, "expfit CHECKPOINT 1\n");
% Solve for polynomial coefficients
A = hankel( y(1:N-deg),y(N-deg:N) );
fprintf(stderr, "expfit CHECKPOINT 2\n");
s = - A(:,1:deg) \ A(:,deg+1);
fprintf(stderr, "expfit CHECKPOINT 3\n");
".

This is what I get during execution:
"
octave:2> expfit(16, 0, 1, rand(1, 32))
expfit CHECKPOINT 1
expfit CHECKPOINT 2
Segmentation fault (core dumped)
".

So the original code line producing segfault is:

"
s = - A(:,1:deg) \ A(:,deg+1);
".

IMO it's some kind of problem with underlying linear algebra library.

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?60620>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/