<< MBv1.2.m

<<AMBRE.m

<<MBnum.m

<<barnesroutines.m

invariants = {p1^2 -> s};

rep=MBrepr[{1},{PR[k1, 0, n1]*PR[k1 + p1, 0,n2]*PR[k2 - k1, 0, n3]*PR[k2, 0, n4]*PR[k4, 0, n5] PR[k2 + k4 + p1, 0, n6]*PR[k3, 0, n7]*PR[k3 - k2 - k4, 0, n8] PR[k3 + p1, 0, n9]},{k4,k3,k2,k1}];

SetOptions[MBnum, Analytical -> True, Numerical -> False, ShowMBrep -> True];

MBanalytic = MBnum[rep, 0, {s -> -1}, {n1 -> 1, n2 -> 1, n3 -> 1, n4 -> 1, n5 -> 1,n6->1,n7->1,n8->1,n9->1}, 4];

Print["before"];

DeleteFile["SE4loop_MBbefore.m"]
Save["SE4loop_MBbefore.m",MBanalytic]

MBintegrate[MBanalytic,{s -> -1}] // AbsoluteTiming

  after = Process[MBanalytic, Range[10]];
  
  Print["after"]
  
  MBintegrate[MBmerge[after],{s -> -1}] // AbsoluteTiming