[SciPy-user] slow integrals
emsellem at obs.univ-lyon1.fr
Thu Jun 8 02:55:33 CDT 2006
> What do you use for your quadrature routine in your C version?
I am using NAG routines. d01amf for the +/- infinity one and a direct
Gaussian quadrature for the inner integral.
> Looking at quadrature(), though, it looks like it wraps the passed function
> in an inefficient value-at-a-time wrapper (vec_func). Try removing that
> wrapper from the source, so it will call your function with an array of
> points to evaluate. Or raise a TypeError in your function if passed a scalar.
not sure I see what you mean here. Do you suggest I modify the scipy
I will of course have a more detailed look at quad_explain
More information about the SciPy-user