[Octave-bug-tracker] [bug #58790] [MXE Octave] Windows 32bit build closes unexpectedly while running the test suite

[Markus Mützel]

Follow-up Comment #11, bug #58790 (project octave):

I was able to condense the test case down to the following:

# from rgb2ind's BISTs
rgb = rand (1000, 1000, 3);
[ind, map] = rgb2ind (rgb);

# from gmres's BISTs
dim = 100;
A = spdiags ([[1./(2:2:2*(dim-1)) 0]; 1./(1:2:2*dim-1); ...
[0 1./(2:2:2*(dim-1))]]', -1:1, dim, dim);
A = A'*A;
b = (1:dim).'/dim;
[x, flag] = gmres (@(x) A*x, b, dim, 1e-10, 1e6,...
                   @(x) diag (diag (A)) \ x, [], []);


Sometimes, I have to execute it twice. But most of the time, Octave closes on
first execution.

If I remove the call to rgb2ind, I can run the code more than 100 times
without the GUI closing.

The following is basically what rgb2ind is doing and it also triggers the
bug:

# from rgb2ind
rgb = rand (1000, 1000, 3);
R = rgb(:,:,1,:);
G = rgb(:,:,2,:);
B = rgb(:,:,3,:);
x = reshape (1:numel (R), size (R));

map    = unique ([R(:) G(:) B(:)], "rows");
[~, x] = ismember ([R(:) G(:) B(:)], map, "rows");

# from gmres's BISTs
dim = 100;
A = spdiags ([[1./(2:2:2*(dim-1)) 0]; 1./(1:2:2*dim-1); ...
[0 1./(2:2:2*(dim-1))]]', -1:1, dim, dim);
A = A'*A;
b = (1:dim).'/dim;
[x, flag] = gmres (@(x) A*x, b, dim, 1e-10, 1e6,...
                   @(x) diag (diag (A)) \ x, [], []);


All variables involved are dense double matrices. I'm not sure where to look
next.

    _______________________________________________________

Reply to this item at:

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

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