flexiblesusy is hosted by Hepforge, IPPP Durham
FlexibleSUSY
fferr.h
Go to the documentation of this file.
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 /
#define ffrota
Definition: externals.h:191
#define ffxdbd
Definition: externals.h:239
#define zfflog
Definition: externals.h:290
#define ffcxr
Definition: externals.h:143
#define ffcxs3
Definition: externals.h:144
#define ffcb0
Definition: externals.h:107
#define ffxdb11
Definition: externals.h:237
#define ffT_lin
Definition: externals.h:204
#define ffdcxr
Definition: externals.h:152
#define ffcrr
Definition: externals.h:138
#define ffxa0
Definition: externals.h:211
#define eta
Definition: externals.h:94
#define ffxli2
Definition: externals.h:250
#define ffxc0p
Definition: externals.h:225
#define ffcs3
Definition: externals.h:140
#define ffd0c
Definition: externals.h:147
#define zxfflg
Definition: externals.h:291
#define nffet1
Definition: externals.h:277
#define ffT13
Definition: externals.h:200
#define ffxdbp
Definition: externals.h:240
#define ffccyz
Definition: externals.h:118
#define ffzli2
Definition: externals.h:258
#define ffxd0e
Definition: externals.h:231
#define ffxd0a
Definition: externals.h:228
#define ffrot4
Definition: externals.h:190
#define ffwbeta
Definition: externals.h:209
#define ffxb0p
Definition: externals.h:213
#define ffcxyz
Definition: externals.h:146
#define ffcc0p
Definition: externals.h:116
#define ffdcrr
Definition: externals.h:150
#define ffxdb1
Definition: externals.h:238
#define ffS3n
Definition: externals.h:195
#define ffxdir
Definition: externals.h:241
#define ffxc0i
Definition: externals.h:221
#define ffdcc0
Definition: externals.h:149
#define ffcoot
Definition: externals.h:135
#define ffxd0p
Definition: externals.h:234
#define ffzzdl
Definition: externals.h:260
#define ffrot3
Definition: externals.h:189
#define ffxc0
Definition: externals.h:223
#define ffdcs
Definition: externals.h:151
#define ffcs4
Definition: externals.h:141
#define ffxd0
Definition: externals.h:232
#define ffxe0
Definition: externals.h:244
#define ffxc0j
Definition: externals.h:222
#define ffcb0p
Definition: externals.h:108
#define ffS2
Definition: externals.h:194
#define nffeta
Definition: externals.h:278
#define ffcc0
Definition: externals.h:115
#define ffxxyz
Definition: externals.h:254
#define ffroot
Definition: externals.h:187
#define ffxb0
Definition: externals.h:212
#define ffdxc0
Definition: externals.h:174
without even the implied warranty of !MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU !General Public License for more details !You should have received a copy of the GNU General Public License !along with FlexibleSUSY If not
Target to(const Source &arg)
Converts an object of type Source to an object of type Target.
Definition: slhaea.h:50
std::complex< T > phase(const std::complex< T > &z) noexcept
Definition: ckm.cpp:236
double C0(double m1, double m2, double m3) noexcept
(arguments are interpreted as unsquared)
double D0(double m1, double m2, double m3, double m4) noexcept
(arguments are interpreted as unsquared)
double B0(double m1, double m2, double scale) noexcept
(arguments are interpreted as unsquared)
const double al
Definition: consts.hpp:28
std::complex< double > f(double tau) noexcept
double * end(GSL_vector &v)
iterator to end of GSL_vector
Definition: gsl_vector.cpp:254
constexpr double pi
Definition: mathdefs.hpp:16
Complex< T > log(const Complex< T > &z) noexcept
Definition: complex.hpp:54
generate run handle
data
Definition: scan_HSSUSY.m:46
#define Im
Definition: types.h:13