[Octave-bug-tracker] [bug #36665] ppval fails when pp.dim > 1 and xi is an array

[Rik]

Follow-up Comment #5, bug #36665 (project octave):

Just a small efficiency question.

In the patch, a column vector is created using the (:) syntax.  This is then
tiled up using repmat to a larger array using repmat ([1, prod(d)]).  Finally,
a transpose of the matrix is performed to turn the multiple column vectors
into row vectors.

Transpose is very fast.  Still, wouldn't it be more efficient to transpose the
smaller object, the column vector, rather than the larger object, the complete
matrix?


Current Code:
dx = (xi - x(idx))(:);
dx = repmat (dx, [1, prod(d)])';

Proposed Code:
dx = (xi - x(idx))(:)';
dx = repmat (dx, [prod(d)], 1);



    _______________________________________________________

Reply to this item at:

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

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