Check on Inversion with even-odd-preconditiong (Defect #389)
Using eo-prec. in the inversion requires the reconstruction of the odd solution by means of the even solution.
Apparently, the analytic expression
x_o = M_inv D x_e - M_inv b_o
has to be replaced by
x_o = -(M_inv D x_e + M_inv b_o)
to give exaclty the same results, which has to be!
The change in the code was applied in c9b6ede9.
Why is that?
|blocked by CL2QCD - Defect #387: Broken geometric conventions in dslash_eo?||New||09 Jan 2013|
I recalculated the analytic expression, and it is actually
x_o = - M_inv D x_e + M_inv b_o
However, this does not solve the problem that one has to calculate
x_o = - M_inv D x_e - M_inv b_o
This is true for tm as well as pure wilson fermions.
It seems as if the odd source has to be taken negative, therefore I think this could be related to #387.
- Status changed from New to In Progress
- % Done changed from 0 to 70