[SciPy-dev] How to speed up the computation of triple integrals
david.huard at gmail.com
Mon Nov 6 10:28:20 CST 2006
Have you tried vectorizing ?
You could define
sign_r = array([1,-1,-1,1,1,-1,-1,1])
sign_s = ...
sign_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.
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>:
> 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.
> Scipy-dev mailing list
> Scipy-dev at scipy.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Scipy-dev