intersect - rows - bug

[Levente Torok]

Hi All,

I found a bug in the new intersect implementation. (I did the same in my code, 
too).

octave:63> a
a =

   3   3

octave:64> b
b =

   1   1
   2   2
   3   3
   4   4

octave:65> [I,ili,iri] = intersect( a(:,1), b(:,1) )
I =  3
ili =  1
iri =  3

which is correct, however

octave:66> [I,ili,iri] = intersect_jaroslav( a,b, 'rows' )
I =

   3   3

ili =  1
iri =  2

which is not correct.

The reason is in intersect.m line #71:

    if (nargout > 1)
      ia = ja(ic(ii));                  ## a(ia) == c
      ib = jb(ic(ii+1) - length (a));   ## b(ib) == c      
    endif

Here the "length(a)" doest not work for the case of 'rows'.
In this case it should be replaced by size( a, 1 ).


Levente

-- 
Blogger of http://fapuma.blogspot.com