[Octave-bug-tracker] [bug #62468] quadgk.m should be able to do array-valued integration

[Nicholas Jankowski]

Follow-up Comment #19, bug #62468 (project octave):

ok, actually looking at _quadgkv_eval_array_ and the two level loop that is
probably slowing it down the most.  is it possible to redefine f(x) = f(x)(:)
before passing to the subfunction such that f(t(:,i,j)) could be a single
vectorized call? or at least one loop level done in a shot? Since there's a
requirement that size(out)=size(in) for arbitrary vector size of in, it seems
like that should be possible? but maybe I'm just not thinking of enough things
that would break it.


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?62468>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/