Preamble

<<"~/Desktop/azurite/code/Azurite.wl";
TemporaryDirectory="/Users/alessandro/tmp/";
SingularDirectory="/Applications/Singular.app/Contents/bin/";

Example

Massless Pentagon Box

LoopMomenta={l1,l2};
ExternalMomenta={k1,k2,k3,k4};
Propagators={l1^2,(l1-k1)^2,(l1-k1-k2)^2,(l1-k1-k2-k3)^2,(l2+k1+k2+k3)^2,(l2+k1+k2+k3+k4)^2,l2^2,(l1+l2)^2,(l1+k4)^2/2,(l2+k1)^2/2,(l2+k2)^2/2};
Kinematics={k1^2->0,k2^2->0,k3^2->0,k4^2->0,k1 k2->s12/2,k1 k3->s13/2,k1 k4->s14/2,k2 k3->s23/2,k2 k4->s24/2,k3 k4->(-s12-s13-s14-s23-s24)/2};
Numerics={s12->1,s13->10037,s14->1201,s23->7477,s24->9311};
Preparation[];

FeynmanGraph[Range[8]]

MIs=FindAllMIs[Range[8]]

FeynmanGraph[#[[1]]]&/@MIs

AzuriteForm2Int/@MIs

Graphs

Ind1=MIs[[1,1]]

Ind2=MIs[[-1,1]]

FeynmanGraph@Ind1

FeynmanGraph@Ind2

FeynmanGraph/@CompleteGlobalSymmetryRep[Ind1]

2^8-Length[%]

FeynmanGraph/@CompleteGlobalSymmetryRep[Ind2]

PropagatorSymmetry[Ind1]

MomentaSymmetry[Ind1]

PropagatorSymmetry[Ind2]

MomentaSymmetry[Ind2]

IBPs on maximal cut

Massless Triple Box

LoopMomenta={l1,l2,l3};
ExternalMomenta={k1,k2,k4};
Propagators={l1^2,(l1-k1)^2,(l1-k1-k2)^2,l3^2,(-l3-k1-k2)^2,(l1+l3)^2,(l2-l3)^2,l2^2,(l2-k4)^2,(l2+k1+k2)^2,
(l1+k4)^2,(l2+k1)^2,(l3+k1)^2,(l3+k4)^2,(l1+l2)^2};
Kinematics={k1^2->0,k2^2->0,k4^2->0,k1 k2->s/2,k2 k4->(-s-t)/2,k1 k4->t/2}; 
Numerics={s->1,t->-6};
Preparation[];

{MIsMc,IBPsMC}=IntegralRed[Range[10]];

MIsMc

Int[1,1,1,1,1,1,1,1,1,1,-2,0,0,0,0]/.IBPsMC

{MIsMc,IBPsMC}=IntegralRed[Range[10],HighestPower->6,WorkingPower->6];

Int[1,1,1,1,1,1,1,1,1,1,-2,0,0,0,0]/.IBPsMC

Int[1,1,1,1,1,1,1,1,1,1,-6,0,0,0,0]/.IBPsMC

Int[1,1,1,1,1,1,1,1,1,1,-1,-1,-1,-1,-2]/.IBPsMC