[SciPy-dev] How to speed up the computation of triple integrals

Nils Wagner nwagner at iam.uni-stuttgart.de
Tue Nov 7 02:25:25 CST 2006


David Huard wrote:
> Have you tried vectorizing ?
>
> You could define
>
> sign_r = array([1,-1,-1,1,1,-1,-1,1])
> sign_s = ...
> sign_t = ...
>
> def h(r,s,t):
>     return .125 * (1 + sign_r * r)  *  (1 + sign_s * s) * (1 + sign_t
> * t)
>
> and similarly for hp.
>
> I don't know how much speed up you'd get but I guess its worth a try.
>
> David
>
> As I understand it, this is not a topic for scipy-dev, but rather for
> scipy-user or numpy-discussion.
>
> 2006/11/6, Nils Wagner <nwagner at iam.uni-stuttgart.de
> <mailto:nwagner at iam.uni-stuttgart.de>>:
>
>     Hi all,
>
>     Is there a way to speed up the computation of triple integrals ?
>
>     Attached is a small script illustrating the task.
>
>     Any pointer would be appreciated.
>
>     Nils
>
>
>
>     _______________________________________________
>     Scipy-dev mailing list
>     Scipy-dev at scipy.org <mailto:Scipy-dev at scipy.org>
>     http://projects.scipy.org/mailman/listinfo/scipy-dev
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev
>   
Thank you for your hint. I will give it a try.

Nils



More information about the Scipy-dev mailing list