MathGroup Archive 2011

[Date Index] [Thread Index] [Author Index]

Search the Archive

Speeding Up Known Integrations

  • To: mathgroup at
  • Subject: [mg120967] Speeding Up Known Integrations
  • From: Sean Violante <sean.violante07 at>
  • Date: Fri, 19 Aug 2011 06:34:46 -0400 (EDT)
  • Delivered-to:

I have a large collection of iterated integrals to evaluate.
These are basically sums of terms of the form.

  t_1^n1 E^{ k1 t_1}]    t_2^n2 E^{ k2 t_2}  t_3^n3 E^{ k3 t_3}

Is there any way of speeding up this evaluation?  I was wondering 
whether gathering up terms in the above form, and then applying the 
integral to each would work out faster, or is this effectively what   
Integrate is doing anyway?

If faster, how would you suggest gathering up terms?
I have tried a few simple things eg  multiplying by exp[n (k+c)W0t2]
( where n is the largest negative power of Exp[ k t] and Exp [c t] ie so 
expression  becomes polynomial)
and then using CoefficientList [ %,{W0t2, e^k W0t2, e^c W0t2} ]
However, I never quite manage to pull out all the exponential terms

Thanks for your help


below is an example ( where Integrate1 will be replaced by integrate).  
It takes about 100 seconds on my machine.
    1/StdF1[t] \[ExponentialE]^(-t \[Kappa] +
       2 W0t2 \[Kappa]) (VC0 + \[ExponentialE]^(-c W0t2)
          VCC + \[ExponentialE]^(-W0t2 \[Kappa])
          VCK)^2 Integrate1[\[ExponentialE]^(
        W0t1 \[Kappa]) (VC0 + \[ExponentialE]^(-c W0t1)
            VCC + \[ExponentialE]^(-W0t1 \[Kappa])
            VCK) \[Rho] ((\[ExponentialE]^(-t \[Kappa] +
            W0t1 \[Kappa]) (VpC0 + \[ExponentialE]^(-c W0t1)
               VpCC) \[Kappa] Subscript[a, 1] Subscript[\[Xi], 2])/(
          c StdF1[t] - \[Kappa] StdF1[t]) + (\[ExponentialE]^(-c t +
            c W0t1) (VpC0 + \[ExponentialE]^(-c W0t1)
               VpCC) \[Kappa] Subscript[a, 1] Subscript[\[Xi],
           2])/(-c StdF1[t] + \[Kappa] StdF1[
             t]) + (\[ExponentialE]^(-c t +
            c W0t1) (VpC0 + \[ExponentialE]^(-c W0t1) VpCC) Subscript[
           a, 2] Subscript[\[Xi], 2])/StdF1[t]), {W0t1, 0,
        W0t2}] Subscript[a, 1] Subscript[\[Xi], 1], {W0t2, 0, t}] \!
\*SubsuperscriptBox[\(\[Xi]\), \(1\), \(2\)] /.

  • Prev by Date: Re: Getting the Date and Time From a Time Server
  • Next by Date: Re: Inverse Interpolation
  • Previous by thread: Re: Memory when reding in data
  • Next by thread: Re: Speeding Up Known Integrations