The interpolation bounce back boundary condition was wrongly implemented because the equilibrium calculation was done by looking at the 'deviation_only' property of the lb method equilibrium. However, here we need to apply the equilibrium on the PDF values and not on moments etc.