Bartosz Kostrzewa
Bartosz Kostrzewa
> In any case I think we need to investigate the force computation now which will inevitably involve going through the derivation... Well, first I will go through the code...
One thing that perplexes me a little is the lack of any mention of beta in the computation of the force terms. The gauge field is stored independently of beta...
> the force of the smeared part should be independent of the monomials at all, its just the chain rule, isn't it? The output here ``` [DEBUG] Comparison of force...
I changed the code a bit more but I still have my doubts... I first compute the numerical derivative as explained before (the gauge field is the smeared one). Then...
Reducing rho subtanitally seems to indicate that the numerical derivative is wrong after all... actually no.. forget what I just said!! Just had a parameter not upated. For rho 0.0001:...
> the force of the smeared part should be independent of the monomials at all, its just the chain rule, isn't it? Hi, on the train ride to Luxembourg I...
> This is true, but the point is that all of the information has already been encoded within the single input force vector. The modifications to that vector made by...
> [...] and one level of stout smearing doesn't work yet, we can first try to work that single level out. That was the plan but in a sufficiently general...
:) Merry Christmas Carsten and Albert :)
I give up... I've checked the calculation and I think I now match P&M everywhere so...