test: extra: enable valgrind (and coverage) for forcer tests

Add two more easier forcer tests for valgrind, the larger tests take too long

Author: jodavies

check/extra/forcer.frm

@@ -19,7 +19,6 @@ L F =
 B ep;
 P;
 .end
-#pend_if valgrind?
 assert succeeded?
 assert result("F") =~ expr("
        + ep^-3 * ( 1/24 )
@@ -141,6 +140,7 @@ ModuleOption noparallel;
 B ep;
 P;
 .end
+# Too slow, but is supposed to pass.
 #pend_if valgrind?
 assert succeeded?
 assert result("F1") =~ expr("
@@ -269,6 +269,7 @@ ModuleOption noparallel;
 B ep;
 P;
 .end
+# Too slow, but is supposed to pass.
 #pend_if valgrind?
 assert succeeded?
 assert result("F1") =~ expr("
@@ -289,3 +290,219 @@ assert result("F2") =~ expr("
          z3 - 3705/8*z3^2
 ")
 *--#] Forcer_1-expand : 
+*--#[ Forcer_2 :
+* An easier version of the above, for checking with valgrind
+* timeout = 60 seconds.
+#include- forcer.h
+CF f,f1,f2,f3;
+V p2,p3;
+S x3;
+
+* Give 1 or -1. n1 is not used.
+Table randomsign(n1?);
+Fill randomsign() = random_(2)*2-3;
+
+* Zip two functions as:
+*   zip(f1,f2(p1,...,pN),f3(q1,...,qN)) -> f1(p1,q1,...,pN,qN),
+* for N >= 1.
+Table zip(f1?(?a1),f2?(p2?,?a2),f3?(p3?,?a3));
+Fill zip() =
+  + thetap_(nargs_(?a2,?a3)) * zip(f1(?a1,p2,p3),f2(?a2),f3(?a3))
+  + delta_(nargs_(?a2,?a3))  * f1(?a1,p2,p3)
+;
+
+* Element-wise multiplication as:
+*   emul(f1,f2(p1,...,pN),f3(a1,...,aN)) -> f1(p1*a1,...,pN*aN)
+* for N >= 1.
+Table emul(f1?(?a1),f2?(p2?,?a2),f3?(x3?,?a3));
+Fill emul() =
+  + thetap_(nargs_(?a2,?a3)) * emul(f1(?a1,p2*x3),f2(?a2),f3(?a3))
+  + delta_(nargs_(?a2,?a3))  * f1(?a1,p2*x3)
+;
+
+L F1 =
+  #do i=1,2
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,0,0)
+  #enddo
+;
+
+L F2 =
+  #do i=1,2
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,-1,0,0)
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,-1,0)
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,0,-1)
+  #enddo
+;
+
+id Zno1(n1?,...,n14?) =
+  +vx(-Q,p4,p5)
+  *vx(p3,-p4,p10)
+  *vx(p2,-p3,p9)
+  *vx(p1,-p2,p11)
+  *vx(-p5,p6,-p11)
+  *vx(-p6,p7,-p10)
+  *vx(-p7,p8,-p9)
+  *vx(-p1,-p8,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /p2.p4^n12/Q.p2^n13/Q.p3^n14
+;
+
+id Zno2(n1?,...,n14?) =
+  +vx(-Q,p4,p5)
+  *vx(p3,-p4,p11)
+  *vx(p6,p7,p10)
+  *vx(p2,-p3,-p10)
+  *vx(p1,-p2,p9)
+  *vx(-p5,-p6,-p9)
+  *vx(-p7,p8,-p11)
+  *vx(-p1,-p8,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /Q.p2^n12/p1.p4^n13/Q.p3^n14
+;
+
+id Zno3(n1?,...,n14?) =
+  +vx(-Q,p3,p4)
+  *vx(p6,p8,p10)
+  *vx(p5,-p10,p11)
+  *vx(p1,-p3,-p5)
+  *vx(-p4,-p8,p9)
+  *vx(p7,-p9,-p11)
+  *vx(p2,-p6,-p7)
+  *vx(-p1,-p2,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /Q.p6^n12/Q.p8^n13/p3.p6^n14
+;
+
+* Make a random permutation of the loop momenta. The result should be the same.
+multiply f1(p1,...,p11);
+multiply ranperm_(f2,p1,...,p11);
+multiply f3(<randomsign(1)>,...,<randomsign(11)>);
+id f2(?a)*f3(?b) = emul(f2,f2(?a),f3(?b));
+id f1(?a)*f2(?b) = zip(f1,f1(?a),f2(?b));
+id f1(?a) = replace_(?a);
+
+ModuleOption noparallel;
+.sort:input;
+
+#call Forcer(msbarexpand=4)
+B ep;
+P;
+.end
+assert succeeded?
+assert result("F1") =~ expr("
+       + ep^-1 * (  - 15*z5 )
+
+       + 21/2*z7 - 75/2*z6 + 55*z5 - 51*z3^2
+")
+assert result("F2") =~ expr("
+       + 161/2*z7 + 20*z5 + 18*z3^2
+")
+*--#] Forcer_2 : 
+*--#[ Forcer_2-expand :
+* An easier version of the above, for checking with valgrind
+* timeout = 60 seconds.
+#include- forcer.h
+CF f,f1,f2,f3;
+V p2,p3;
+S x3;
+
+* Give 1 or -1. n1 is not used.
+Table randomsign(n1?);
+Fill randomsign() = random_(2)*2-3;
+
+* Zip two functions as:
+*   zip(f1,f2(p1,...,pN),f3(q1,...,qN)) -> f1(p1,q1,...,pN,qN),
+* for N >= 1.
+Table zip(f1?(?a1),f2?(p2?,?a2),f3?(p3?,?a3));
+Fill zip() =
+  + thetap_(nargs_(?a2,?a3)) * zip(f1(?a1,p2,p3),f2(?a2),f3(?a3))
+  + delta_(nargs_(?a2,?a3))  * f1(?a1,p2,p3)
+;
+
+* Element-wise multiplication as:
+*   emul(f1,f2(p1,...,pN),f3(a1,...,aN)) -> f1(p1*a1,...,pN*aN)
+* for N >= 1.
+Table emul(f1?(?a1),f2?(p2?,?a2),f3?(x3?,?a3));
+Fill emul() =
+  + thetap_(nargs_(?a2,?a3)) * emul(f1(?a1,p2*x3),f2(?a2),f3(?a3))
+  + delta_(nargs_(?a2,?a3))  * f1(?a1,p2*x3)
+;
+
+L F1 =
+  #do i=1,2
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,0,0)
+  #enddo
+;
+
+L F2 =
+  #do i=1,2
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,-1,0,0)
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,-1,0)
+    + Zno`i'(1,1,1,1,1,1,1,1,1,1,1,0,0,-1)
+  #enddo
+;
+
+id Zno1(n1?,...,n14?) =
+  +vx(-Q,p4,p5)
+  *vx(p3,-p4,p10)
+  *vx(p2,-p3,p9)
+  *vx(p1,-p2,p11)
+  *vx(-p5,p6,-p11)
+  *vx(-p6,p7,-p10)
+  *vx(-p7,p8,-p9)
+  *vx(-p1,-p8,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /p2.p4^n12/Q.p2^n13/Q.p3^n14
+;
+
+id Zno2(n1?,...,n14?) =
+  +vx(-Q,p4,p5)
+  *vx(p3,-p4,p11)
+  *vx(p6,p7,p10)
+  *vx(p2,-p3,-p10)
+  *vx(p1,-p2,p9)
+  *vx(-p5,-p6,-p9)
+  *vx(-p7,p8,-p11)
+  *vx(-p1,-p8,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /Q.p2^n12/p1.p4^n13/Q.p3^n14
+;
+
+id Zno3(n1?,...,n14?) =
+  +vx(-Q,p3,p4)
+  *vx(p6,p8,p10)
+  *vx(p5,-p10,p11)
+  *vx(p1,-p3,-p5)
+  *vx(-p4,-p8,p9)
+  *vx(p7,-p9,-p11)
+  *vx(p2,-p6,-p7)
+  *vx(-p1,-p2,Q)
+  /<p1.p1^n1>/.../<p11.p11^n11>
+  /Q.p6^n12/Q.p8^n13/p3.p6^n14
+;
+
+* Make a random permutation of the loop momenta. The result should be the same.
+multiply f1(p1,...,p11);
+multiply ranperm_(f2,p1,...,p11);
+multiply f3(<randomsign(1)>,...,<randomsign(11)>);
+id f2(?a)*f3(?b) = emul(f2,f2(?a),f3(?b));
+id f1(?a)*f2(?b) = zip(f1,f1(?a),f2(?b));
+id f1(?a) = replace_(?a);
+
+ModuleOption noparallel;
+.sort:input;
+
+#call Forcer(msbarexpand=4,polyratfunexpand=15)
+B ep;
+P;
+.end
+assert succeeded?
+assert result("F1") =~ expr("
+       + ep^-1 * (  - 15*z5 )
+
+       + 21/2*z7 - 75/2*z6 + 55*z5 - 51*z3^2
+")
+assert result("F2") =~ expr("
+       + 161/2*z7 + 20*z5 + 18*z3^2
+")
+*--#] Forcer_2-expand :