herwig is hosted by Hepforge, IPPP Durham
Herwig  7.2.1
fferr.h
1  character*80 e1
2  parameter (e1="ffca0: minimum value complex logarit"//
3  + "hm gives problem, change mu.")
4  character*80 e2
5  parameter (e2="ffxa0: minimum value real logarithm "//
6  + "gives problem, change mu.")
7  character*80 e3
8  parameter (e3="ffcb0: minimum value complex logarit"//
9  + "hm gives problem, change mu.")
10  character*80 e4
11  parameter (e4="ffxb0: minimum value real logarithm "//
12  + "gives problem, change mu.")
13  character*80 e5
14  parameter (e5="ffcb0p: cannot handle complex k^2 yet")
15  character*80 e6
16  parameter (e6="ffcb0p: minimum value complex log giv"//
17  + "es problem in unequal masses.")
18  character*80 e7
19  parameter (e7="ffxb0p: divergence for k->0, m1=m2=0.")
20  character*80 e8
21  parameter (e8="ffxb0p: minimum value real log gives "//
22  + "problem in equal masses.")
23  character*80 e9
24  parameter (e9="ffxb0p: minimum value real log gives "//
25  + "problem in unequal masses.")
26  character*80 e10
27  parameter (e10="ffcc0p: cannot handle two spacelike m"//
28  + "omenta and one zero.")
29  character*80 e11
30  parameter (e11="ffxc0p: cannot handle two spacelike m"//
31  + "omenta and one zero.")
32  character*80 e12
33  parameter (e12="ffcs3: illegal code for isoort(1) (s"//
34  + "hould not occur)")
35  character*80 e13
36  parameter (e13="ffcs3: illegal code for isoort(2) (s"//
37  + "hould not occur)")
38  character*80 e14
39  parameter (e14="ffcs3: imaginary part wrong, will be"//
40  + " improved later")
41  character*80 e15
42  parameter (e15="ffcs3: isoort = -1,0 not yet ready")
43  character*80 e16
44  parameter (e16="ffcs3: illegal combination in isoort"//
45  + " (should not occur)")
46  character*80 e17
47  parameter (e17="ffcxs3: illegal code for isoort(1) (s"//
48  + "hould not occur)")
49  character*80 e18
50  parameter (e18="ffcxs3: illegal code for isoort(2) (s"//
51  + "hould not occur)")
52  character*80 e19
53  parameter (e19="ffcs4: imaginary part is wrong (shou"//
54  + "ld be updated)")
55  character*80 e20
56  parameter (e20="ffdcrr: Taylor expansion in 1/x not y"//
57  + "et ready")
58  character*80 e21
59  parameter (e21="ffdcxr: imaginary part is wrong")
60  character*80 e22
61  parameter (e22="ffdcxr: Taylor expansion in 1/x not y"//
62  + "et ready")
63  character*80 e23
64  parameter (e23="ffcrr: minimum value complex log cau"//
65  + "ses correction term to be wrong.")
66  character*80 e24
67  parameter (e24="ffcxr: minimum value real log causes"//
68  + " correction term to be wrong.")
69  character*80 e25
70  parameter (e25="ffcrr: illegal code for iclas1 (shou"//
71  + "ld not occur)")
72  character*80 e26
73  parameter (e26="ffcxr: illegal code for iclas1 (shou"//
74  + "ld not occur)")
75  character*80 e27
76  parameter (e27="ffcrr: illegal code for iclas2 (shou"//
77  + "ld not occur)")
78  character*80 e28
79  parameter (e28="ffcxr: illegal code for iclas2 (shou"//
80  + "ld not occur)")
81  character*80 e29
82  parameter (e29="ffxli2: argument too large (should no"//
83  + "t occur)")
84  character*80 e30
85  parameter (e30="ffzli2: argument too large (should no"//
86  + "t occur)")
87  character*80 e31
88  parameter (e31="ffzzdl: imaginary part dilog is undef"//
89  + "ined for real x > 1.")
90  character*80 e32
91  parameter (e32="nffeta: eta is not defined for real n"//
92  + "egative numbers a,b, ab.")
93  character*80 e33
94  parameter (e33="nffet1: eta is not defined for real n"//
95  + "egative numbers a,b, ab.")
96  character*80 e34
97  parameter (e34="ffcota: illegal flag (should not occu"//
98  + "r)")
99  character*80 e35
100  parameter (e35="ffrota: illegal flag (should not occu"//
101  + "r)")
102  character*80 e36
103  parameter (e36="ffccyz: I took the wrong value for ca"//
104  + "lpha... (should not occur)")
105  character*80 e37
106  parameter (e37="ffxxyz: I took the wrong value for al"//
107  + "pha... (should not occur)")
108  character*80 e38
109  parameter (e38="ffcoot: a=0, trying to find two roots"//
110  + " of a linear equation ...")
111  character*80 e39
112  parameter (e39="ffroot: a=0, trying to find two roots"//
113  + " of a linear equation ...")
114  character*80 e40
115  parameter (e40="ffrot3: all three external masses zer"//
116  + "o !")
117  character*80 e41
118  parameter (e41="ffxc0: lambda(p1,p2,p3) < 0, unphysi"//
119  + "cal configuration")
120  character*80 e42
121  parameter (e42="ffxc0: cannot handle this case (p1,p"//
122  + "2,p3 dependent, on threshold)")
123  character*80 e43
124  parameter (e43="ffcxs3: illegal code for isoort(1) (s"//
125  + "hould not occur)")
126  character*80 e44
127  parameter (e44="ffxd0: lambda(p1,p2,p3,p4) < 0, unph"//
128  + "ysical configuration")
129  character*80 e45
130  parameter (e45="ffxd0: cannot handle this case (p1,p"//
131  + "2,p3 dependent, on threshold)")
132  character*80 e46
133  parameter (e46="ffxd0p: correction terms for Ai <0 in"//
134  + "finite (mass zero?)")
135  character*80 e47
136  parameter (e47="ffcxyz: p_i^2 = 0 (should not occur)")
137  character*80 e48
138  parameter (e48="ffeta: answer not consistent with no"//
139  + "rmal result (old)")
140  character*80 e49
141  parameter (e49="ffcc0: cannot handle complex externa"//
142  + "l momenta or im > 0")
143  character*80 e50
144  parameter (e50="ffcd0: cannot handle complex externa"//
145  + "l momenta.")
146  character*80 e51
147  parameter (e51="zfflog: imaginary part undefined for "//
148  + "real z < 0.")
149  character*80 e52
150  parameter (e52="zxfflg: imaginary part undefined for "//
151  + "x < 0.")
152  character*80 e53
153  parameter (e53="ffcs3: eta changes within (0,1), add"//
154  + " sophisticated terms...")
155  character*80 e54
156  parameter (e54="ffrot4: cannot find any physical vert"//
157  + "ex to apply transformation.")
158  character*80 e55
159  parameter (e55="fftra0: too many vectors parallel, p_"//
160  + "1.p_7 or p_2.p_7 is zero.")
161  character*80 e56
162  parameter (e56="zfflog: tiny imaginary part in confli"//
163  + "ct with ieps prescription.")
164  character*80 e57
165  parameter (e57="ffxe0: lambda(p1,p2,p3,p4,p5) < 0, u"//
166  + "nphysical")
167  character*80 e58
168  parameter (e58="ffxc0j: IR divergent C0 with lambda(p"//
169  + "1,p2,p3)=0.")
170  character*80 e59
171  parameter (e59="ffxc0i: IR divergent C0 with lambda2=0.")
172  character*80 e60
173  parameter (e60="ffxc0j: IR divergent C0 obtained from"//
174  + " D0 is singular. Contact author.")
175  character*80 e61
176  parameter (e61="ffxd0p: IR divergent D0 with lambda2=0.")
177  character*80 e62
178  parameter (e62="ffxc0p: I never expected complex root"//
179  + "s in an IR divergent diagram.")
180  character*80 e63
181  parameter (e63="ffxd0p: can only handle one IR diverg"//
182  + "ence per 3point function")
183  character*80 e64
184  parameter (e64="ffxd0p: cannot handle a threshold in"//
185  + " (3,4), rotated wrongly.")
186  character*80 e65
187  parameter (e65="ffcxr: IR divergence but iclass!=3. "//
188  + " should not occur.")
189  character*80 e66
190  parameter (e66="ffcxs3: different imaginary signs sho"//
191  + "uld not occur for ipole=3.")
192  character*80 e67
193  parameter (e67="ffxdbd: I cannot use this algorithm f"//
194  + "or a linear IR divergence")
195  character*80 e68
196  parameter (e68="ffxd0: cannot find a proj. transform"//
197  + "ation; try another permutation.")
198  character*80 e69
199  parameter (e69="ff5ind: could not find independent mo"//
200  + "menta (should not occur).")
201  character*80 e70
202  parameter (e70="ffxdna: lambda(pi,pj,pk) < 0, unphysi"//
203  + "cal configuration")
204  character*80 e71
205  parameter (e71="ffxdna: cannot handle lambda(pi,pj,pk"//
206  + ") = 0, dependent momenta.")
207  character*80 e72
208  parameter (e72="ffxd0e: could not find a stable root;"//
209  + " please try another permutation")
210  character*80 e73
211  parameter (e73="ffxdir: cannot handle a linearly dive"//
212  + "rgent four point function (yet)")
213  character*80 e74
214  parameter (e74="ffxdbd: IR divergent B0' without cuto"//
215  + "ff in /ffregul/")
216  character*80 e75
217  parameter (e75="ffdcxr: dyz=0, should not occur")
218  character*80 e76
219  parameter (e76="ffdcrr: cdwz=0, but iepsz!=iepsz and "//
220  + "significant")
221  character*80 e77
222  parameter (e77="ffdcrr: cdyz=0, should not occur")
223  character*80 e78
224  parameter (e78="ffdcc0: imaginary part wrong")
225  character*80 e79
226  parameter (e79="ffdcs: cannot handle isoort=0")
227  character*80 e80
228  parameter (e80="ffdcs: mixed up iep's, 2*pi^2 wrong "//
229  + "somewhere")
230  character*80 e81
231  parameter (e81="ffdcs: wrong value for isoort")
232  character*80 e82
233  parameter (e82="ffdxc0: imaginary part Ai < 0 terms unc"//
234  + "ertain")
235  character*80 e83
236  parameter (e83="ffxc0j: sorry, complex roots not yet "//
237  + "supported here")
238  character*80 e84
239  parameter (e84="ffxc0p: imaginary part Ai < 0 terms unc"//
240  + "ertain")
241  character*80 e85
242  parameter (e85="ffxd0a: t3 = t4, don''t know what to do")
243  character*80 e86
244  parameter (e86="ffxdbp: cannot compute derivative, la"//
245  + "m=0")
246  character*80 e87
247  parameter (e87="ffxdi: dependent momenta not yet sup"//
248  + "ported (boundary of phase space)")
249  character*80 e88
250  parameter (e88="ffxxyz: xk = 0 not yet implemented")
251  character*80 e92
252  parameter (e92="ffxc1: cannot invert matrix with zer"//
253  + "o determinant.")
254  character*80 e93
255  parameter (e93="ffze0: Im(m^2) > 0")
256  character*80 e94
257  parameter (e94="ffze0: Im(p^2) != 0")
258  character*80 e95
259  parameter (e95="ffzf0: Im(m^2) > 0")
260  character*80 e96
261  parameter (e96="ffzf0: Im(p^2) != 0")
262  character*80 e97
263  parameter (e97="ffxc0j: ill-defined IR-divergent C0 "//
264  + "for massless charged particles.")
265  character*80 e98
266  parameter (e98="ffxdbd: ill-defined IR-divergent D0 "//
267  + "for massless charged particles.")
268  character*80 e100
269  parameter (e100="ffrcvr: probably underflow, I do"//
270  + " not know where or how severe.")
271  character*80 e101
272  parameter (e101="ffxdb1: case not defined")
273  character*80 e102
274  parameter (e102="ffxdb11: case not defined")
275  character*80 e103
276  parameter (e103="ffd0c: cannot handle this case")
277  character*80 e104
278  parameter (e104="ffwbeta: prefactor 1/(SV-TU) = 1/0 "//
279  + "for all y")
280  character*80 e105
281  parameter (e105="ffT_lin: prefactor 1/(SV-TU) = 1/0 "//
282  + "for all y")
283  character*80 e99
284  parameter (e99="ffT13: prefactor 1/(SV-TU) = 1/0 "//
285  + "for all y")
286  character*80 e89
287  parameter (e89="ffS2: log(0) singularity")
288  character*80 e90
289  parameter (e90="ffS3n: end-point singularity")
290  character*80 e91
291  parameter (e91="ffS3n: log(0) singularity")
292  character*80 error(105)
293  data error / e1,e2,e3,e4,e5,e6,e7,e8,e9,
294  + e10,e11,e12,e13,e14,e15,e16,e17,e18,e19,
295  + e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,
296  + e30,e31,e32,e33,e34,e35,e36,e37,e38,e39,
297  + e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,
298  + e50,e51,e52,e53,e54,e55,e56,e57,e58,e59,
299  + e60,e61,e62,e63,e64,e65,e66,e67,e68,e69,
300  + e70,e71,e72,e73,e74,e75,e76,e77,e78,e79,
301  + e80,e81,e82,e83,e84,e85,e86,e87,e88,e89,
302  + e90,e91,e92,e93,e94,e95,e96,e97,e98,e99,
303  + e100,e101,e102,e103,e104,e105 /
Generator< Zero > zero(double low, double up)
Construct a zero density.