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

[Michael Leitner]

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

Nicholas, that seems correct, however this made me see some other problems:
first, if waypoints are given, they would have to be reversed also. This would
probably make quite ugly code at this point, but I ask myself: is it necessary
at all to have always a<b?  Apart from line 193, it seems to me that all the
rest can equally well handle reversed subintervals (perhaps the isinf(a)
(analogously isinf(b)) transformation implicitly assumes that if isinf(a) it
is minus infinity, this would have to be looked into). A further problem: line
192 should test also whether a and b are real, not only waypoints. 

And h0 is always only set and never used. This can definitely go.

And two documentation bugs: only when no waypoints and the limits are real are
given, the interval is divided into 10 subintervals, and if the limits or any
waypoint is complex, it is a piecewise linear (or straight) path, not only
continuous.


    _______________________________________________________

Reply to this item at:

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

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