1%----------------------------------------------------------------------
2% File:      xcolor.tst
3% Purpose:   Tests for xcolor package.
4% Author:    A.Kryukov
5% E-address: kryukov@npi.msu.su
6% Vertion:   2.2.1
7% Release:   Aug. 17, 1994
8%----------------------------------------------------------------------
9% Revision:  14/03/91  Start
10%            17/08/94  C0a,C0b,C0c,C0d
11%----------------------------------------------------------------------
12
13%load xcolor$
14%in "xcolor.red"$
15
16% on time$
17
18SUdim n$
19
20        % Set the order of the SU(n) group, i.e. SU(n).
21                % n=3 is a default value.
22
23spTT 1/2$
24
25       % Set the normalization coefficient A:
26		%     Sp(TiTj) = A*Delta(i,j)
27                % A=1/2 is a default value.
28
29                % QG(inQuark,outQuark,gluon) - quark-gluon vertex.
30                % G3(gluon1,gluon2,gluon3)   - three-gluon vertex.
31
32%----------------------------------------------------------------------
33%               Run time (ms) without GC, compiled vertion.
34%----------------------------------------------------------------------
35% Test   | SUN3/60 |i386SX-25|DEC APX   |HP 735    |
36%        |         |RAM=2M   |3000/128MB|          |
37%        |         |         |R=20MB    |R=8MB     |
38%----------------------------------------------------------------------
39% c1     | 323     | 2800    |34        |40        |
40% c2     | 102     | 170     |0         |0         |
41% c3     | 102     | 160     |0         |10        |
42%----------------------------------------------------------------------
43% c4a    | 323     | 1270    |17        |20        |
44% c4b    | 51      | 110     |-         |10        |
45% c4c    | 135     | 160     |17        |0         |
46%----------------------------------------------------------------------
47% c0a    |         |         |17        |20        |
48% c0b    |         |         |0         |0         |
49% c0c    |         |         |17        |20        |
50% c0d    |         |         |17        |0         |
51%----------------------------------------------------------------------
52% cz     | 187     | 990     |0         |10        |
53% c0     | 374     | 330     |          |20        |
54%----------------------------------------------------------------------
55% c91    | 255     | 1480    |17        |10        |
56% c92    | 323     | 1480    |17        |20        |
57% c93    | 2397    | 2970    |68        |130       |
58% c94    | 18411   | 35838   |495       |750       |
59% c95    | 139009  | 262956  |3485      |5670      |
60%----------------------------------------------------------------------
61
62%----------------------------------------------------------------------
63% Test 1. See test run in: A.Kryukov & A.Rodionov
64%               Program "COLOR" for computing the group-theoretical
65%               weight of Feynman diagrams in non-abelian theories.
66%               Comp. Phys. Comm., 48(1988), 327-334.
67%----------------------------------------------------------------------
68
69%----------------------------------------------------------------------
70%         1     4                  1  4                  1
71% c1.   --*--<--*--     c2.   --<--*--*--     c3.   --<--*--<--
72%      |  :     :  |         |     :  :  |         |     :     |
73%      |  :     :  |         |     :  :  |         |     *4    |
74%      |  :     :  |         |   ..:..   |         |    : :    |
75%      |  :     :  |         |  :  :     |         |   :   :   |
76%       --*-->--*--           --*--*-->--           ---*->-*---
77%         2     3               2  3                   2   3
78%----------------------------------------------------------------------
79%
80%         4      2
81%        N  - 2*N  + 1
82% C1 := ---------------
83%             4*N
84%
85%
86%            2
87%         - N  + 1
88% C2 := -----------
89%           4*N
90%
91%
92%            2
93%        N*(N  - 1)
94% C3 := ------------
95%            4
96%
97%----------------------------------------------------------------------
98
99c1 := QG(e4,e1,e5)*QG(e1,e2,e5)*QG(e2,e3,e6)*QG(e3,e4,e6);
100
101
102        4      2
103       n  - 2*n  + 1
104c1 := ---------------
105            4*n
106
107c2 := QG(e4,e1,e5)*QG(e1,e2,e6)*QG(e2,e3,e5)*QG(e3,e4,e6);
108
109
110           2
111        - n  + 1
112c2 := -----------
113          4*n
114
115c3 := QG(e3,e1,e4)*QG(e1,e2,e5)*QG(e2,e3,e6)*G3(e4,e6,e5);
116
117
118           2
119       n*(n  - 1)
120c3 := ------------
121           4
122
123
124%----------------------------------------------------------------------
125%
126% c4a. --*--<--*--      c4b. --*--<--*--      c4c. --*--<--*--
127%     |  :     :  |         |  :     :  |         |  :     :  |
128%     |   ..*..   |         |  :     :  |         |  *...  :  |
129%     |     :     |         |  *.....*  |         |  :   : :  |
130%     |   ..*..   |         |  :     :  |         |  *...:.   |
131%     |  :     :  |         |  :     :  |         |  :   :    |
132%      --*-->--*--           --*-->--*--           --*---*->--
133%
134%----------------------------------------------------------------------
135%
136%          2   2
137%         N *(N  - 1)
138% C4A := -------------
139%              8
140%
141%
142%          2   2
143%         N *(N  - 1)
144% C4B := -------------
145%              8
146%
147%
148% C4C := 0
149%
150%----------------------------------------------------------------------
151
152c4a := QG(e4,e1,e7)*QG(e1,e2,e9)*QG(e2,e3,e8)*QG(e3,e4,e5)
153      *G3(e5,e6,e7)*G3(e6,e8,e9);
154
155
156         2   2
157        n *(n  - 1)
158c4a := -------------
159             8
160
161c4b := c4a;
162
163
164         2   2
165        n *(n  - 1)
166c4b := -------------
167             8
168
169c4c := QG(e4,e1,e7)*QG(e1,e2,e9)*QG(e2,e3,e5)*QG(e3,e4,e8)
170      *G3(e5,e6,e7)*G3(e6,e8,e9);
171
172
173c4c := 0
174
175
176%----------------------------------------------------------------------
177% Test 2.
178%----------------------------------------------------------------------
179%        ......1......
180%       :             :
181%       :    3   6    :
182%      1*.3..*-6-*..4.*2
183%       :    |   |    :  = N^2*(N^2-1)/8
184% C0a:  2    V   |    5
185%       :    |   |    :
186%       :    7   9    :
187%       :    |   |    :
188%        ....*-8-*....
189%            4   5
190%
191%        ......1......
192%       :             :
193%       :    3   6    :
194%      1*.3..*-6-*..4.*2
195%       :    |   |    :  = 0
196% C0b:  2    V   |    5
197%       :    |   |    :
198%       :     \ /     :
199%       :      X      :
200%       :     / \
201%       :    7   9    :
202%       :    |   |    :
203%        ....*-8-*....
204%            4   5
205%
206%-------------------------------------------------------
207
208C0a :=  g3(e1,e3,e2)*g3(e1,e5,e4)
209       *qg(e6,e7,e3)*qg(e7,e8,e2)*qg(e8,e9,e5)*qg(e9,e6,e4);
210
211
212         2   2
213        n *(n  - 1)
214c0a := -------------
215             8
216
217
218C0b :=  g3(e1,e3,e2)*g3(e1,e5,e4)
219       *qg(e6,e9,e3)*qg(e9,e8,e5)*qg(e8,e7,e2)*qg(e7,e6,e4);
220
221
222c0b := 0
223
224
225%--------------------------------------------------------
226%
227%        ......1......
228%       :             :
229%       :    3        :
230%      1*.3..*-6--   -*2
231%       :    |    \ / |  = (N^2-1)/8
232% C0c:  2    V     X  V
233%       :    |    / \ 4
234%       :    |   |   -*6
235%       :    7   9    :
236%       :    |   |    5
237%        ....*-8-*....
238%            4   5
239%
240%        ......1......
241%       :             :
242%       :    3        :
243%      1*.3..*-6------*2
244%       :    |        |  = -(N^2-1)^2/8
245% C0d:  2    V        A
246%       :    |        4
247%       :    |   |----*6
248%       :    7   9    :
249%       :    |   |    5
250%        ....*-8-*....
251%            4   5
252%
253%----------------------------------------------------------------------
254
255C0c := g3(e1,e3,e2)*qg(e9,e4,e1)*qg(e6,e7,e3)
256      *qg(e7,e8,e2)*qg(e8,e9,e5)*qg(e4,e6,e5);
257
258
259         2
260        n  - 1
261c0c := --------
262          8
263
264
265C0d := g3(e1,e3,e2)*qg(e4,e6,e1)*qg(e6,e7,e3)
266      *qg(e7,e8,e2)*qg(e8,e9,e5)*qg(e9,e4,e5);
267
268
269            4      2
270         - n  + 2*n  - 1
271c0d := ------------------
272               8
273
274
275%----------------------------------------------------------------------
276%            1                     1
277% cz.   .....*.....     c0.   .....*.....
278%      :     :     :         :     :e2   :
279%      :     :     :         :e3   *4    :e1
280%      :     :     :         :    : :    :
281%      :     :     :         :   :e4 :e6 :
282%       .....*.....           ...*...*...
283%            2                   2e5 3
284%----------------------------------------------------------------------
285%
286%           2
287% CZ := N*(N  - 1)
288%
289%
290%        2   2
291%       N *(N  - 1)
292% C0 := -----------
293%            2
294%
295%----------------------------------------------------------------------
296
297cz := G3(e1,e2,e3)*G3(e1,e3,e2);
298
299
300          2
301cz := n*(n  - 1)
302
303c0 := G3(e1,e2,e3)*G3(e3,e4,e5)*G3(e5,e6,e1)*G3(e2,e6,e4);
304
305
306        2   2
307       n *(n  - 1)
308c0 := -------------
309            2
310
311
312%$END$
313%----------------------------------------------------------------------
314% Test 3.
315%----------------------------------------------------------------------
316
317%----------------------------------------------------------------------
318%
319%                 1
320% c9n.   .....*.......*.....
321%       :     :   2   :     :
322%       :     *.......*     :
323%       :     :       :     :
324%       :     :  ...  :     :
325%       :     :       :     :
326%       :     :   n   :     :
327%        .....*.......*.....
328%
329%----------------------------------------------------------------------
330%
331% C91 := 0
332%
333%
334%         2   2
335% C92 := N *(N  - 1)
336%
337%
338%          3   2
339%         N *(N  - 1)
340% C93 := -------------
341%              4
342%
343%
344%          2   4       2
345%         N *(N  + 11*N  - 12)
346% C94 := ----------------------
347%                  8
348%
349%
350%          3   4       2
351%         N *(N  + 15*N  - 16)
352% C95 := ----------------------
353%                  16
354%
355%----------------------------------------------------------------------
356
357c91 := G3(e3,e1,e1)*G3(e3,e2,e2);
358
359
360c91 := 0
361
362c92 := G3(e3,e4,e1)*G3(e3,e2,e5)*G3(e6,e1,e4)*G3(e6,e5,e2);
363
364
365        2   2
366c92 := n *(n  - 1)
367
368c93 := G3(e3,e4,e1)*G3(e3,e2,e5)*G3(e6,e7,e4)*G3(e6,e5,e8)
369      *G3(e9,e1,e7)*G3(e9,e8,e2);
370
371
372         3   2
373        n *(n  - 1)
374c93 := -------------
375             4
376
377c94 := G3(e3,e4,e1)*G3(e3,e2,e5)*G3(e6,e7,e4)*G3(e6,e5,e8)
378      *G3(e9,e10,e7)*G3(e9,e8,e11)*G3(e12,e1,e10)*G3(e12,e11,e2);
379
380
381         2   4       2
382        n *(n  + 11*n  - 12)
383c94 := ----------------------
384                 8
385
386c95 := G3(e3,e4,e1)*G3(e3,e2,e5)*G3(e6,e7,e4)*G3(e6,e5,e8)
387      *G3(e9,e10,e7)*G3(e9,e8,e11)*G3(e12,e13,e10)*G3(e12,e11,e14)
388      *G3(e15,e1,e13)*G3(e15,e14,e2);
389
390
391         3   4       2
392        n *(n  + 15*n  - 16)
393c95 := ----------------------
394                 16
395
396
397%----------------------------------------------------------------------
398$
399
400END$
401
402Tested on x86_64-pc-windows CSL
403Time (counter 1): 16 ms
404
405End of Lisp run after 0.01+0.04 seconds
406real 0.22
407user 0.03
408sys 0.06
409