diff --git a/src/DCPSE_op/DCPSE_op.hpp b/src/DCPSE_op/DCPSE_op.hpp
index abe53792e41c51160f3f055dd15141c0ddd2455f..9ca49ba1b8ab1a55b28e1ca25550725c591f1631 100644
--- a/src/DCPSE_op/DCPSE_op.hpp
+++ b/src/DCPSE_op/DCPSE_op.hpp
@@ -576,9 +576,9 @@ public:
                 auto k = dcp[i].getIndexNN(key,j);
 
 
-                cols[p_map. template getProp<0>(k)*Sys_eqs::nvar + comp] += coeff_dc * coeff / dcp[i].getEpsilonPrefactor(key);
+                cols[p_map. template getProp<0>(k)*Sys_eqs::nvar + comp] += o1.value(k)[i]*coeff_dc * coeff / dcp[i].getEpsilonPrefactor(key);
 
-                cols[p_map. template getProp<0>(key)*Sys_eqs::nvar + comp] += dcp[i].getSign() * coeff_dc * coeff / dcp[i].getEpsilonPrefactor(key);
+                cols[p_map. template getProp<0>(key)*Sys_eqs::nvar + comp] += o1.value(key)[i]*dcp[i].getSign() * coeff_dc * coeff / dcp[i].getEpsilonPrefactor(key);
             }
         }
     }
diff --git a/src/DCPSE_op/DCPSE_op_test2.cpp b/src/DCPSE_op/DCPSE_op_test2.cpp
index 029d1316a4cd981c6fc134b63cce2aebc3fb541b..dd5e6c0a99e88cfcd941262144805ce6fa9b0eaa 100644
--- a/src/DCPSE_op/DCPSE_op_test2.cpp
+++ b/src/DCPSE_op/DCPSE_op_test2.cpp
@@ -196,7 +196,7 @@ BOOST_AUTO_TEST_SUITE(dcpse_op_suite_tests2)
             auto Stokes2 = Adv(V[1],V_star[1])-nu*Lap(V_star[1]);
             Solver.impose(Stokes1,bulk,RHS[0],vx);
             Solver.impose(Stokes2,bulk,RHS[1],vy);
-            Solver.impose(V_star[0], up_p,1,vx);
+            Solver.impose(V_star[0], up_p,1.0,vx);
             Solver.impose(V_star[1], up_p,0,vy);
             Solver.impose(V_star[0], r_p, 0,vx);
             Solver.impose(V_star[1], r_p, 0,vy);