[Octave-bug-tracker] [bug #30481] cquad integrator, Octave and Oct-File sources

[David Bateman]

Follow-up Comment #7, bug #30481 (project octave):

But cquad is calling the user function with a vector, while quadpack, by the
way its written is calling the user function with scalar arguments.. If the
user function is vectorrized, this is a huge win for cquad, at least with
Octave. A simple test to show what I mean with an oscillatory function :

octave:1> f = @(x) sin(1 ./ x);
octave:2> tic; quad(f,5e-4,1); toc
Elapsed time is 0.279176 seconds.
octave:3> tic; cquad(f,5e-4,1);toc
Elapsed time is 0.146284 seconds.

D.



    _______________________________________________________

Reply to this item at:

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

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