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

[Nicholas Jankowski]

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

playing with test cases including ndim arrays, noticed that neither octave nor
matlab's quadv can handle something like 

quadv(@(x) cat(3,x,x), 0, 1)

but matlab's integral can and this quadgkv can. i suspect because of reshaping
to a single column for processinng.  that is going to limit any large data
processing to numel=intmax, no?  I don't know if that's something we hould be
concerned about. 

(I started testing and noticed that a trivial integrand got really slow for
n>50M, so I didn't push it to intmax)


    _______________________________________________________

Reply to this item at:

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

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