[Octave-bug-tracker] [bug #46160] image package: normxcorr2 returns complex values sometimes

[Hartmut]

Follow-up Comment #2, bug #46160 (project octave):

I think that I don't like the idea of eliminating small negative values with a
find command, as you suggested, because it costs calculation time without
really doing something. Couldn't we try to avoid the numerically instable
subtraction in line 64 at all? 

Why must there be a minus sign at all in line 64? For 2d images, shouldn't it
be enough to calculate b in the same way that a is calculated one line below?
Something similar to this:


 b = sumsq (b(:)) .* some_size_term;


But since I do not fully understand what the original term does with Nd
images, I hesitate to suggest an easy fix, here.

But if there is no more elegant way to solve this, I should probably accept
the find command instead of keeping the current behavior :(


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?46160>

_______________________________________________
  Nachricht gesendet von/durch Savannah
  http://savannah.gnu.org/