now testing if mobility is 0 on the bottom wall

stochasticHydroTools · May 9, 2024 · 79fe6a7 · 79fe6a7
1 parent 3c5002e
commit 79fe6a7
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/tests/test_wall_mobility.py b/tests/test_wall_mobility.py
@@ -62,6 +62,8 @@ def test_self_mobility(Solver, periodicity, ref_file):
     interp_diags = [np.diag(matrix) for matrix in Minterp]
     ref_diags = [np.diag(matrix)[0:3] for matrix in refM] # only take diagonal elements from forces
 
+    assert np.all(np.diag(allM[0]) == [0,0,0]), "Self mobility is not zero on the wall at z=0"
+
     assert np.allclose(interp_diags, ref_diags, atol=3e-2), "Self mobility does not match reference"
 
 if __name__ == "__main__":