<<FIESTA_1.0.0.m; (* everything performed on a dual-core laptop; most of the time on the box is due to opening and closing the databases - hundreds of files on disk *)
FIESTA, version 1.0.0
Box:
PrimarySectorCoefficients={2,2,0,0};
SDEvaluate[UF[{k},{-k^2,-(k+p1)^2,-(k+p1+p2)^2,-(k+p1+p2+p4)^2},{p1^2->0,p2^2->0,p4^2->0,p1 p2->-S/2,p2 p4->-T/2,p1 p4->(S+T)/2,S->3,T->1}],{1,1,1,1},0]
FIESTA, version 1.0.0
UsingC: True
NumberOfLinks: 2
UsingQLink: True
IntegrationCut: 0
IfCut: 0.
Strategy: STRATEGY_S
External integration ready! Use CIntegrate to perform calls
External integration ready! Use CIntegrate to perform calls
QLink created (June 2008 version) ! You can read information on QOpen, QRead, QRemoveDatabase, QRepair, QClose, QList, QSize, QPut, QGet, QSafeGet, QCheck and QRemove
Integration has to be performed up to order 0
Sector decomposition.......0.0312 seconds; 6 sectors.
Variable substitution..........5.1406 seconds.
Decomposing ep-independent term..........0 seconds
Pole resolution..........0.0469 seconds; 20 terms.
Expression construction..........0.0312 seconds.
Replacing variables..........0.0938 seconds.
Epsilon expansion..........0.0469 seconds.
Expanding, making integration string..........0 seconds.
Terms of order -2: 4 (1-fold integrals).
Numerical integration: 4 parts; 2 links;
Integrating......!....3.3593750 seconds; returned answer: 1.333333
Integration of order -2: 1.333333
(1.333333)/ep^2
Terms of order -1: 6 (2-fold integrals).
Numerical integration: 6 parts; 2 links;
Integrating........!..2.1875000 seconds; returned answer: -2.065735 ± 7.*^-6
Integration of order -1: -2.065735 ± 7.*^-6
(1.333333)/ep^2 + (-0.732402 ± 7.*^-6)/ep
Terms of order 0: 14 (2-fold integrals).
Numerical integration: 6 parts; 2 links;
Integrating..........6.0781250 seconds; returned answer: -3.417389 ± 0.000017
Integration of order 0: -3.417389 ± 0.000017
(1.333333)/ep^2 + (-0.732402 ± 7.*^-6)/ep + (-4.386502 ± 0.000018)
DoubleBox:
PrimarySectorCoefficients={4,2,0,0,1,0,0};
SDEvaluate[UF[{k,l},{-k^2,-(k+p1)^2,-(k+p1+p2)^2,-l^2,-(l-k)^2,-(l+p1+p2)^2,-(l+p1+p2+p4)^2},{p1^2->0,p2^2->0,p4^2->0,p1 p2->-S/2,p2 p4->-T/2,p1 p4->(S+T)/2,S->3,T->1}],{1,1,1,1,1,1,1},0]
FIESTA, version 1.0.0
UsingC: True
NumberOfLinks: 2
UsingQLink: True
IntegrationCut: 0
IfCut: 0.
Strategy: STRATEGY_S
External integration ready! Use CIntegrate to perform calls
External integration ready! Use CIntegrate to perform calls
QLink created (June 2008 version) ! You can read information on QOpen, QRead, QRemoveDatabase, QRepair, QClose, QList, QSize, QPut, QGet, QSafeGet, QCheck and QRemove
Integration has to be performed up to order 0
Sector decomposition.....2.2344 seconds; 169 sectors.
Variable substitution..........6.125 seconds.
Decomposing ep-independent term..........0.3594 seconds
Pole resolution..........1.9844 seconds; 1472 terms.
Expression construction..........1.9688 seconds.
Replacing variables..........2.6406 seconds.
Epsilon expansion..........6.6719 seconds.
Expanding, making integration string..........0.0312 seconds.
Terms of order -4: 50 (2-fold integrals).
Numerical integration: 50 parts; 2 links;
Integrating.......!...11.6875000 seconds; returned answer: 0.222223
Integration of order -4: 0.222223
(0.444446)/ep^4
Terms of order -3: 220 (3-fold integrals).
Numerical integration: 106 parts; 2 links;
Integrating....!......80.9218750 seconds; returned answer: -0.849763 ± 5.*^-6
Integration of order -3: -0.849763 ± 5.*^-6
(0.444446)/ep^4 + (-0.366188 ± 0.00001)/ep^3
Terms of order -2: 548 (4-fold integrals).
Numerical integration: 155 parts; 2 links;
Integrating..!........232.6875000 seconds; returned answer: 0.002971 ± 0.000055
Integration of order -2: 0.002971 ± 0.000055
(0.444446)/ep^4 + (-0.366188 ± 0.00001)/ep^3 + (-2.741575 ± 0.000114)/ep^2
Terms of order -1: 901 (5-fold integrals).
Numerical integration: 169 parts; 2 links;
Integrating.!.........861.2031250 seconds; returned answer: 0.755831 ± 0.000298
Integration of order -1: 0.755831 ± 0.000298
(0.444446)/ep^4 + (-0.366188 ± 0.00001)/ep^3 + (-2.741575 ± 0.000114)/ep^2 + (-4.498939 ± 0.000683)/ep
Terms of order 0: 1070 (6-fold integrals).
Numerical integration: 169 parts; 2 links;
Integrating..........2131.9062500 seconds; returned answer: 0.427623 ± 0.001475
Integration of order 0: 0.427623 ± 0.001475
(0.444446)/ep^4 + (-0.366188 ± 0.00001)/ep^3 + (-2.741575 ± 0.000114)/ep^2 + (-4.498939 ± 0.000683)/ep + (-2.924529 ± 0.003499)
An example to make strategy X fail: SDEvaluate[UF[{k,l},{-k^2,-(k-q)^2,-l^2,-(l-q)^2,-(k-l)^2,-v k,-v l},{q^2->-QQ,v^2->vv,q*v->0,QQ->1,vv->1}],{1,1,1,1,1,1,1},0]