[Octave-bug-tracker] [bug #39314] convn(x, y, "valid") produces results that are different from Matlab

[Carnë Draug]

Follow-up Comment #14, bug #39314 (project octave):

I'm not sure that it is possible to account for the "order" problem the way
N-dimensional convolution is being done. convn reduces the problem to a series
of 2 dimensional convolutions, whose results are added.

You will also notice that that there's no precision errors when limited to 2
dimensions.

I see now that I have missed the markup on my previous comment. The example
was:


octave> a = rand (1, 27);
octave> b = rand (1, 27);
octave> c1 = sum (a .* b);
octave> c2 = sum (a(end:-1:1) .* b(end:-1:1));
octave> c1 - c2
ans =  1.7764e-15
-verbatim

    _______________________________________________________

Reply to this item at:

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

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/