update

tests/qcd/scalar.py

@@ -119,7 +119,7 @@
 A1 = g.qcd.scalar.action.fourier_mass_term(sqrt_mass)
 A1.assert_gradient_error(rng, U_mom, U_mom, 1e-3, 1e-8)
 r = A1.draw(U_mom, rng)
-assert abs(A1(U_mom) - r) < 1e-10
+assert abs(A1(U_mom)/r - 1) < 1e-10
 
 # group defect would be triggered if sqrt_mass does not have sqrt_mass[k] = sqrt_mass[-k]
 eps = g.group.defect(U_mom[0])
@@ -177,8 +177,6 @@ def __slap(dst, src):
 slap = g.matrix_operator(mat=lap, inv_mat=cg(lap), accept_list=True)
 slap2 = slap * slap
 
-# TODO: need stencil version of __slap
-
 A2 = g.qcd.scalar.action.general_mass_term(M=slap2, sqrt_M=slap)
 
 A2.assert_gradient_error(rng, U_mom, U_mom, 1e-3, 1e-8)