BOSS 7.1.0
BESIII Offline Software System
Loading...
Searching...
No Matches
eemmg-lib-new/src/pjfry.cpp
Go to the documentation of this file.
1/*
2 * pjfry.cpp - interface functions
3 *
4 * this file is part of PJFry library
5 * Copyright 2011 Valery Yundin
6 */
7
8#include "common.h"
9#include "cache.h"
10#include "kinem.h"
11
12#include "pjfry.h"
13
14/* ------------------------------------------------------------
15 * ------------------------------------------------------------
16 * PENTAGON Integral section
17 * ------------------------------------------------------------
18 * ------------------------------------------------------------
19*/
20
21// Masses shifted to the left to comply with LoopTools notation
22
23pj_complex PJFry::E0v0(const double p1, const double p2, const double p3, const double p4, const double p5,
24 const double s12, const double s23, const double s34, const double s45, const double s15,
25 const double m5, const double m1, const double m2, const double m3, const double m4)
26{
27 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
28 return ICache::getE(0, kinem);
29}
30
31pj_complex PJFry::E0v0(const double p1, const double p2, const double p3, const double p4, const double p5,
32 const double s12, const double s23, const double s34, const double s45, const double s15,
33 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
34{
35 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
36 return ICache::getE(abs(ep), kinem);
37}
38
39
40pj_complex PJFry::E0v1(const int i,
41 const double p1, const double p2, const double p3, const double p4, const double p5,
42 const double s12, const double s23, const double s34, const double s45, const double s15,
43 const double m5, const double m1, const double m2, const double m3, const double m4)
44{
45 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
46 return -ICache::getE(0, i, kinem);
47}
48
49pj_complex PJFry::E0v1(const int i,
50 const double p1, const double p2, const double p3, const double p4, const double p5,
51 const double s12, const double s23, const double s34, const double s45, const double s15,
52 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
53{
54 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
55 return -ICache::getE(abs(ep), i, kinem);
56}
57
58
59pj_complex PJFry::E0v2(const int i, const int j,
60 const double p1, const double p2, const double p3, const double p4, const double p5,
61 const double s12, const double s23, const double s34, const double s45, const double s15,
62 const double m5, const double m1, const double m2, const double m3, const double m4)
63{
64 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
65 return ICache::getE(0, i, j, kinem);
66}
67
68pj_complex PJFry::E0v2(const int i, const int j,
69 const double p1, const double p2, const double p3, const double p4, const double p5,
70 const double s12, const double s23, const double s34, const double s45, const double s15,
71 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
72{
73 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
74 return ICache::getE(abs(ep), i, j, kinem);
75}
76
77
78pj_complex PJFry::E0v3(const int i, const int j, const int k,
79 const double p1, const double p2, const double p3, const double p4, const double p5,
80 const double s12, const double s23, const double s34, const double s45, const double s15,
81 const double m5, const double m1, const double m2, const double m3, const double m4)
82{
83 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
84 return -ICache::getE(0, i, j, k, kinem);
85}
86
87pj_complex PJFry::E0v3(const int i, const int j, const int k,
88 const double p1, const double p2, const double p3, const double p4, const double p5,
89 const double s12, const double s23, const double s34, const double s45, const double s15,
90 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
91{
92 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
93 return -ICache::getE(abs(ep), i, j, k, kinem);
94}
95
96
97pj_complex PJFry::E0v4(const int i, const int j, const int k, const int l,
98 const double p1, const double p2, const double p3, const double p4, const double p5,
99 const double s12, const double s23, const double s34, const double s45, const double s15,
100 const double m5, const double m1, const double m2, const double m3, const double m4)
101{
102 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
103 return ICache::getE(0, i, j, k, l, kinem);
104}
105
106pj_complex PJFry::E0v4(const int i, const int j, const int k, const int l,
107 const double p1, const double p2, const double p3, const double p4, const double p5,
108 const double s12, const double s23, const double s34, const double s45, const double s15,
109 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
110{
111 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
112 return ICache::getE(abs(ep), i, j, k, l, kinem);
113}
114
115
116pj_complex PJFry::E0v5(const int i, const int j, const int k, const int l, const int m,
117 const double p1, const double p2, const double p3, const double p4, const double p5,
118 const double s12, const double s23, const double s34, const double s45, const double s15,
119 const double m5, const double m1, const double m2, const double m3, const double m4)
120{
121 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
122 return -ICache::getE(0, i, j, k, l, m, kinem);
123}
124
125pj_complex PJFry::E0v5(const int i, const int j, const int k, const int l, const int m,
126 const double p1, const double p2, const double p3, const double p4, const double p5,
127 const double s12, const double s23, const double s34, const double s45, const double s15,
128 const double m5, const double m1, const double m2, const double m3, const double m4, const int ep)
129{
130 Kinem5 kinem(p1, p2, p3, p4, p5, s12, s23, s34, s45, s15, m1, m2, m3, m4, m5);
131 return -ICache::getE(abs(ep), i, j, k, l, m, kinem);
132}
133
134/* ------------------------------------------------------------
135 * ------------------------------------------------------------
136 * BOX Integral section
137 * ------------------------------------------------------------
138 * ------------------------------------------------------------
139*/
140
141// No mass shift for scalar 4-poconst int D0
142
143pj_complex PJFry::D0v0(const double p1, const double p2, const double p3, const double p4,
144 const double s12, const double s23,
145 const double m1, const double m2, const double m3, const double m4)
146{
147 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
148 return ICache::getD(0, kinem);
149}
150
151pj_complex PJFry::D0v0(const double p1, const double p2, const double p3, const double p4,
152 const double s12, const double s23,
153 const double m1, const double m2, const double m3, const double m4, const int ep)
154{
155 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
156 return ICache::getD(abs(ep), kinem);
157}
158
159// Masses shifted to the left to comply with LoopTools notation
160
161pj_complex PJFry::D0v1(const int i,
162 const double p1, const double p2, const double p3, const double p4,
163 const double s12, const double s23,
164 const double m4, const double m1, const double m2, const double m3)
165{
166 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
167 return -ICache::getD(0, i, kinem);
168}
169
170pj_complex PJFry::D0v1(const int i,
171 const double p1, const double p2, const double p3, const double p4,
172 const double s12, const double s23,
173 const double m4, const double m1, const double m2, const double m3, const int ep)
174{
175 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
176 return -ICache::getD(abs(ep), i, kinem);
177}
178
179
180pj_complex PJFry::D0v2(const int i, const int j,
181 const double p1, const double p2, const double p3, const double p4,
182 const double s12, const double s23,
183 const double m4, const double m1, const double m2, const double m3)
184{
185 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
186 return ICache::getD(0, i, j, kinem);
187}
188
189pj_complex PJFry::D0v2(const int i, const int j,
190 const double p1, const double p2, const double p3, const double p4,
191 const double s12, const double s23,
192 const double m4, const double m1, const double m2, const double m3, const int ep)
193{
194 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
195 return ICache::getD(abs(ep), i, j, kinem);
196}
197
198
199pj_complex PJFry::D0v3(const int i, const int j, const int k,
200 const double p1, const double p2, const double p3, const double p4,
201 const double s12, const double s23,
202 const double m4, const double m1, const double m2, const double m3)
203{
204 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
205 return -ICache::getD(0, i, j, k, kinem);
206}
207
208pj_complex PJFry::D0v3(const int i, const int j, const int k,
209 const double p1, const double p2, const double p3, const double p4,
210 const double s12, const double s23,
211 const double m4, const double m1, const double m2, const double m3, const int ep)
212{
213 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
214 return -ICache::getD(abs(ep), i, j, k, kinem);
215}
216
217
218pj_complex PJFry::D0v4(const int i, const int j, const int k, const int l,
219 const double p1, const double p2, const double p3, const double p4,
220 const double s12, const double s23,
221 const double m4, const double m1, const double m2, const double m3)
222{
223 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
224 return ICache::getD(0, i, j, k, l, kinem);
225}
226
227pj_complex PJFry::D0v4(const int i, const int j, const int k, const int l,
228 const double p1, const double p2, const double p3, const double p4,
229 const double s12, const double s23,
230 const double m4, const double m1, const double m2, const double m3, const int ep)
231{
232 Kinem4 kinem(p1, p2, p3, p4, s12, s23, m1, m2, m3, m4);
233 return ICache::getD(abs(ep), i, j, k, l, kinem);
234}
235
236
237/* ------------------------------------------------------------
238 * ------------------------------------------------------------
239 * TRIANGLE Integral section
240 * ------------------------------------------------------------
241 * ------------------------------------------------------------
242*/
243
244// No mass shift for scalar 3-poconst int C0
245
246pj_complex PJFry::C0v0(const double p1, const double p2, const double p3,
247 const double m1, const double m2, const double m3)
248{
249 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
250 return ICache::getC(0, kinem);
251}
252
253pj_complex PJFry::C0v0(const double p1, const double p2, const double p3,
254 const double m1, const double m2, const double m3, const int ep)
255{
256 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
257 return ICache::getC(abs(ep), kinem);
258}
259
260// Masses shifted to the left to comply with LoopTools notation
261
262pj_complex PJFry::C0v1(const int i,
263 const double p1, const double p2, const double p3,
264 const double m3, const double m1, const double m2)
265{
266 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
267 return -ICache::getC(0, i, kinem);
268}
269
270pj_complex PJFry::C0v1(const int i,
271 const double p1, const double p2, const double p3,
272 const double m3, const double m1, const double m2, const int ep)
273{
274 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
275 return -ICache::getC(abs(ep), i, kinem);
276}
277
278
279pj_complex PJFry::C0v2(const int i, const int j,
280 const double p1, const double p2, const double p3,
281 const double m3, const double m1, const double m2)
282{
283 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
284 return ICache::getC(0, i, j, kinem);
285}
286
287pj_complex PJFry::C0v2(const int i, const int j,
288 const double p1, const double p2, const double p3,
289 const double m3, const double m1, const double m2, const int ep)
290{
291 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
292 return ICache::getC(abs(ep), i, j, kinem);
293}
294
295
296pj_complex PJFry::C0v3(const int i, const int j, const int k,
297 const double p1, const double p2, const double p3,
298 const double m3, const double m1, const double m2)
299{
300 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
301 return -ICache::getC(0, i, j, k, kinem);
302}
303
304pj_complex PJFry::C0v3(const int i, const int j, const int k,
305 const double p1, const double p2, const double p3,
306 const double m3, const double m1, const double m2, const int ep)
307{
308 Kinem3 kinem(p1, p2, p3, m1, m2, m3);
309 return -ICache::getC(abs(ep), i, j, k, kinem);
310}
311
312/* ------------------------------------------------------------
313 * ------------------------------------------------------------
314 * BUBBLE Integral section
315 * ------------------------------------------------------------
316 * ------------------------------------------------------------
317*/
318
319// Masses shifted to the left to comply with LoopTools notation
320
321pj_complex PJFry::B0v0(const double p1, const double m2, const double m1)
322{
323 Kinem2 kinem(p1, m1, m2);
324 return ICache::getB(0, kinem);
325}
326
327pj_complex PJFry::B0v0(const double p1, const double m2, const double m1, const int ep)
328{
329 Kinem2 kinem(p1, m1, m2);
330 return ICache::getB(abs(ep), kinem);
331}
332
333
334pj_complex PJFry::B0v1(const int i, const double p1, const double m2, const double m1)
335{
336 Kinem2 kinem(p1, m1, m2);
337 return -ICache::getB(0, i, kinem);
338}
339
340pj_complex PJFry::B0v1(const int i, const double p1, const double m2, const double m1, const int ep)
341{
342 Kinem2 kinem(p1, m1, m2);
343 return -ICache::getB(abs(ep), i, kinem);
344}
345
346
347pj_complex PJFry::B0v2(const int i, const int j, const double p1, const double m2, const double m1)
348{
349 Kinem2 kinem(p1, m1, m2);
350 return ICache::getB(0, i, j, kinem);
351}
352
353pj_complex PJFry::B0v2(const int i, const int j, const double p1, const double m2, const double m1, const int ep)
354{
355 Kinem2 kinem(p1, m1, m2);
356 return ICache::getB(abs(ep), i, j, kinem);
357}
358
359/* ------------------------------------------------------------
360* ------------------------------------------------------------
361* TADPOLE Integral section
362* ------------------------------------------------------------
363* ------------------------------------------------------------
364*/
365
366
367pj_complex PJFry::A0v0(const double m1)
368{
369 return ICache::getA(0, Kinem1(m1));
370}
371
372pj_complex PJFry::A0v0(const double m1, const int ep)
373{
374 return ICache::getA(abs(ep), Kinem1(m1));
375}
376
377/* ------------------------------------------------------------
378 * ------------------------------------------------------------
379 * MISC section
380 * ------------------------------------------------------------
381 * ------------------------------------------------------------
382 */
383
384
385double PJFry::GetMu2()
386{
387 return ICache::getMu2();
388}
389
390
391double PJFry::SetMu2(const double mu2)
392{
393 return ICache::setMu2(mu2);
394}
395
396void PJFry::ClearCache()
397{
398 // clear only coefficient cache and minor cache
399 // keep scalar integrals cache
400 ICache::ClearCC();
401 MCache::Clear();
402}
403
404// =============================================================
405// =============================================================
406// =============================================================
407
408double pgetmusq_()
409{
410 return PJFry::GetMu2();
411}
412
413void psetmusq_(double *mu2)
414{
415 PJFry::SetMu2(*mu2);
416}
417
418void pclearcache_()
419{
420 PJFry::ClearCache();
421}
422
423#ifdef USE_F2C
424// f2c,g77,ifort calling convention section (result is the first parameter)
425void pe0_(pj_complex *rslt,
426 double *p1, double *p2, double *p3, double *p4, double *p5,
427 double *s12, double *s23, double *s34, double *s45, double *s15,
428 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
429{
430 *rslt=PJFry::E0v0(*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
431}
432
433void pe0i_(pj_complex *rslt,
434 int *i,
435 double *p1, double *p2, double *p3, double *p4, double *p5,
436 double *s12, double *s23, double *s34, double *s45, double *s15,
437 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
438{
439 *rslt=PJFry::E0v1(*i,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
440}
441
442void pe0ij_(pj_complex *rslt,
443 int *i, int *j,
444 double *p1, double *p2, double *p3, double *p4, double *p5,
445 double *s12, double *s23, double *s34, double *s45, double *s15,
446 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
447{
448 *rslt=PJFry::E0v2(*i,*j,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
449}
450
451void pe0ijk_(pj_complex *rslt,
452 int *i, int *j, int *k,
453 double *p1, double *p2, double *p3, double *p4, double *p5,
454 double *s12, double *s23, double *s34, double *s45, double *s15,
455 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
456{
457 *rslt=PJFry::E0v3(*i,*j,*k,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
458}
459
460void pe0ijkl_(pj_complex *rslt,
461 int *i, int *j, int *k, int *l,
462 double *p1, double *p2, double *p3, double *p4, double *p5,
463 double *s12, double *s23, double *s34, double *s45, double *s15,
464 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
465{
466 *rslt=PJFry::E0v4(*i,*j,*k,*l,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
467}
468
469void pe0ijklm_(pj_complex *rslt,
470 int *i, int *j, int *k, int *l, int *m,
471 double *p1, double *p2, double *p3, double *p4, double *p5,
472 double *s12, double *s23, double *s34, double *s45, double *s15,
473 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
474{
475 *rslt=PJFry::E0v5(*i,*j,*k,*l,*m,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
476}
477
478void pd0_(pj_complex *rslt,
479 double *p1, double *p2, double *p3, double *p4,
480 double *s12, double *s23,
481 double *m1, double *m2, double *m3, double *m4, int *ep)
482{
483 *rslt=PJFry::D0v0(*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
484}
485
486void pd0i_(pj_complex *rslt,
487 int *i,
488 double *p1, double *p2, double *p3, double *p4,
489 double *s12, double *s23,
490 double *m1, double *m2, double *m3, double *m4, int *ep)
491{
492 *rslt=PJFry::D0v1(*i,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
493}
494
495void pd0ij_(pj_complex *rslt,
496 int *i, int *j,
497 double *p1, double *p2, double *p3, double *p4,
498 double *s12, double *s23,
499 double *m1, double *m2, double *m3, double *m4, int *ep)
500{
501 *rslt=PJFry::D0v2(*i,*j,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
502}
503
504void pd0ijk_(pj_complex *rslt,
505 int *i, int *j, int *k,
506 double *p1, double *p2, double *p3, double *p4,
507 double *s12, double *s23,
508 double *m1, double *m2, double *m3, double *m4, int *ep)
509{
510 *rslt=PJFry::D0v3(*i,*j,*k,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
511}
512
513void pd0ijkl_(pj_complex *rslt,
514 int *i, int *j, int *k, int *l,
515 double *p1, double *p2, double *p3, double *p4,
516 double *s12, double *s23,
517 double *m1, double *m2, double *m3, double *m4, int *ep)
518{
519 *rslt=PJFry::D0v4(*i,*j,*k,*l,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
520}
521
522void pc0_(pj_complex *rslt,
523 double *p1, double *p2, double *p3,
524 double *m1, double *m2, double *m3, int *ep)
525{
526 *rslt=PJFry::C0v0(*p1,*p2,*p3,*m1,*m2,*m3,*ep);
527}
528
529void pc0i_(pj_complex *rslt,
530 int *i,
531 double *p1, double *p2, double *p3,
532 double *m1, double *m2, double *m3, int *ep)
533{
534 *rslt=PJFry::C0v1(*i,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
535}
536
537void pc0ij_(pj_complex *rslt,
538 int *i, int *j,
539 double *p1, double *p2, double *p3,
540 double *m1, double *m2, double *m3, int *ep)
541{
542 *rslt=PJFry::C0v2(*i,*j,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
543}
544
545void pc0ijk_(pj_complex *rslt,
546 int *i, int *j, int *k,
547 double *p1, double *p2, double *p3,
548 double *m1, double *m2, double *m3, int *ep)
549{
550 *rslt=PJFry::C0v3(*i,*j,*k,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
551}
552
553void pb0_(pj_complex *rslt,
554 double *p1, double *m1, double *m2, int *ep)
555{
556 *rslt=PJFry::B0v0(*p1,*m1,*m2,*ep);
557}
558
559void pb0i_(pj_complex *rslt,
560 int *i,
561 double *p1, double *m1, double *m2, int *ep)
562{
563 *rslt=PJFry::B0v1(*i,*p1,*m1,*m2,*ep);
564}
565
566void pb0ij_(pj_complex *rslt,
567 int *i, int *j,
568 double *p1, double *m1, double *m2, int *ep)
569{
570 *rslt=PJFry::B0v2(*i,*j,*p1,*m1,*m2,*ep);
571}
572
573void pa0_(pj_complex *rslt, double *m1, int *ep)
574{
575 *rslt=PJFry::A0v0(*m1,*ep);
576}
577#else
578// GNU calling convention section (result is a complex return value)
579pj_complex pe0_( double *p1, double *p2, double *p3, double *p4, double *p5,
580 double *s12, double *s23, double *s34, double *s45, double *s15,
581 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
582{
583 return PJFry::E0v0(*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
584}
585
586pj_complex pe0i_(int *i,
587 double *p1, double *p2, double *p3, double *p4, double *p5,
588 double *s12, double *s23, double *s34, double *s45, double *s15,
589 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
590{
591 return PJFry::E0v1(*i,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
592}
593
594pj_complex pe0ij_(int *i, int *j,
595 double *p1, double *p2, double *p3, double *p4, double *p5,
596 double *s12, double *s23, double *s34, double *s45, double *s15,
597 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
598{
599 return PJFry::E0v2(*i,*j,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
600}
601
602pj_complex pe0ijk_(int *i, int *j, int *k,
603 double *p1, double *p2, double *p3, double *p4, double *p5,
604 double *s12, double *s23, double *s34, double *s45, double *s15,
605 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
606{
607 return PJFry::E0v3(*i,*j,*k,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
608}
609
610pj_complex pe0ijkl_(int *i, int *j, int *k, int *l,
611 double *p1, double *p2, double *p3, double *p4, double *p5,
612 double *s12, double *s23, double *s34, double *s45, double *s15,
613 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
614{
615 return PJFry::E0v4(*i,*j,*k,*l,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
616}
617
618pj_complex pe0ijklm_(int *i, int *j, int *k, int *l, int *m,
619 double *p1, double *p2, double *p3, double *p4, double *p5,
620 double *s12, double *s23, double *s34, double *s45, double *s15,
621 double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
622{
623 return PJFry::E0v5(*i,*j,*k,*l,*m,*p1,*p2,*p3,*p4,*p5,*s12,*s23,*s34,*s45,*s15,*m1,*m2,*m3,*m4,*m5,*ep);
624}
625
626pj_complex pd0_(double *p1, double *p2, double *p3, double *p4,
627 double *s12, double *s23,
628 double *m1, double *m2, double *m3, double *m4, int *ep)
629{
630 return PJFry::D0v0(*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
631}
632
633pj_complex pd0i_(int *i,
634 double *p1, double *p2, double *p3, double *p4,
635 double *s12, double *s23,
636 double *m1, double *m2, double *m3, double *m4, int *ep)
637{
638 return PJFry::D0v1(*i,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
639}
640
641pj_complex pd0ij_(int *i, int *j,
642 double *p1, double *p2, double *p3, double *p4,
643 double *s12, double *s23,
644 double *m1, double *m2, double *m3, double *m4, int *ep)
645{
646 return PJFry::D0v2(*i,*j,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
647}
648
649pj_complex pd0ijk_(int *i, int *j, int *k,
650 double *p1, double *p2, double *p3, double *p4,
651 double *s12, double *s23,
652 double *m1, double *m2, double *m3, double *m4, int *ep)
653{
654 return PJFry::D0v3(*i,*j,*k,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
655}
656
657pj_complex pd0ijkl_(int *i, int *j, int *k, int *l,
658 double *p1, double *p2, double *p3, double *p4,
659 double *s12, double *s23,
660 double *m1, double *m2, double *m3, double *m4, int *ep)
661{
662 return PJFry::D0v4(*i,*j,*k,*l,*p1,*p2,*p3,*p4,*s12,*s23,*m1,*m2,*m3,*m4,*ep);
663}
664
665pj_complex pc0_(double *p1, double *p2, double *p3,
666 double *m1, double *m2, double *m3, int *ep)
667{
668 return PJFry::C0v0(*p1,*p2,*p3,*m1,*m2,*m3,*ep);
669}
670
671pj_complex pc0i_(int *i,
672 double *p1, double *p2, double *p3,
673 double *m1, double *m2, double *m3, int *ep)
674{
675 return PJFry::C0v1(*i,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
676}
677
678pj_complex pc0ij_(int *i, int *j,
679 double *p1, double *p2, double *p3,
680 double *m1, double *m2, double *m3, int *ep)
681{
682 return PJFry::C0v2(*i,*j,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
683}
684
685pj_complex pc0ijk_(int *i, int *j, int *k,
686 double *p1, double *p2, double *p3,
687 double *m1, double *m2, double *m3, int *ep)
688{
689 return PJFry::C0v3(*i,*j,*k,*p1,*p2,*p3,*m1,*m2,*m3,*ep);
690}
691
692pj_complex pb0_( double *p1, double *m1, double *m2, int *ep)
693{
694 return PJFry::B0v0(*p1,*m1,*m2,*ep);
695}
696
697pj_complex pb0i_(int *i,
698 double *p1, double *m1, double *m2, int *ep)
699{
700 return PJFry::B0v1(*i,*p1,*m1,*m2,*ep);
701}
702
703pj_complex pb0ij_(int *i, int *j,
704 double *p1, double *m1, double *m2, int *ep)
705{
706 return PJFry::B0v2(*i,*j,*p1,*m1,*m2,*ep);
707}
708
709pj_complex pa0_( double *m1, int *ep)
710{
711 return PJFry::A0v0(*m1,*ep);
712}
713#endif
ICache::setMu2
static double setMu2(const double newmu2)
Definition: eemmg-lib-new/src/cache.cpp:25
ICache::getMu2
static double getMu2()
Definition: eemmg-lib-new/src/cache.cpp:20
ICache::getD
static ncomplex getD(int ep, const Kinem4 &kin)
Definition: eemmg-lib-new/src/cache.h:119
ICache::getA
static ncomplex getA(int ep, const Kinem1 &kin)
Definition: eemmg-lib-new/src/cache.h:137
ICache::getB
static ncomplex getB(int ep, const Kinem2 &kin)
Definition: eemmg-lib-new/src/cache.h:132
ICache::getC
static ncomplex getC(int ep, const Kinem3 &kin)
Definition: eemmg-lib-new/src/cache.h:126
ICache::getE
static ncomplex getE(int ep, const Kinem5 &kin)
Definition: eemmg-lib-new/src/cache.cpp:233
ICache::ClearCC
static void ClearCC()
Definition: eemmg-lib-new/src/cache.cpp:47
MCache::Clear
static void Clear()
Definition: eemmg-lib-new/src/cache.cpp:69
pc0_
pj_complex pc0_(double *p1, double *p2, double *p3, double *m1, double *m2, double *m3, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:665
pb0ij_
pj_complex pb0ij_(int *i, int *j, double *p1, double *m1, double *m2, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:703
pd0i_
pj_complex pd0i_(int *i, double *p1, double *p2, double *p3, double *p4, double *s12, double *s23, double *m1, double *m2, double *m3, double *m4, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:633
pe0ijklm_
pj_complex pe0ijklm_(int *i, int *j, int *k, int *l, int *m, double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:618
pc0i_
pj_complex pc0i_(int *i, double *p1, double *p2, double *p3, double *m1, double *m2, double *m3, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:671
pgetmusq_
double pgetmusq_()
Definition: eemmg-lib-new/src/pjfry.cpp:408
pb0_
pj_complex pb0_(double *p1, double *m1, double *m2, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:692
pc0ij_
pj_complex pc0ij_(int *i, int *j, double *p1, double *p2, double *p3, double *m1, double *m2, double *m3, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:678
pd0ijkl_
pj_complex pd0ijkl_(int *i, int *j, int *k, int *l, double *p1, double *p2, double *p3, double *p4, double *s12, double *s23, double *m1, double *m2, double *m3, double *m4, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:657
pclearcache_
void pclearcache_()
Definition: eemmg-lib-new/src/pjfry.cpp:418
pc0ijk_
pj_complex pc0ijk_(int *i, int *j, int *k, double *p1, double *p2, double *p3, double *m1, double *m2, double *m3, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:685
pd0ijk_
pj_complex pd0ijk_(int *i, int *j, int *k, double *p1, double *p2, double *p3, double *p4, double *s12, double *s23, double *m1, double *m2, double *m3, double *m4, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:649
pa0_
pj_complex pa0_(double *m1, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:709
pe0ijkl_
pj_complex pe0ijkl_(int *i, int *j, int *k, int *l, double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:610
psetmusq_
void psetmusq_(double *mu2)
Definition: eemmg-lib-new/src/pjfry.cpp:413
pe0i_
pj_complex pe0i_(int *i, double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:586
pe0ij_
pj_complex pe0ij_(int *i, int *j, double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:594
pe0_
pj_complex pe0_(double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:579
pd0_
pj_complex pd0_(double *p1, double *p2, double *p3, double *p4, double *s12, double *s23, double *m1, double *m2, double *m3, double *m4, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:626
pe0ijk_
pj_complex pe0ijk_(int *i, int *j, int *k, double *p1, double *p2, double *p3, double *p4, double *p5, double *s12, double *s23, double *s34, double *s45, double *s15, double *m1, double *m2, double *m3, double *m4, double *m5, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:602
pb0i_
pj_complex pb0i_(int *i, double *p1, double *m1, double *m2, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:697
pd0ij_
pj_complex pd0ij_(int *i, int *j, double *p1, double *p2, double *p3, double *p4, double *s12, double *s23, double *m1, double *m2, double *m3, double *m4, int *ep)
Definition: eemmg-lib-new/src/pjfry.cpp:641
m3
double double double double double * m3
Definition: qcdloop1.h:76
p4
double double double * p4
Definition: qcdloop1.h:77
s12
double double double double * s12
Definition: qcdloop1.h:77
s23
double double double double double * s23
Definition: qcdloop1.h:77
m4
double double double double double double double double double * m4
Definition: qcdloop1.h:77
p2
double * p2
Definition: qcdloop1.h:76
ep
double int * ep
Definition: qcdloop1.h:74
mu2
double * mu2
Definition: qcdloop1.h:74
m1
double * m1
Definition: qcdloop1.h:75
p3
double double * p3
Definition: qcdloop1.h:76
m2
double double * m2
Definition: qcdloop1.h:75