diff --git a/src/Dofs/ConstraintHandler.jl b/src/Dofs/ConstraintHandler.jl
index 5f9f321812..47d5b6363b 100644
--- a/src/Dofs/ConstraintHandler.jl
+++ b/src/Dofs/ConstraintHandler.jl
@@ -590,7 +590,7 @@ conditions specified in `ch` such that `K \\ rhs` gives the expected solution.
 !!! note
     `apply!(K, rhs, ch)` essentially calculates
     ```julia
-    rhs[free] = rhs[free] - K[constrained, constrained] * a[constrained]
+    rhs[free] = rhs[free] - K[free, constrained] * a[constrained]
     ```
     where `a[constrained]` are the inhomogeneities. Consequently, the sign of `rhs` matters
     (in contrast with `apply_zero!`).