2 parameter (e1=
"ffca0: minimum value complex logarit"
3 +
"hm gives problem, change mu.")
5 parameter (e2="
ffxa0: minimum value real logarithm "
6 + "gives problem, change mu.")
8 parameter (e3="
ffcb0: minimum value complex logarit"
9 + "hm gives problem, change mu.")
11 parameter (e4="
ffxb0: minimum value real logarithm "
12 + "gives problem, change mu.")
16 parameter (e6="
ffcb0p: minimum value complex
log giv"
17 + "es problem in unequal masses.")
19 parameter (e7="
ffxb0p: divergence for k->0, m1=m2=0.")
21 parameter (e8="
ffxb0p: minimum value real
log gives "
22 + "problem in equal masses.")
24 parameter (e9="
ffxb0p: minimum value real
log gives "
25 + "problem in unequal masses.")
28 + "omenta and one zero.")
31 + "omenta and one zero.")
33 parameter (e12="
ffcs3: illegal code for isoort(1) (s"
36 parameter (e13="
ffcs3: illegal code for isoort(2) (s"
39 parameter (e14="
ffcs3: imaginary part wrong, will be"
42 parameter (e15="
ffcs3: isoort = -1,0
not yet ready")
44 parameter (e16="
ffcs3: illegal combination in isoort"
45 + " (should
not occur)")
47 parameter (e17="
ffcxs3: illegal code for isoort(1) (s"
50 parameter (e18="
ffcxs3: illegal code for isoort(2) (s"
53 parameter (e19="
ffcs4: imaginary part is wrong (shou"
56 parameter (e20="
ffdcrr: Taylor expansion in 1/x
not y"
59 parameter (e21="
ffdcxr: imaginary part is wrong")
61 parameter (e22="
ffdcxr: Taylor expansion in 1/x
not y"
64 parameter (e23="
ffcrr: minimum value complex
log cau"
65 + "ses correction term
to be wrong.")
67 parameter (e24="
ffcxr: minimum value real
log causes"
68 + " correction term
to be wrong.")
70 parameter (e25="
ffcrr: illegal code for iclas1 (shou"
73 parameter (e26="
ffcxr: illegal code for iclas1 (shou"
76 parameter (e27="
ffcrr: illegal code for iclas2 (shou"
79 parameter (e28="
ffcxr: illegal code for iclas2 (shou"
82 parameter (e29="
ffxli2: argument too large (should no"
85 parameter (e30="
ffzli2: argument too large (should no"
88 parameter (e31="
ffzzdl: imaginary part dilog is undef"
89 + "ined for real x > 1.")
92 + "egative numbers a,b, ab.")
95 + "egative numbers a,b, ab.")
97 parameter (e34="ffcota: illegal flag (should
not occu"
100 parameter (e35="
ffrota: illegal flag (should
not occu"
103 parameter (e36="
ffccyz: I took the wrong value for ca"
104 + "lpha... (should
not occur)")
106 parameter (e37="
ffxxyz: I took the wrong value for
al"
107 + "pha... (should
not occur)")
109 parameter (e38="
ffcoot: a=0, trying
to find two roots"
110 + " of a linear equation ...")
112 parameter (e39="
ffroot: a=0, trying
to find two roots"
113 + " of a linear equation ...")
115 parameter (e40="
ffrot3: all three external masses zer"
118 parameter (e41="
ffxc0: lambda(p1,p2,p3) < 0, unphysi"
119 + "cal configuration")
121 parameter (e42="
ffxc0: cannot
handle this case (p1,p"
122 + "2,p3 dependent, on threshold)")
124 parameter (e43="
ffcxs3: illegal code for isoort(1) (s"
125 + "hould
not occur)")
127 parameter (e44="
ffxd0: lambda(p1,p2,p3,p4) < 0, unph"
128 + "ysical configuration")
130 parameter (e45="
ffxd0: cannot
handle this case (p1,p"
131 + "2,p3 dependent, on threshold)")
133 parameter (e46="
ffxd0p: correction terms for Ai <0 in"
134 + "finite (mass zero?)")
136 parameter (e47="
ffcxyz: p_i^2 = 0 (should
not occur)")
138 parameter (e48="ffeta: answer
not consistent with no"
139 + "rmal result (old)")
141 parameter (e49="
ffcc0: cannot
handle complex externa"
142 + "l momenta or im > 0")
144 parameter (e50="ffcd0: cannot
handle complex externa"
147 parameter (e51="
zfflog: imaginary part undefined for "
150 parameter (e52="
zxfflg: imaginary part undefined for "
153 parameter (e53="
ffcs3:
eta changes within (0,1), add"
154 + " sophisticated terms...")
156 parameter (e54="
ffrot4: cannot find any physical vert"
157 + "ex
to apply transformation.")
159 parameter (e55="fftra0: too many vectors parallel, p_"
160 + "1.p_7 or p_2.p_7 is zero.")
162 parameter (e56="
zfflog: tiny imaginary part in confli"
163 + "ct with ieps prescription.")
165 parameter (e57="
ffxe0: lambda(p1,p2,p3,p4,p5) < 0, u"
168 parameter (e58="
ffxc0j: IR divergent
C0 with lambda(p"
171 parameter (e59="
ffxc0i: IR divergent
C0 with lambda2=0.")
173 parameter (e60="
ffxc0j: IR divergent
C0 obtained from"
174 + "
D0 is singular. Contact author.")
176 parameter (e61="
ffxd0p: IR divergent
D0 with lambda2=0.")
178 parameter (e62="
ffxc0p: I never expected complex root"
179 + "s in an IR divergent diagram.")
182 + "ence per 3point function")
185 + " (3,4), rotated wrongly.")
187 parameter (e65="
ffcxr: IR divergence but iclass!=3. "
188 + " should
not occur.")
190 parameter (e66="
ffcxs3: different imaginary signs sho"
191 + "uld
not occur for ipole=3.")
193 parameter (e67="
ffxdbd: I cannot use this algorithm
f"
194 + "or a linear IR divergence")
196 parameter (e68="
ffxd0: cannot find a proj. transform"
197 + "ation; try another permutation.")
199 parameter (e69="ff5ind: could
not find independent mo"
200 + "menta (should
not occur).")
202 parameter (e70="ffxdna: lambda(
pi,pj,pk) < 0, unphysi"
203 + "cal configuration")
205 parameter (e71="ffxdna: cannot
handle lambda(
pi,pj,pk"
206 + ") = 0, dependent momenta.")
208 parameter (e72="
ffxd0e: could
not find a stable root;"
209 + " please try another permutation")
212 + "rgent four point function (yet)")
214 parameter (e74="
ffxdbd: IR divergent
B0' without cuto"
217 parameter (e75="
ffdcxr: dyz=0, should
not occur")
219 parameter (e76="
ffdcrr: cdwz=0, but iepsz!=iepsz and "
222 parameter (e77="
ffdcrr: cdyz=0, should
not occur")
224 parameter (e78="
ffdcc0: imaginary part wrong")
228 parameter (e80="
ffdcs: mixed
up iep's, 2*
pi^2 wrong "
231 parameter (e81="
ffdcs: wrong value for isoort")
233 parameter (e82="
ffdxc0: imaginary part Ai < 0 terms unc"
236 parameter (e83="
ffxc0j: sorry, complex roots
not yet "
239 parameter (e84="
ffxc0p: imaginary part Ai < 0 terms unc"
242 parameter (e85="
ffxd0a: t3 = t4, don''t know what
to do")
244 parameter (e86="
ffxdbp: cannot compute derivative, la"
247 parameter (e87="ffxdi: dependent momenta
not yet sup"
248 + "ported (boundary of
phase space)")
250 parameter (e88="
ffxxyz: xk = 0
not yet implemented")
252 parameter (e92="ffxc1: cannot invert matrix with zer"
255 parameter (e93="ffze0:
Im(m^2) > 0")
257 parameter (e94="ffze0:
Im(p^2) != 0")
259 parameter (e95="ffzf0:
Im(m^2) > 0")
261 parameter (e96="ffzf0:
Im(p^2) != 0")
263 parameter (e97="
ffxc0j: ill-defined IR-divergent
C0 "
264 + "for massless charged particles.")
266 parameter (e98="
ffxdbd: ill-defined IR-divergent
D0 "
267 + "for massless charged particles.")
269 parameter (e100="ffrcvr: probably underflow, I do"
270 + "
not know where or how severe.")
272 parameter (e101="
ffxdb1: case
not defined")
278 parameter (e104="
ffwbeta: prefactor 1/(SV-TU) = 1/0 "
281 parameter (e105="
ffT_lin: prefactor 1/(SV-TU) = 1/0 "
284 parameter (e99="
ffT13: prefactor 1/(SV-TU) = 1/0 "
287 parameter (e89="
ffS2:
log(0) singularity")
289 parameter (e90="
ffS3n:
end-point singularity")
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 /
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.
std::complex< T > phase(const std::complex< T > &z) noexcept
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)
std::complex< double > f(double tau) noexcept
double * end(GSL_vector &v)
iterator to end of GSL_vector
Complex< T > log(const Complex< T > &z) noexcept