Merged fermionmatrix kernels including dslash broken? (Unit Test #323)
With cdf48d884ad7aa9d3bd068382ee0ec00af45b296, the merged fermion kernels including dslash seem to be broken.
One can see that in the failing test for the dslash + m_tm_inverse_sitediagonal(_minus) kernels. If one sets mu=0, the first two of these should give the same result as a pure dslash!
The other kernel, dslash + gamma5, passes the tests, however, it seems strange that both tests give the same result, where in the one the argument to the dslash is even and odd in the other. Thus, they should use different gaugelinks and give different results!
|precedes CL2QCD - Defect #416: Merged spinor kernels broken||Feedback||27 Jan 2013|
|follows CL2QCD - Feature #412: Make own test module for merged fermion kernel||Done||24 Jan 2013|
The fermionmatrix gamma_5*Aee (and similar Aee_minus) has (in the pure Wilson case) a structure like
Hence, merging gamma5 with dslash would result in
- dslash + gamma5 merged
- gamma5 * 1
which would yield no benefit as gamma5 is still called alone!
I removed the dslash+gamma5 kernel.
Now remaining are the merged kernels of the dslash with the diagonal matrices in the twisted mass case.
Actually, these test fail (at least partly) on the CPU! Hence, there seems to be a more fundamental problem!
- Assignee set to Christopher Pinke