Commit dfe8782c authored by ulrich_y's avatar ulrich_y

Added relevant chen's

parent 6435d95e
......@@ -13,6 +13,8 @@ GetChenGs[]:=Block[Evaluate[ToExpression["F"<>ToString[#]]&/@Range[40]],
}]
]
GetChenFFGs[]:=Uncompress["1:eJzlXduOHDUQ7V4SiY/gha/hhYco4gcC2oiRQlZig0D8POxcti/ucp1T1e61XZOH0WRdl1O3Y7tnkv3x16ePn394GIbnd8Mw/Hx6/vb55S/D8/jy8tNpOL99//Lyz2+/P37drJzf/rv46fXtWf7D09+Pf17XT/+9/Dmd31KSorvbemLk+lNRYbb/y+mPx+eLZYTuu3V0KMTZsBGHHIwN8cJZArtA5EdCWQZPi4sulKAu5UtcgByILm5Cs7m0XOeZ+fDlr+dlm6h5W6spmX2Y3H3/8vLx07fT09dPXy4oTg85vHP0N0CLpRG0dC6Uw2BnWz9fqQTGXBJLzkmt+gF1Db5ONQCpNluI1tjFhucI3zRx8+LKXoK2MwUDsYtSm+4NXcrbYJeTed66WXs2+PxmjqasT2ox9kUvzEOf3vovVsGoakagnKDzJ9nrsoPffCxWykeVoat5dIlELcwlr+sSdUUoG/Dj61t5XK3TD05Y4EiETmWGQ6GDN+QnMosMLX86K4gXDdn2Nirp7CoerJUsIHklq/IuslyTy3iBPbw6ND6TQQ9ZlPVIvNF/RIcXqu+ivGkB3k0hLC6kAytzjTpdPy8hG7NM3sbqUpsKSFSSlxFtAL65Jt2mkwfDXcUTpQW7iWhk9lOENl4upgH8REa2AYqeyORt5Io+5QA5UYp+/qmrgIyi6A0WKBVSrOQyMxIlTGRkG6CEiUzehgoUlDCREW24Ssgoit5gCVOh1ArP22tmTCfea0dEDY7AOxXzKSgHM9+JDK/nk6u0hsWT3ueUIu/NMxF+azwuNDucJjNQy51yKrJRR4RI5dJiKduCt3W9cW7JoOudXki929vm+DFXBuTbQ3DGjuQNKYhRB1GZl5U2XtMHAvlK8xPU3z2QntZOowpXMPcYhgs1XGkpKu86rANL1j9VBqHF6rzQa3Xq0V0AamuFxl6dpyfVN/jG1zCo9NQCsiYT56KTZuE2mWI4/k1Ac6SujRFuaFwbHs3Wx7CRkXuL8Xo/YUmeWiPBi4yUZFGQtjiSglQ/OhSB200tgZBibooUCKJIs1gYRY9bT+rd1jwAUZEozbzjpfmsj1QImfMkldJEjlGyNlKywYv+mkGQIMjNLAg4YRK05pBVBG61DEpCirl1cpAgk8BE8NoxohCT5fOLK8uEInALs5wIKeZGNssMRc+CTAJHbwIJReAWJjARypvjb/Jgl9pjjIGX/0weQKG+BcBrsxkoHO26NQ2Kc997FT14PXNR+Fyhm0QTtONcsVZHH+Q7FPNuqTzz6aMQm9FNljRzetumQsinpx13nMhmddRmhhNZmSe6tT85Zpq+09Bilm4fk4eLN2aR+d2w69iOLl4QIo1EmtUJo+s6VSbDKMTXAMldJJl76HzwBmWeBWmL4Lw+LbvabjZO5WINPy/ov9Hr2u4Yd9xMx3V2KEk3Tu5qwzynEceTVvS4NXL6PmsegApv8ZpWxxdvniSsFV1uPX2wVvS7HR1u16TnUfTgLde3hDUPQF/fJppWx8lOxfsl9gigXa4ghDUPQF9BEk2jY4q6aTCMtXIAHW3kPkH4TdoYi5qAotaKhVtsvsoe1DzWXJPoPwxCTUcSd+YmPZqaHXvHKU292bFjTFJNF+Ri7e8/6yObrrbeqLqnrL/nlfu2mJ7jvaMiH7BvhUvCHbXDzr2/64CrlDkyjYen7OrMFK+iLVJxaNqtT7HUUw3+YYX/Ur7zrs1doQ0346UVBHotBM2hPK+FOHMjMDcyPkfWp6+2hCZyjGubSGkG55wo9tZCyJwvMYQmcowTk0gpBhUoIBJG0+cYlInR3OE4P12KJhg5RtMF2dWDbnMuiLBbKVXFNRe0IRaRcKEUaTDfXhsSoqSgW1eLUKrQNSz9RozqsP4uvzyl9Rxb0OLt5NhwAQcts2Gf6jq4w8sXhU5DUWd11ui7UrUpMQz91aI68Vcaz8vSf6OGLsSUTnKan3Xk+disJ/rib7iVolyoy1nQr3u3ksu/1WizNqgTspUXWgLKZ3zzfZNkeuevY8vj0L7ZAIBvZjAr3tXhI1xAYN6DRLmgkOIRdd108RqskWYyPoOkyZTR6Wo0Y0al7KHRQr0dCo4Lq/8eDNpvDfRWvf8kVtmUK4NqMl3WY0CLSJtMrHYaqY3KkbDqQ9rGQLY5fA0PWv2hOniA5sdp8wUC1dWmA5wtc5JfzxtBhww6DI8hDywlYKSUdyd2DYoYKSF3SiRLAdnMubbWMhE6wFkGcrKuGDHOB6EDnGmI0XxchITDNUSs6wBnGuKRQMyf8yeSEJ057DCg5A8wdAA7qMlryBOKUjkfNeXKcFM3lW6hk3e2iEA0bOZJRgnBUdLK8WSZe3X1jzaiRRWzYG5uihRqzNJS/N5rWEeXrH+qDEKLQSkwMt0FoLa6NDYOznvkpIisl7tj6NYQDt9he9ZU7CtjbbuUuQ15YMmdTSlZ3ekNZoDg6FRgzZwF0PMozBJZAIPggGXOwmIubO6ATR6C28Vh1otlwloNFylS7nkMiF7LWXJFY03pDs7v+cZj/I5196HeUWnLsmKk+O+oCfw7S6+xViluUKKOTMr3RMB3R7ZRibUqibpuXP6LFLwfcdeem9RgGkhKCblTUc8Cmhmh+hi1roTcqahHBjUALJtmJoeL32PJFY2WKaSlOFyYlW1TWTZrQURauKtXJpf9nS6VfzncfVhBS+ZnjkixBi0ux769xnV40QIQZhRyjEqEoUkvAsG9PZn9DzD2eDo="]/.z->y/.xchen->x
GetMuoneGs[]:=DeleteDuplicates[
Cases[{
Import["https://arxiv.org/src/1709.07435v2/anc/arXiv-fam1.m"],
......@@ -79,6 +81,7 @@ MakeFile[getter_, name_, vars_]:=Block[{gs,filename},
Switch[Last[$CommandLine],
"checks/test-chen.f90", MakeFile[GetChenGs , "chen" , {x, y}],
"checks/test-chenff.f90", MakeFile[GetChenFFGs , "chenFF" , {x, y}],
"checks/test-muone.f90", MakeFile[GetMuoneGs , "muone" , {x, y}],
"checks/test-muoneNP.f90", MakeFile[GetMuoneNPGs, "muoneNP", {w, z}],
_, Print["Availble are chen, muone, muonenp"];
......
......@@ -502,6 +502,9 @@ gpl: build/mathlink.o libgpl.a
checks/test-chen.f90: checks/generate.m
@echo "MATH generate \$@"
@\$(MATH) -script \$< \$@
checks/test-chenFF.f90: checks/generate.m
@echo "MATH generate \$@"
@\$(MATH) -script \$< \$@
checks/test-muone.f90: checks/generate.m
@echo "MATH generate \$@"
@\$(MATH) -script \$< \$@
......@@ -516,6 +519,7 @@ echo -n "test: \$(objects) " >> makefile
echo -n "build/test-chenref.o " >> makefile
$HAVE_GINAC && echo -n "build/ginac.o " >> makefile
$HAVE_GINAC && $HAVE_MCC && echo -n "build/test-chen.o " >> makefile
$HAVE_GINAC && $HAVE_MCC && echo -n "build/test-chenFF.o " >> makefile
$HAVE_GINAC && $HAVE_MCC && echo -n "build/test-muone.o " >> makefile
$HAVE_GINAC && $HAVE_MCC && echo -n "build/test-muoneNP.o " >> makefile
......
......@@ -295,12 +295,14 @@ CONTAINS
subroutine do_timing_tests(n)
use gtestchen , only: inichen =>args
use gtestchenff , only: inichenff =>args
use gtestmuone , only: inimuone =>args
use gtestmuonenp, only: inimuonenp=>args
implicit none
integer, intent(in) :: n
integer i
complex(kind=prec) :: cargs( 1399,5,n)
complex(kind=prec) :: fargs( 540,5,n)
complex(kind=prec) :: pargs( 198,5,n)
complex(kind=prec) :: nargs( 1733,5,n)
real(kind=prec) :: z, x, y, w
......@@ -310,7 +312,8 @@ CONTAINS
do i=1,n
z = ran2(ranseed) / 2.
x = ran2(ranseed)*(1-z) + z
cargs(:,:,i) = inichen(cmplx(x), cmplx(z))
cargs(:,:,i) = inichen (cmplx(x), cmplx(z))
fargs(:,:,i) = inichenff(cmplx(x), cmplx(z))
w = ran2(ranseed) ! 0<w<1
z = ran2(ranseed) * (sqrt(1-w+w**2)-sqrt(w)) + sqrt(w)
......@@ -322,9 +325,12 @@ CONTAINS
enddo
cargs(1181,:,:)=cargs(1181,:,:)/0
fargs(367,:,:)=fargs(367,:,:)/0
open(unit=9, file="stats.txt")
write(9,*) "Chen form factor"
call do_one_speed_test(fargs,9,"Chen FF")
write(9,*) "Chen"
call do_one_speed_test(cargs,9,"Chen")
write(9,*) "MUonE-planar"
......
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