Re: Speeding up simple Mathematica expressions?
- To: mathgroup at smc.vnet.net
- Subject: [mg63254] Re: Speeding up simple Mathematica expressions?
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Tue, 20 Dec 2005 23:35:38 -0500 (EST)
- Organization: The University of Western Australia
- References: <do8ioc$rvd$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <do8ioc$rvd$1 at smc.vnet.net>, AES <siegman at stanford.edu>
wrote:
> I'd appreciate some practical advice on speeding up some simple function
> evaluations.
>
> I'm evaluating a series of functions of which a typical example is
>
> f[a_, x_] := Sum[
> Exp[-(Pi a)^2 n^2 -
> ((x - n Sqrt[1 - (Pi^2 a^4)])/a)^2],
> {n, -Infinity, Infinity}];
Essentially, this question appeared in the In and Out column of The
Mathematica Journal 9(4) under the heading "Sum of Gaussians". After
completing the square of the exponent, such sums can be expressed in
closed form in terms of EllipticTheta functions. See
http://functions.wolfram.com/09.03.06.0001.01
For this particular example,
f[a_, x_] = a Sqrt[Pi] E^(-(Pi a x)^2) *
EllipticTheta[3, -Pi Sqrt[1 - a^4 Pi^2] x, E^(-(Pi a)^2)]
> (The function is essentially a set of narrow gaussian peaks located at x
> ? n Sqrt[1 - (Pi a^2)^2] ? n , with the peak amplitudes dropping off
> rapidly with increasing x.)
>
> Despite being a fairly simple function, this evaluates very slowly on my
> iBook G4 -- takes a long time to make a plot of say f[0.1, x] for 0 <
> x < 3. What can or should I do to speed this up?
Using the EllipticTheta form, evaluation is immediate.
> c) Since the individual peaks have very little overlap for a < 0.2,
> maybe I can truncate the series to a small range of n?
Indeed, this will also work fine for such parameter values.
Cheers,
Paul
_______________________________________________________________________
Paul Abbott Phone: 61 8 6488 2734
School of Physics, M013 Fax: +61 8 6488 1014
The University of Western Australia (CRICOS Provider No 00126G)
AUSTRALIA http://physics.uwa.edu.au/~paul