diff --git a/tests/test_conmech/regression/test_density_update.py b/tests/test_conmech/regression/test_density_update.py
index 7b908e16..2e5c6f9b 100644
--- a/tests/test_conmech/regression/test_density_update.py
+++ b/tests/test_conmech/regression/test_density_update.py
@@ -240,7 +240,6 @@ def friction_bound(u_nu: float) -> float:
     return test_suites
 
 
-# TODO find out why "schur" method gives different result in first test scenario
 @pytest.mark.parametrize("setup, density_func, expected_displacement_vector", generate_test_suits())
 def test_direct_solver(solving_method, setup, density_func, expected_displacement_vector):
     runner = NonHomogenousSolver(setup, solving_method)