Commit cacf0fa4 authored by ulrich_y's avatar ulrich_y
Browse files

Bug-fix branch cut if specified

In (60) of VW the branch cut was wrong (also in
Mathematica). I misunderstood the meaning of

if (a.back() > 0) {
    new_pending_integrals.push_back(-scale);
    result += I*Pi;
} else {
    new_pending_integrals.push_back(scale);
    result -= I*Pi;
}

in GiNaC
parent 7840656d
......@@ -127,7 +127,7 @@ CONTAINS
! + G_flat(p(2:size(p)), p(1)) * log(-sub_ieps(g(1)))
!TODO
res = pending_integral(p,2,[g(1)]) - pending_integral(p,2,[izero]) &
+ G_flat(p(2:size(p)), p(1)) * conjg(log(-g(1)%c))
+ G_flat(p(2:size(p)), p(1)) * (log(g(1)%c) + p(1)%i0 * pi * (0.,1.))
return
end if
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment