Re: Why does Limit sometimes just return what I input?
- To: mathgroup at smc.vnet.net
- Subject: [mg44605] Re: Why does Limit sometimes just return what I input?
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Mon, 17 Nov 2003 03:38:40 -0500 (EST)
- Organization: The University of Western Australia
- References: <bp28u9$351$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <bp28u9$351$1 at smc.vnet.net>, jmonk <jmonk at hep.man.ac.uk>
wrote:
> I am trying to find the limit of some gamma matrix traces when a small
> parameter \[Epsilon]->0. I am using the package Tracer.m
Note that I picked up an error/omission in this package when I loaded it
in. In the SortLine routine there is the following:
If[ list === rlist, ,
If[ Sort[list] === Sort[rlist],
Now, since the (extended) syntax of If is
If[condition, t, f, u] gives u if condition evaluates to neither True
nor False.
then the meaning associated with the repeated comma is rather important.
I note that the pakcage at
http://library.wolfram.com/infocenter/MathSource/2987/
has a revision date of 1991.
> and I define 4 lines:
>
> Spur[l1, l2];
> line11 = G[l1, p1, p3, G5, n, {ν}, (p1 + k), {μ}];
> line12 = G[l1, p1, p3, G5, n, {μ}, (p3 - k), {ν}];
> line21 = G[l2, p2, p4, G5, np, {λ}, (p2 - k), {μ}];
> line22 = G[l2, p2, p4, G5, np, {μ}, (p4 + k), {λ}];
Your message came through with some very strange characters. Perhaps you
could re-post in a readable format?
> The vectors n and np are well defined by
>
> OnShell[on,{n, xx_,\[Chi]1. xx + λ\_1**\[Chi2p.xx]},
> {np, xx_, lim[ζ1 . xx + λ\_\[Zeta]2p.¹*¶xx}];
>
> with \[Chi]1, \[Chi]2p, ζ1and\[Zeta]2p similarly defined. The momenta
> p1...p4 also have their dot products defined. n and np depend on a small
> parameter \[Epsilon], and it turns out that the highest order terms in the
> lines line11...line22 are of order 1/\[Epsilon].
>
> I want to take the limit e.g.
>
> Limit[line11*\[Epsilon],\[Epsilon]->0,Direction->-1]
Instead of taking the limit (using Limit), it is usually preferable to
perform a series expansion using Series.
> and then contract the 2 lines on the \[Mu] Lorentz index, but usually when
> I try to take the limit it just returns the limit unevaluated.
> Sometimes, however, it does return the limit. I can use the function
> Residue to get the answer, but I want to know why Limit does not work.
Well, Residue (or Series) may be more appropriate anyway. In my
experience, Limit is not "as solid" as Series.
> The expressions for the lines are quite long - there are 338 terms, so I
> wonder if it is too long for mathematica to cope with unless all of the
> computer resources are available to it. I am running Mathematica 5.0 on
> Mac Os X, G4 1Ghz 256 Mb.
>
> The full script I have is too long to include, but I could forward it if
> requested.
Perhaps you could post the relevant input expressions or give a URL to
where the Notebook can be obtained.
Cheers,
Paul
--
Paul Abbott Phone: +61 8 9380 2734
School of Physics, M013 Fax: +61 8 9380 1014
The University of Western Australia (CRICOS Provider No 00126G)
35 Stirling Highway
Crawley WA 6009 mailto:paul at physics.uwa.edu.au
AUSTRALIA http://physics.uwa.edu.au/~paul