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

[Carnë Draug]

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

I think the cause of the different values comes from the order the values are
added. Looking at the code for inner convolution and comparing to other, it
seems that it's being done in reverse order so machine precision errors get
in. Kind of like in the following example:


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/