Le dimanche 30 décembre 2012 à 16:47 +0400, Happyman a écrit : > Actually > These two functions namely F1 and F2 are really exponential and Bessel > functions respectively. But I can not change its analytic form.. > > I mean is there way to get more quickly the result? > Let's say above mentioned two functions, both of them, one function but the > dimension I showed should not be changed.. > > > What do you think here whether the problem is with 3 dimension or....??? > Thanks in advance for your answer!
The question was: does F1 and F2 change depending on i,j or k ? The answer seems to be yes. I don't think any improvement with vectorisation is possible while using integrate.quad or similar. This function does adaptative meshing of the integration interval. Singular points and/or discontinuities may arise at different points for the several integrands, so that vectorisation would even slow the computation! Maybe you may have a look at romberg's method, which support vector-form. _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion