etmc/tmLQCD — pull request #165 by urbach (2012-10-10 19:23:32+00:00)
Title: Non-degenerate twisted doublet with clover term

This is providing an implementation for a non-degenerate twisted doublet with clover term. A new input monomial was introduced:

BeginMonomial NDCLOVER
Timescale = 1
StildeMin = 0.01225
StildeMax = 3.5
LocNormConst = 3.3775885577830275786
PrecisionPtilde = 1e-05
DegreeOfMDPolynomial = 48
PrecisionHfinal = 1e-10
ComputeEVFreq = 1
2KappaEpsBar = 0.0935
2Kappamubar = 0.1105
kappa = 0.170
CSW = 1.0
RootsFile = "clover_roots.dat"
EndMonomial

which is very similar to the NDPOLY monomial, but allows for CSW. It is probably not yet the most efficient implementation, however, one gets acceptance. I put this here as a pull request, because now it needs some testing! There is a sample input file with roots in the sample directory.

For instance
* a comparison with the mass degenerate clover monomial would be highly welcome by setting EpsBar = 0 (with high statistics)
* a long statistics comparison with NDPOLY by setting CSW=0 would be useful

Of course I have already checked for one trajectory that in these cases one gets identical results, up to some roundoff, of course.

What I have not yet clarified is whether it is already possible to have two or more polynomial monomials. I tried to unify everything such that this works, however, that was not my main purpose. I would also like to be able to simulate other exponents that 1/2 in order to maybe do the n-th root trick for preconditioning (not like staggered rooting, of course).