[SciPy-user] slow integrals

Eric Emsellem 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 
routines themselves?
I will of course have a more detailed look at quad_explain

thanks
Eric



More information about the SciPy-user mailing list