diff --git a/src/common_operators/homogeneousdata_operator.jl b/src/common_operators/homogeneousdata_operator.jl
index e26124ec..f8c5113b 100644
--- a/src/common_operators/homogeneousdata_operator.jl
+++ b/src/common_operators/homogeneousdata_operator.jl
@@ -174,7 +174,7 @@ function apply_penalties!(A, b, sol, O::HomogeneousData{UT}, SC::SolverConfigura
             penalty = O.parameters[:penalty]
             AE = A.entries
             for dof in bdofs
-                AE[dof, dof] = penalty
+                rawupdateindex!(AE, (a, b) -> b, penalty, dof, dof, 1)
             end
             flush!(AE)
         end
diff --git a/src/common_operators/interpolateboundarydata_operator.jl b/src/common_operators/interpolateboundarydata_operator.jl
index 23e4f0e5..df3cee71 100644
--- a/src/common_operators/interpolateboundarydata_operator.jl
+++ b/src/common_operators/interpolateboundarydata_operator.jl
@@ -204,7 +204,7 @@ function apply_penalties!(A, b, sol, O::InterpolateBoundaryData{UT}, SC::SolverC
         if assemble_matrix
             AE = A.entries
             for dof in bdofs
-                AE[dof, dof] = penalty
+                rawupdateindex!(AE, (a, b) -> b, penalty, dof, dof, 1)
             end
             flush!(AE)
         end