[SciPy-User] Double integration with integrate.trapz

[bhanukiran perabathini]

Hi all,

- My task is to perform a double integral (with two free variables)
and make a filled-contour plot of the logarithm of its absolute
square.

- I'm kind of sure that the code has no major mistakes. Nevertheless,
I've badly got to accelerate my double-integration algorithm.

- I bet there must be ways to increase the speed substantially. My
guess is that the over use of scipy.vectorize is the culprit.
   Here I'm attaching my code and the integral, any suggestions on how
to improve the speed would be deeply appreciated.

- I've checked out integrate.dblquad and it doesn't seem to handle
complex types easily. But, no matter what, I primary intension is to
improve on the trapz method.

Thanks a tonne for every second of yours,
- Bhanukiran,
FOSSEE,
IITB, India
-------------- next part --------------
A non-text attachment was scrubbed...
Name: the_integral.pdf
Type: application/pdf
Size: 31627 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20110226/88063950/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dbl_integral.py
Type: text/x-python-script
Size: 1320 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20110226/88063950/attachment.bin>